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Chaptar 1. Haterottructuraa for High Parformanca Davicaa 


Chapter 1. Heterostructures for High Performance 
Devices 

Academic and Research Staff 
Professor Clifton G. Fonstad. Jr. 

Visiting Scientists and Research Affiliates 

Anton Failia,' Sheila Prasad,^ Jae-Jin Lee,’ Jong Tae Park* 

Graduate Students 

Rajni J. Aggarwal, Thomas P.E. Broekaert, Geoffrey F. Burns, Woo-Young Choi, Isako Hoshino, Paul S. 
Martin, Lung-Han Peng, Yakov Royter, Krishna V. Shenoy, Richard A. Singer, Jurgen H. Smet, James C. 
VIcek 

Technical and Support Staff 

Kelley S. Donovan, Angela R. Odoardi, Richard R. Perilli 


1.1 Introduction 

The broad objective of our research effort is to 
develop lll-V quantum heterostructures for high 
performance electronic, optoelectronic, and pho¬ 
tonic devices for applications in high speed optical 
communications and signal processing. To this 
end, we are developing: (1) new, higher perfor¬ 
mance materials systems including InP-based 
InGaAIAs heterostructures and <111> oriented 
strained layer superlattices; (2) novel approaches 
to integrate laser diodes on VLSI-level electronic 
integrated circuits; (3) a new family of quantum- 
well-base, tunnel-barrier n-n-n transistors and 
near- and far-infrared optoelectronic devices; and 
(4) new damage-free in situ processing techniques 
for fabricating advanced quantum structure and 
embedded heterostructures. 

The following sections describe our progress 
during the past year in the above research areas. 
Our group works closely with Professors Hermann 
A. Haus, James G. Fujimoto, and Erich P. Ippen to 
develop the optical device application, character¬ 
ization, and modeling aspects of this program, and 
with Professor Sylvia T. Ceyer to develop new In 
situ processing techniques. 


1.2 Computer Controlled 
Growth of Lett ice-Matched 
InGaAIAs Heterostructures on 
InP 

Sponsors 

Charles S. Draper Laboratories 
Contract DL-H-418483 
DARPA/NCIPT 

Joint Sen/ices Electronics Program 
Contract DAAL03-89-C-0001 
Contract DAAL03-92-C-0001 

Project Staff 

James C. VIcek, Professor Clifton G. Fonstad, Jr. 

In electronic and optical semiconductor devices, 
the need for both graded-composition and hyper- 
abrupt metallurgical junctions frequently arises. 
The use of graded-composition junctions allows 
for precise control of the confinement of charge 
carriers and optical fields through spatially varying 
bandgaps and refractive indices, respectively. On 
the other hand, hyper-abrupt junctions may be 


^ Centro Studi e Laboratori Telecomunicazoni, Torino, Italy. 

2 Northeastern University, Boston, Massachusetts. 

’ Electronics and Telecommunications Research Institute (ETRI), Daejeon, Korea. 
* Incheon Unversity of Korea, Incheon, Korea. 
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da»ired in some devices, most notably quantum 
effect devices, where a sharp interface is desired to 
reduce fluctuations in the confinement energy 
which might arise from graded interfaces due to 
shutter transients. 

(n this work, we have grown the quaternary 
InxGayAli-yAsi-x, lattice-matched to InP substrates. 
The need to lattice-match, such as to an InP sub¬ 
strate, further constrains the ratio of constituent 
fluxes. 

We have implemented a computer-automated 
MBE control system whose aim is to provide 
precise control of the constituent and dopant 
fluxes such as is necessary to achieve graded- 
composition lattice-matched alloys and maintain 
uniform lattice-matched compositions in the pres¬ 
ence of shutter events. Key to this system are 
accurate models of the characteristics of the 
effusion cells, including not only the static flux 
versus temperature relationship, but also the time- 
dependent behavior of the cell in the presence of 
changes in either setpoint and/or shutter status. 
These models approximate the temporal response 
of the effusion cells to a setpoint change with a 
single-pole system function, with the location of 
the pole being determined by a direct measure¬ 
ment of the flux profile generated by a step change 
in the cell setpoint temperature. 

In the case of graded-composition layers, the 
thermal lag time of the cell requires modification of 
the time-dependent effusion cell setpoint temper¬ 
ature profile, most notably to eliminate the over- 
and undershoot effects at the endpoints of the 
graded layer which would result if the effects of 
cell response time were neglected. With such 
modifications, we have achieved linearly graded 
flux profiles whose deviation from the desired 
profile can be held within one percent. 

To remove flux transients which arise from the 
operation of the cell shutters, we have structured 
the control system so that it seeks to maintain the 
flux in the effusion cell at a constant value, irre¬ 
spective of the shutter status. In practice, this 
entails maintaining the setpoint temperature of the 
cell at a lower value with the shutter closed than 
with it open and performing an exponential ramp 
sequence upon change of shutter status. Using 
this technique, we have reduced cell shutter tran¬ 
sients from as high as 30 percent to one percent or 
less. 


1.3 Molecular Beam Epitaxy of 
InGaAIAs Strained-Layer 
Heterostructures on 111 GaAs 
and InP 

Sponsors 

DARPA/NCIPT 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
Contract DAAL03-92-C-0001 

Project Staff 

Richard A. Singer, Professor Clifton G. Fonstad, Jr. 

The issues associated with the molecular beam 
epitaxial (MBE) deposition of inx(GayAh_y)i_xAs 
upon (111)B InP have been extensively explored. 
As a result, a reproducible growth algorithm has 
been developed whose execution has resulted in 
bulk n and p-type InGaAs and InGaAIAs of com¬ 
parable crystallinity to simultaneously grown (100) 
oriented epilayers. 

When epitaxial layers are pseudomorphically 
deposited upon a (111) zincblende substrate, a 
piezoelectric field is generated due to the ionic 
character of the chemical bond between the con¬ 
stituent elements and the lack of inversion sym¬ 
metry in the (111) orientation. Recently, this 
strain-generated electric field has been demon¬ 
strated in a GaAs based p-i-n modulator whose 
absorption edge was in the vicinity of 0.9 /im. 

Growth on (111) InP is especially difficult because 
of the volatility of the InP surface and the prob¬ 
lems attendant to (111) surface chemistry. Not 
surprisingly, initial attempts to grow epitaxial layer 
of InGa^ on (111)B InP according to the proven 
methods of (100) orientations were unsuccessful. 
Moreover, reflection high-energy electron diffrac¬ 
tion (RHEED) patterns, which usually provide 
essential in situ information about the condition of 
a growing surface, are unclear in this material 
system. ^ a result, the most critical step in the 
MBE growth process—the synchronization of the 
desorption of the substrate's native oxide layer 
with the initiation of growth—is performed ‘in the 
dark.' 

To address these issues, an alternative growth 
algorithm, called quasi-migration enhanced 
epitaxy, was developed and named for the unusu¬ 
ally low intensity of As overpressure, low substrate 
temperature, and slow growth rate which comprise 
its parameter set. It was learned that an intense 
group III overpressure, expressed in a rapid rate of 
growth, is neede d to stabilize the (111) growth 
front; a y/l9 x .y/l9 surface reconstruction that is 
rich in group III atoms. This is quite different from 
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the case of (100) materials whose c(2x8) surface 
reconstructions are As stabilized. Another impor¬ 
tant modification was the elimination of the tradi¬ 
tional 5: 1: 1 H2SO4H2O2-H20 etching step in the 
pre-growth substrate preparation procedure. As a 
result, the thickness of the residual oxide layer was 
reduced which, in turn, allowed growth to be suc¬ 
cessfully initiated without the need of a risky high 
temperature desorption cycle. 


1.4 Monolithic Fabrication of 
Strain-free GaAIAs Laser 
Diodes on Silicon Substrates 

Sponsors 

DAR PA/NCI PT 

IBM Corporation Fellowship 

Project Staff 

Geoffrey F. Burns, Anton Failla, Professor Clifton 
G. Fonstad, Jr. 

Prospects for monolithic integration of lll-V elec¬ 
trical and optical devices with Si circuits have 
fueled vigorous research in direct heteroepitaxial 
growth of these compounds on Si substrates. 
Developing lll-V optical sources useful for VLSI 
optical interconnects is one goal which has been 
pursued by several groups. Along these lines we 
have focused upon the laser diode due to its 
capacity for high speed modulation and efficient 
electrical to optical power conversion. 

In this work molecular beam epitaxy of (AI,Ga)As 
heterostructures directly on silicon substrates has 
been investigated for use in fabricating semicon¬ 
ductor lasers for monolithic optoelectronic inte¬ 
gration. A base laser process was developed and 
evaluated using homoepitaxy of (AI,Ga)As on 
GaAs substrates, incorporating a graded-index 
separate confinement heterostructure (GRINSCH) 
quantum well active region. Etched ridge wave¬ 
guide lasers were fabricated from these layers, 
achieving continuous wave threshold currents of 
19 mA and differential quantum efficiency of 66%. 
This same laser process was adopted to 
GaAs-on-Si growth initiation layers to produce 
monolithic (AI,Ga)As laser on silicon substrates. 
High dislocation density and residual thermal 
tension in the laser epitaxial layer degraded 
threshold currents and slightly degraded external 
quantum efficiency, with threshold currents as low 


as 40 mA and differential efficiencies as high as 60 
percent. In addition, the same problems precluded 
continuous wave operation for these devices. 

To address residual thermal tension, a technique 
for monolithic epitaxial separation and reattach¬ 
ment was developed to produce strain-free 
GaAs-on-Si. Rapid thermal annealing of these 
layers was also investigated for alleviating defect 
density imposed by the heteroepitaxial growth 
process. Photoluminescence studies demonstrated 
that epitaxial separation, along with rapid thermal 
annealing 600*C, produced strain-free material. In 
addition, these layers could be annealed at temper¬ 
atures up to 800‘C without reintroducing thermal 
stress and demonstrated further material improve¬ 
ment as indicated by photoluminescence intensity. 
Epitaxial separation and annealing was also 
applied to (AI,Ga)As laser epitaxy layers to 
produce the first strain-free lasers fabricated 
monolithically on Si substrates, an important 
achievement for advancing heteroepitaxy as a 
viable optoelectronic integration strategy. 

1.5 Integration of Vertical 
Cavity Surface Emitting Lasers 
on GaAs Integrated Circuits 

Sponsors 

OARPA/NCIPT 
IBM Corporation Fellowship 
National Science Foundation Fellowship 
Vitesse Semiconductor 

Project Staff 

Krishna V. Shenoy, Geoffrey F. Burns, Professor 
Clifton G. Fonstad, Jr. in collaboration with J. 
Mikkelson* 

There exists a long history of interest in the inte¬ 
gration of high performance optoelectronic 
devices, particularly lasers and detectors, on high 
density electronic integrated circuits, and consider¬ 
able attention has been focused on GaAs-on-Si 
processes because silicon has traditionally been 
the only source of VLSI circuitry. However, with 
the recent availability of GaAs-based VLSI, realized 
using refractory-metal gate processes, we feel that 
this picture has changed significantly. It is the 
purpose of this program to demonstrate that the 
growth of optoelectronic heterostructures on par¬ 
tially processed GaAs-VLSI wafers is a viable alter¬ 
native to GaAs-on-Si for applications requiring the 


s Vitesse Semiconductor, Camarillo, California. 
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Other applications for these devices include pure 
FM laser diodes with reduced AM noise, two 
wavelength integrated but independently con¬ 
trolled laser diodes and other non-laser devices 
like tunable narrow bandwidth filters and light 
modulators. 

1.9 Design and Fabrication of 
Distributed Feedback (DFB) 
InGaAIAs Laser Diodes Grown 
by Molecular Beam Epitaxy 

Sponsors 

DARPA/NCIPT 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
Contract DAAL03-92-C-0001 

Project Staff 

Woo-Young Choi, Professor Clifton G. Fonstad, Jr. 

It is essential that laser diodes for optical fiber 
communication applications operate with a single 
oscillation frequency. Typically, schemes of distrib¬ 
uted feedback (DFB) or distributed Bragg reflector 
(OBR) are used in which gratings above or below 
the active layer (DFB) or outside the active layer 
(DBR) perform the act of frequency selection. To 
make such a device structure, either epitaxial 
growth is initiated on a corrugated substrate, or 
regrowth is performed on the epitaxial material 
onto which gratings have been formed. 

In this project, we are in the pursuit of a DFB 
device structure in which gratings are formed 
entirely after the complete epitaxial growth. Our 
motivation for this approach stems from the fact 
that we are working on laser diodes grown by the 
conventional solid-source Molecular Beam Epitaxy 
technology and, due to its growth kinetics, suc¬ 
cessful MBE growth on the corrugated surface is 
very difficult to achieve. 

We have proposed a ridge stripe structure in which 
gratings are made on the side walls of the ridge as 
well as on the bottom channels next to the ridge. 
According to our initial calculation, there is 
enough coupling between the optical wave and 
gratings to result in single mode selection. 

To realize such a device structure, we are collab¬ 
orating with Professor Henry I. Smith's group at 
MIT and utilizing their x-ray lithography tech¬ 
nology. Currently, we are in the process of com¬ 


bining x-ray lithography technology with the 
ridge-stripe laser diode fabrication technology. 

1.10 Laser Diode Modeling and 
Design for Narrow Linewidth 
Operation 

Sponsor 

Charles S. Draper Laboratories 
Project Staff 

Yakov Royter, Professor Clifton G. Fonstad, Jr. in 
collaboration with J.H. Hopps’ 

Many of the current laser applications require 
narrow linewidth output characteristics. Unfortu¬ 
nately, semiconductor lasers, possessing such 
advantages as small size, direct modulation capa¬ 
bility, and integrability with other optical and elec¬ 
tronic devices, have broad linewidths in 
comparison to most other kinds of lasers. The rel¬ 
atively broad linewidths are due to the low Q of 
the cavity and to the intrinsic effects specific to the 
semiconductor gain medium. Since most of the 
linewidth reduction methods have involved the 
increase of the cavity Q, our goal is to investigate 
and to help in the reduction of the intrinsic line- 
width broadening effects. Reduction of these 
effects, in conjunction with the increase of the 
cavity Q, would produce semiconductor lasers 
with narrowest possible linewidths. 

The first step in our research is to investigate the¬ 
oretically the intrinsic linewidth broadening effects. 
In particular, we concentrate on the linewidth 
broadening parameter, a, the largest intrinsic line- 
width broadening factor, which describes the cou¬ 
pling of gain and refractive index fluctuations in 
the gain medium via carrier density fluctuations. 
Also, due to their apparent promise over conven¬ 
tional diode lasers, we consider only quantum-well 
semiconductor lasers. We begin the investigation 
of the a parameter by first calculating the carrier 
energy bands of the active region semiconductor 
material, using both parabolic band and k • p 
approximations. From the energy bands, we can 
obtain gain and refractive index profiles and the 
linewidth broadening parameter. Thus, the energy 
band calculations are the first major task in our 
investigation. 

To test our models, we plan to compare the calcu¬ 
lated laser characteristics, such as threshold 
current, gain, refractive index and linewidth, with 


f Charles S. Draper Laboratories, Cambridge, Massachusetts. 
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be the superior method for optimization because of 
its speed and convenience. 

Small-signal S parameter measurements were 
made by on-wafer probing of the devices using 
the Hewlett-Packard 851 OB automatic network 
analyzer and the Cascade Microtech microwave 
probe station for 39 different bias points. The bias 
dependence of elements of the small-signal equiv¬ 
alent circuit was determined and it was found that 
five of the elements of the circuit were highly bias 
dependent. These results are being used for the 
large-signal modeling of the HBT which is in pro¬ 
gress. Measurements are also in progress to deter¬ 
mine the third-order intermodulation product 
(IPS). The optimized model for high IPS which is 
related to device linearity is to be determined. 

1.13 Analysis of Three-Terminal 
n-n-n Quantum Well Base, 
Tunnel-Barrier Transistors 

Sponsor 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
Contract DAAL03-92-C-0001 

Project Staff 

Thomas P.E. Broekaert. Professor Clifton G. 
Fonstad, Jr. 

Three-terminal n-n-n quantum-well base, tunnel 
barrier transistors, consisting of a resonant 
tunneling structure with an ohmic contact to the 
quantum well have been fabricated. Although an 
ohmic contact to the quantum well was success¬ 
fully achieved by heavily doping the quantum well, 
it was found that all such structures exhibit strong 
leakage currents from base to emitter. 

The conductance of surface quantum wells in such 
transistor structures has been measured success¬ 
fully, and the dependence of sheet carrier concen¬ 
tration, sheet conductance, and carrier mobility as 
a function of temperature has been measured. No 
transistor action was seen for the tested temper¬ 
ature range of 10 to 300 K. The lack of transistor 
action is tentatively attributed to fast hot carrier 
relaxation, resulting in large base currents and 
small collector currents. Possible mechanisms for 
the hot carrier relaxation are impurity scattering, 
polar optical phonon scattering, and electron- 
electron scattering. 

Several non-idealities that hamper the active 
behavior of the resonant tunneling transistor have 
been analyzed. First, electron-polar optical phonon 
and electron-electron scattering from second to 


first subband in the quantum well are believed to 
be the dominant mechanisms at low temperatures 
that limit the ultimate gain of the resonant 
tunneling transistor. Second, due to the fact that 
this device only has a single dipole layer, it follows 
that the emitter accumulation region is not com¬ 
pletely shielded from the collector potential, 
which results in a low output impedance and a 
low voltage gain. 

A novel device structure is introduced, in which 
the hot carrier relaxation, which is prevalent in 
InGaAIAs resonant tunneling structures, is 
reduced. This reduction in hot carrier relaxation 
rate is obtained by having a quantum well in 
which the quantum-well level for base potential 
modulation and the quantum-well level for emitter 
to collector tunneling are located at different 
points in the Brillouin zone, e.g., direct vs. indirect 
valley quantum well. A possible material system for 
this novel device structure would be 
GaAs/AIAs/Ge/AIAs/GaAs for emitter/bar- 
rier/base/barrier/collector respectively. This mate¬ 
rial system would also improve the performance of 
the bipolar quantum-well resonant tunneling tran¬ 
sistor and the hot electron transistor. 

1.14 Applications of AlAs 
Etch-Stop Layers in 
InGaAIAs/lnP Heterostructure 
Electronics and Optoelectronics 

Sponsor 

Joint Senrices Electronics Program 
Contract DAAL03-89-C-0001 
Contract DAAL03-92-C-0001 

Project Staff 

Thomas P.E. Broekaert, Woo-Young Choi, Pro¬ 
fessor Clifton G. Fonstad, Jr. 

We have used wet chemical etching solutions that 
allow the selective etching of InP lattice-matched 
InGaAIAs quaternary compounds using thin 
pseudomorphic AlAs layers as etch stops to 
explore new device concepts. 

We developed the etchants used here last year, 
and they typically consist of succinic acid (or other 
dicarboxylic acids), ammonia, and hydrogen 
peroxide. The etchant is well buffered and can be 
used over a wide pH range including from 4.2 to 
7.0 by varying the amount of ammonia added. In 
addition, the etchant is compatible with Cr/Au 
contact metallization and standard positive 
photoresists. Typically, the InGaAs etch rate is 
about 100 nm/min. At a pH of 4.2, the etch rate 
of InGaAs is found to be over 1000 times the etch 
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rate of AlAs, while the etch rate of IrtAIAs is over 
500 times that of the AlAs. A typical AlAs stop 
layer is about 10 monolayers (m.l.) thick (2.73 
nm). Buffered HF can be used to remove the AlAs 
stop layer without etching InGaAIAs to any signif¬ 
icant degree. 

These selective etchants have enabled us, for the 
first time, to measure the conductance and 
mobility of a directly contacted quantum well in a 
resonant tunneling structure, and they open up the 
possibility of fabricating a whole new set of novel 
quantum devices. 

Also, we utilized the selective etchants in fabri¬ 
cating ridge-stripe InGaAIAs/lnP GRIN-SCH 
MOW laser diodes. By inserting a thin 
pseudomorphic AlAs layer at a desired location 
during the epitaxial growth and then selectively 
etching away InGaAIAs materials up to the etch 
stop layer, we were able to fabricate laser devices 
with ridges whose heights are precisely and reli¬ 
ably controlled. 

1.15 Electrical Transport 
Studies in Directly Contacted 
InGaAs Quantum Wells 

Sponsors 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
Contract DAAL03-92-C-0001 
National Science Foundation Fellowship 
U.S. Army Research Office 

Project Staff 

Rajni J. Aggarwal, Thomas P.E. Broekaert, Pro¬ 
fessor Clifton G. Fonstad, Jr. 

Using a recently developed selective etch, we have 
been able to demonstrate direct electrical contact 
to the quantum well of an AlAs/lnAs/InGaAs reso¬ 
nant tunneling structure. This contact serves as 
the base contact in our three-terminal, unipolar 
resonant tunneling transistor (RTT). Our devices 
incorporate high indium content materials in the 
quantum well. In addition to their favorable trans¬ 
port properties, these compounds have low 
Schottky barrier heights, enhancing our ability to 
make ohmic contact to them. We are limited in 
the indium fraction that we can use by the lattice 
mismatch of such compounds with the InP based 
lattice constant of the device. 

We believe that the performance of our device is 
limited by two factors: (1) doping within the 
quantum well and (2) depletion of the well due to 
surface exposure. Neither of these problems has 


been overcome by the inclusion of high indium 
content materials alone. We feel that it is neces¬ 
sary to study the nature of the contact itself in 
order to engineer around the present limitations. 
Specifically, we are interested in elucidating the 
effect, if any, of quantization on Schottky barrier 
height. In addition, we would like to exper¬ 
imentally determine any potential barriers that may 
exist in electron transport from quasi-confined 
(contact) regions to confined (intrinsic base) 
regions within the device. 

We are measuring the Schottky barrier height of a 
variety of quantum well contacts and comparing 
them to those obtained for comparable bulk mate¬ 
rials. Among the experimental variables are 
quantum well and barrier thickness and well 
indium content. In addition, we are performing 
photoluminescence studies in hopes of exper¬ 
imentally observing the shift in energy of the n=1 
quantum state within the well as one barrier is 
removed. 

1.16 Investigation of Infrared 
Intersubband Emission from 
InGaAs/AIAs/lnP Quantum Well 
Heterostructures 

Sponsor 

National Science Foundation Fellowship 

Project Staff 

Jurgen H. Smet, Professor Clifton G. Fonstad, Jr. 
in collaboration with Professor Qing Hu 

Subbands in quasi-two-dimensional semicon¬ 
ductor structures exhibit a large electric dipole 
matrix element, making them attractive as the 
active medium in an infrared source. Triple 
quantum-well structures are proposed in which 
electrons are selectively injected into the second 
subband of a wide center quantum well. After 
relaxation to the first subband, electrons are selec¬ 
tively removed. Two important issues need to be 
addressed: the achievement of population inver¬ 
sion and the confinement of the emitted radiation. 

To avoid fast relaxation through 
electron/LO-phonon scattering, it is desirable to 
design a structure for a lasing frequency smaller 
than the LO-phonon frequency. This requires a 
subband separation of less than 30 meV. This 
requirement fixes the center well width and the 
width of the side wells. The balancing of gain and 
losses determines the required population inver¬ 
sion. The threshold current density is then, in 
essence, this population inversion divided by the 
intersubband relaxation time, provided that popu- 
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Chapter 2. High-Frequency InAIAs/InGaAs 
Metal-lnsulator-Doped Semiconductor Field-Effect 
Transistors (MIDFETs) for Telecommunications 

Academic and Research Staff 
Professor JesOs A. del Alamo 

Visiting Scientists and Research Affiliates 
Dr. Yuji Awano' 

Graduate Students 
Sandeep R. Bahl, Brian R. Bennett 

Undergraduate Students 
Michael H. Leary, Akbar A. Mooiji 

Technical and Support Staff 
Kelley S. Donovan, Angela R. Odoardi 


2.1 Introduction 

Sponsors 

Charles S. Draper Laboratories, Inc. 

Contract DL-H-418488 
Fujitsu Laboratories 
Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
Contract DAAL03-92-C-0001 
Texas Instruments 

The goal of this project is to investigate 
InAIAs/n-^-InGaAs Metal-Insulator Doped channel 
Field-Effect Transistors (MIDFETs) on InP. These 
devices are of great interest for applications in 
long-wavelength lightwave communication 
systems and ultra-high frequency high-power 
microwave telecommunications. 

InAIAs/InGaAs Modulation-Doped Field-Effect 
Transistors (MODFETs) on InP have recently 
emerged as an optimum choice for a variety of 
microwave and photonics applications. This is 


because the outstanding transport properties of 
InGaAs have yielded devices with very low-noise 
and high-frequency characteristics. Unfortunately, 
the low breakdown voltage of InAIAs/InGaAs 
MODFETs on InP (typically less than 5 V) severely 
restricts their use in high-power applications, such 
as large-signal microwave amplification and laser 
driving. It also forces the use of a separate high 
voltage supply to operate the Metal- 
Semiconductor-Metal (MSM) photodetectors in 
InP photonics receivers. 

A device strategy with great potential for power 
handling is the InAIAs/n"^-InGaAs MIDFET fea¬ 
turing an undoped insulator and a thin, heavily- 
doped channel. In this structure, the breakdown 
voltage, Vb, is large and can be engineered using 
pseudomorphic insulators^ and channel quantiza¬ 
tion (shown below). Drain current, lo, can also be 
considerably improved with InAs-rich channels.’ 
The attainment of high power, however, demands 
a large Id x Vb product. InAs-rich channel devices 
unfortunately suffer from a low breakdown voltage 
due to: (1) the reduced channel bandgap, and (2) 


1 Fujitsu Laboratories, Atsugi, Japan. 

3 S.R. Bahl, W.J. Azzam, and J.A. del Alamo, 'Strained-Insulator ln(Ali-.As/n'*^-lno 63 Gao 47 As Heterostructure Field- 
Effect Transistors,* IEEE Tnns. Electron Dev. 38(9): 1986-1992 (1991). 

) S.R. Bahl and J.A. del Alamo, 'An InosjAlo^sAs/n'^-ln.Gsi-xAs Heterostructure Field-Effect Transistor with an 
In-Enriched Channel,' Proceedings of the Second Internetionel Conference on InP end Releted Compounds, 
Denver, Colorado, April 23-25,1990, p. 100. 
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severe gate leakage at the sidewall of the mesa, 
where the gate comes in contact with the heavily- 
doped channel * 

Towards solving this severe problem in this period 
of performance, we have studied the effect of 
quantizing the InGaAs channel. This has resulted 
in a drastically improved breakdown voltage 
through a quantum mechanically engineered 
enlargement of the effective bandgap of the 
channel. We have also solved the isolation 
problem through a selective chemical recessing of 
the edge of the channel on the mesa sidewall. We 
have finally integrated all our understanding devel¬ 
oped over the last three years into a device with a 
strained (InAs-rich) channel, strained (AlAs-rich) 
insulator, quantized channel, optimized channel 
doping, and edge isolated device which has dis¬ 
played unprecedented power and frequency per¬ 
formance. A detailed description of these 
experiments is presented in this report. 


2.2 Quantum-channel 
InAIAs/m-InGaAs MIDFETs 

In previous work, we have shown that enriching 
the InAs mole fraction of lno. 63 Gao. 47 As channel 
results in MIOFETs with superior transport proper¬ 
ties. However this comes at the cost of a severely 
reduced breakdown voltage. Vs, presumably 
through the decrease in the energy gap, Eg.^ A 
method of increasing the effective energy gap in 
the channel is to introduce energy quantization by 
reducing the channel thickness to dimensions 
comparable to the electron wavelength (figure 1 ). 
In fact, it has been shown that in 
lnoB3Gao.47As/lno62Alo4BAs quantum wells the pho¬ 
toluminescence emission wavelength decreases” 
with a reduction in well thickness. 



equilibrium of InAIAs/n'*'-InGaAs HFETs for thick and 
thin channels, showing the increase in effective energy 
gap in the thinner channel. 

In this work, we exploit this effect to enhance the 
breakdown voltage of 

lno62Alo4aAs/n*-lnoB3Gao47As MIDFETs on InP. 
We have doubled Vb by shrinking the lnn. 63 Gao. 47 As 
channel thickness from 350A to lOOA, keeping 
other physical parameters constant. The principle 
behind our work should allow one to better exploit 
the excellent transport properties of InAs-rich 
InGaAs and other promising narrow gap semicon¬ 
ductors like InAs” and InSb. 

A cross section of the device structure is shown in 
figure 2. The wafers were grown on SI-InP by 
MBE in MIT's Riber 2300 system. In an effort to 
keep the channel charge constant, the thickness of 
its undoped portion was varied while the thickness 
of its heavily doped portion was kept constant at 
IOOA. Four wafers were subseouently grown with 
subchannel thicknesses of 250A, IOOA, 50A, and 
OA, i.e., total channel thicknesses of 350A, 200A, 
150A, and IOOA. Devices were fabricated with 
nominal gate lengths of 1 pm and widths of 30 


” S.R. Bahl and J.A. del Alamo, "An lno. 62 Alo. 4 sAs/n'^-ln.Gai-,As Heterostructure Field-Effect Transistor with an 
In-Enriched Channel,’ Proceedings of the Second tnternationel Conference on InP and Related Compounds, 
Denver, Colorado, April 23-2S, 1990, p. 100. 

” D.F. Welch, G.W. Wicks, and L.F. Eastman, ’Optical Properties of GalnAs/AllnAs Single Quantum Wells,’ Appl. 
Phys. Lett. 43(8): 762-764 (1983); W. Stolz, K. Fujiwara, L. Tapfer, H. Oppolzer, and K. Ploog, ’Luminescence of 
lnoB 3 Gso 47 As/lnos 2 Alo 48 As Quantum Well Heterostructures Grown by Molecular Beam Epitaxy,’ Inst. Phys. Conf. 
Ser. 74(3); 139-144 (1988). 

* C.C. Eugster, T.P. Broekaert, J.A. del Alamo, and C.G. Fonstad, 'An InAIAs/lnAs MODFET,’ IEEE Electron Dev. 
Lett 12(12); 707-709 (1992). 
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Figure 2. Cross-section of grown device structures in 
quantum-channel experiment. 


/im. Fabrication is similar to that used in Bahi, 
A 2 zam, and del Alamo (1991).^ 

Our baseline device, 200A channel thickness, had 
a peak transconductance, gmdMrt). of 202 mS/mm 
and a maximum drain current, ld(m«). of 312 
mA/mm. The output conductance, ga. was 5.73 
mS/mm, resulting in a voltage gain. Ay, of 35. 
These are average values over five devices. gm(pMk) 
and ld(m») were measured at Vd. => 4 V, and gd at 
Vds» 4 V and V„ = 0 V. The contact and channel 
sheet resistances, measured by TLM, are 0.37 
Q • mm and 625 O/D, respectively. 

The reverse gate breakdown voltage, Vb, was 
measured with the source and drain grounded, and 
was defined at a reverse gate current of 500 ftA, 
corresponding to about 5% of the peak drain 
current of our baseline device.* Here we focus on 
our main result: the increased breakdown voltage 
in devices with thinner channels and the exper- 



Figura 3. Breakdown voltage, Vb, for typical HFETs 
with channel thicknesses of lOOA, 150A, 200A, and 
350A. 


imental confirmation of the energy quantization 
therein. 

Figure 3 shows typical reverse gate l-V character¬ 
istics of HFETs as a function of channel thickness. 
As shown, Vb increases gradually from 9 V at a 
channel thickness of 350A to 10.6 V at 200A, 11 .9 
V at 150A, and to 19.1 V at lOOA. Average Vb 
measurements over several devices are within 1 V 
of the typical values shown in figure 3. Vb for the 
200A channel is also similar to what we have pre¬ 
viously measured in identical devices grown and 
processed separately.* The drastic improvement in 
breakdown voltage of our quantized-channel 
HFETs is 8 significant merit for high-power appli¬ 
cations. This is particularly so in this material 
system because typical InAIAs/InGaAs MODFET 
breakdown voltages are on the order of 5 V.* 

In order to verify the bandgap enhancement in the 
channel as a result of carrier quantization, we have 
carried out photoluminescence (PL) measurements 
on unprocessed portions of the device samples at 
77 K. The results are shown in figure 4. The 
energy of the peak PL intensity increases from 
0.83 eV for the 350A channel, to 0.86 for 200A, to 
0.88 for 150A, and to 0.92 eV for the lOOA 
channel. The literature reports a temperature- 
independent PL energy shift of about 60 meV over 


7 S.R. Bahl, W.J. Azzam, and J.A. del Alamo, 'Strained-Insulator lnxAli-xAs/n'^-lno. 53 Gao. 47 As Heterostructure Field- 
Effect Transistors,' /£££ Trans. £lactron Dev. 38(9): 1986-1992 (1991). 

• S.R. Baht and J.A. del Alamo, 'An lnos 2 Alo. 48 As/n'^-lnxG 8 t-xAs Heterostructure Field-Effect Transistor with an 
In-Enriched Channel," Proceedings of the Second International Conference on InP and Related Compounds. 
Denver, Colorado, April 23-25,1990, p. 100. 

• P.C. Chao, A.J. Tessmer, K.-H.G. Duh, P. Ho, M.-Y. Kao, P.M. Smith, J.M. Ballingall, S.-M. Liu, and A.A. Jabra, 
'W-band Low-Noise InAIAs/InGaAs Lattice-Matched HEMT’s," l£££ electron Dev. Lett. 11(1): 59-62 (1990); 
Y.-C. Pao, C.K. Nishimoto, R. Majidi-Ahy, J. Archer, N.G. Bechtel, and J.S. Harris, Jr., 'Characterization of 
Surface-Undoped lno 52 Alo 4 BAs/lno 53 Gao 47 As/lnP High Electron Mobility Transistors," l£££ Trans, electron. Dev. 
37(10): 2165-2170 (1990). 
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Wavelength (^m) 


Figure 4. Photoluminescence spectra of the device 
heterostructures, showing an increase in photolum¬ 
inescence energy with decreasing channel thickness. 

the lno 53 Gao 47 As bulk value for 100A thick 
lno52Alo48As/lno53Gao47As quantum wells.^® This 
value of 60 meV is also consistent with simple cal¬ 
culations for a finite square well of thickness 
lOOA. Our slightly larger shift of 90 meV in the 
lOOA channel might be due to depletion at the top 
and possibly bottom interfaces, causing a reduced 
effective well thickness. Additionally, band 
bending in the insulator results in an effective 
stronger potential than a square well (figure 1). 
Both effects will tend to enhance the strength of 
carrier quantization over a simple quantum well. 

Electron quantization, as schematically shown in 
figure 1, also implies a reduction of the effective 


conduction band discontinuity between channel 
and metal gate. This, in return, is expected to 
result in enhanced forward gate leakage current. 
We have experimentally found this to be the case 
at 300 K and 77 providing us in this manner 
with an independent confirmation of the presence 
of quantization in the channel. 

Unfortunately, we have found that a reduction in 
the channel thickness results in the degradation of 
transconductance and peak drain current.^' gtn(pMk) 
decreases from 262 mS/mm to 138 mS/mm, and 
Id(m 4 ») from 451 mA/mm to 208 mA/mm in going 
from a channel thickness of 350A to 100A. 
However, the output conductance g^, improves 
due to the enhanced channel aspect ratio, 
decreasing from 10.8 mS/mm to 1.84 mS/mm. 
This results in an enhancement in the voltage gain, 
Av. from 24 to 75. The degradation in gni(p.ak) and 
ld<muc) results from an increased source resistance 
R., a reduced channel sheet cha je concentration 
Ot, and degraded mobility, fi. From 350A to 100 
A, Rg increases from 1.4 £J • mm to 2.7 n • mm, n, 
drops from 2.38 x 10’* cm-* to 

1.77 X 10’* cm-*, and/I decreases from 
4318 cm*/V-s to 3591 cm*/V-s. n, and (i were 
measured by the Hall-effect. The acknowledged 
poor quality of the reverse InGaAs/lnAIAs interface 
at the back of the channel can be held responsible 
for the mobility reduction.’* As the undoped sub¬ 
channel is thinned down, the reverse interface has 
a larger impact on carrier mobility since the 
channel electrons travel closer to the reverse 
InAIAs/InGaAs interface. There are, however, 
techniques that could mitigate this degradation: 
superlattice buffers to improve mobility’* and the 
Migration-Enhanced Epitaxy (MEE) growth tech¬ 
nique to reduce interface roughness.’* 


^0 D.F. Welch, G.W. Wicks, and LF. Eastman, 'Optical Properties of GalnAs/AllnAs Single Quantum Wells,* Appl. 
Phys. Lett. 43(*.1): 762-764 (1983); W. Stolz, K. Fujiwara, L. Tapfer, H. Oppolzer, and K. Ploog, 'Luminescence of 
lno 53 Gao 47 As/lros 2 Alo 48 As Quantum Well Heterostructures Grown by Molecular Beam Epitaxy,' Inst. Phys. Conf. 
Ser. 74(3): 139-144 (1985). 

’■I S.R. Bahl and J.A. del Alamo, 'A Quantized-channel fno 62 Alp 48 As/n'^-lnos 3 Gao. 47 As HFET with High Breakdown 
Voltage,' Extended Abstracts of the Materials Research Society 1990 Fall Meeting, Boston, Massachusetts, no. 
EA-21, p. 117. 


A.S. Brown, J.A. Henige, and M.J. Delaney, 'Photoluminescence Broadening Mechanisms in High Quality 
GalnAs-AllnAs Quantum Well Structures,' Appl. Phys. Lett. 52(14); 1142-1143 (1988); T. Sajoto, M. Santos, 
J.J. Heremans, M. Shayegan, M. Heiblum, M.V. Weckwerth, and U. Meirav, 'Use of Superlattices to Realize 
Inverted GaAs/AIGaAs Heterojunctions with Low-temperature Mobility of 2x10® cm*/V s,' Appl. Phys. Lett. 
54(9); 840-842 (1989). 

T. Sajoto, M. Santos, J.J. Heremans, M. Shayegan, M. Heiblum, M.V. Weckwerth, and U. Meirav, 'Use of Super¬ 
lattices to Realize Inverted GaAs/AIGaAs Heterojunctions with Low-temperature Mobility of 2 xIO® cm*/V s,' 
Appl. Phys. Lett. 54(9); 840-842 (1989). 

y.C. Chen, P. Bhattacharya, and J. Singh, 'Strained Layer Epitaxy of InGaAs (on GaAs) by MBE and Migration 
Enhanced Epitaxy - Comparison of Growth Modes and Sudace Quality,' J. Cryst. Growth 111(1-4): 228-232 
(1991). 
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Thinning down the subchannel reduces the sheet 
carrier concentration, possibly from backside 
depletion. A reduced channel doping can also 
result in an improvement of breakdown voltage. In 
a separate experiment on similar devices (with 200 
A channel thickness), we examined the impact of 
channel doping on Vb- This experiment indicated 
that a reduction in sheet charge concentration 
from 2.38x 10’* cm“*to 1.77 x 10’* cm”* should 
result in an improvement of Vb by 5 V. Our exper¬ 
imental observation of a 10 V improvement in 
going from a 350A to lOOA channel is evidence 
that quantization is instrumental in drastically 
enhancing the breakdown characteristics of our 
HFETs. 

In conclusion, lno.s2Alo48As/n'*'-lnos3Gao47As 
HFETs have been fabricated with channel thick¬ 
nesses of lOOA, 150A, 200A, and 350A. For 
devices with channel thicknesses of lOOA, the 
breakdown voltage improved twofold over the 350 
A devices. This is postulated to arise partially fronr 
an enlargement of the effective energy gap caused 
by energy quantization introduced from electron 
confinement, as observed by PL. 

2.2.1 Elimination of Mesa-Sidewall 
Gate-Leakage Current by Selective 
Sidewall Recessing 

We have shown in previous work’® that fabrication 
of lnAIAs/n*-lnGaAs MIDFETs by conventional 
mesa isolation results in sidewalls where the 
InGaAs channel is exposed and comes in contact 
with the gate metallization running up the mesa 
(figure 5). The low Schottky barrier height of 
metals on InGaAs’® potentially results in a sidewall 
leakage path from the gate to the channel. 
Sidewall-leakage in InAIAs/InGaAs HFETs results 


in excessive gate-leakage current,’® reduced break¬ 
down voltage,’® and increased sidegating.’* In 
HFETs, sidewall-leakage was also found, by the 
present authors, to worsen with high doping, an 
increased channel thickness, and increased x (x > 
0.53) in the InxGai-xAs channel.’® This latest trend 
is particularly consequential, since HFETs with 
InAs enriched channels have shown excellent 
electron transport properties.’* 

Airbridging and ion implantation have been uti¬ 
lized in InGaAs HFET isolation. Airbridging is 
complex, and implantation demands capital- 
intensive tools. In this paper, we propose and 
demonstrate a new, simpler technique by which 
we selectively recess the exposed InGaAs channel 
into the sidewall. The subsequently e-beamed 
gate metallization does not enter this cavity and 
remains isolated from the channel edge. This one- 
step technique is self aligned to the mesa, requires 
no additional masks, gate-length insensitive, and 
works for mesas in all crystallographic directions 
on the ( 100 ) surface. 

The MBE grown heterostructure, lattice-matched 
to S.l. InP, consists of (bottom to top) a lOOOA 
undoped lno 52 Alo. 48 As buffer layer, an lnos3Gao47As 
channel consisting of a lOOA undoped subchannel 
and a 80A heavily Si doped (No = 6x10’® cm"®) 
transport layer, a 300A undoped lnos 2 Alo 48 As gate 
insulator layer, and an undoped 50A lnos 3 Gao 47 As 
cap. For clarity in SEM imaging, an additional 
heterostructure with a thicker lno. 536 ao 47 As 
channel (120A doped and 300A undoped) was 
processed along with the device sample. 

Devices were fabricated by first chemically etching 
a mesa down to the InP substrate using a 
H 2 S 04 :H 202 ;H 20 1:10:220 etch. Then, before 
removing the mesa-level photoresist mask, the 
wafer was dipped for 45 seconds into a SA:H202 


16 S.R. Bahl and J.A. del Alamo, 'An lno 62 Alo 48 As/n+-lnxGai_*As Heterostructure Field-Effect Transistor with an 
In-Enriched Channel," Proceedings of the Second International Conference on InP and Related Compounds, 
Denver, Colorado, April 23-25, 1990, p. 100. 

16 H.H. Wieder, 'Fermi Level and Surface Barrier of Gaxlni_xAs Alloys,' Appl. Phys. Lett. 3B(3): 170-171 (1981). 

1* Y-J. Chan, D. Pavlidis, and G-l. Ng, 'The Influence of Gate-Feeder/Mesa-Edge Contacting on Sidegating Effects 
in lno 52 Alo 48 As/n+-lno 53 Gao 47 As Heterostructure FET's," l££E Electron Dev. Lett. 12 (7): 360-362 (1991). 

1® S.R. Bahl and J.A. del Alamo, "An lno 52 Alo 48 As/n+-lnxGai_xAs Heterostructure Field-Effect Transistor with an 
In-Enriched Channel," Proceedings of the Second International Conference on InP and Related Compounds. 
Denver, Colorado, April 23-25, 1990, p. 100; S.R. Bahl, M.H. Leary, and J.A. del Alamo, 'Mesa-Sidewall Gate- 
Leakage in InAIAs/InGaAs Heterostructure Field-Effect Transistors," submitted for publication. 

19 S.R. Bahl and J.A. del Alamo, 'An lno 62 Alo 48 As/n*-lnxG 8 i_xAs Heterostructure Field-Effect Transistor with an 
In-Enriched Channel," Proceedings of the Second International Conference on InP and Related Compounds. 
Denver, Colorado, April 23-25, 1990, p. 100; U.K. Mishra, A.S. Brown, and S.E. Rosenbaum, "DC and RF Perfor¬ 
mance of 0.1 ^m Gate Length Alo 48 lno 62 As-Gao 38 lno 62 As Pseudomorphic HEMT's," Proceedings of the Interna¬ 
tional Electron Device Meeting, 1988. pp. 180-183. 
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Figure 5. Three-dimensional perspective of an HFET 
showing the mesa-sidewall gate-leakage path. 


6:1 solution^o to selectively etch the exposed 
portion of the lno. 53 Gao. 47 As channel in a self- 
aligned manner. The SA solution was prepared by 
adding 1 liter H 2 O to 200 g. succinic acid with the 
addition of ammonium hydroxide until the pH was 
5.5. A planar selectivity of 23:1 was measured, 
with the lno 52 Alo 4 BAs etching at 25A/min. For ref¬ 
erence, a portion of the wafer was masked during 
this etch by coating with photoresist using a 


paintbrush. In this portion, mesa etching was per¬ 
formed, but no selective sidewall-recessing was 
carried out. For the ohmic contacts, 2000A of 
AuGe and 6 OOA of Ni were evaporated, lifted off, 
and RTA alloyed at 360° C. For the gate and pad, 
300A of Ti, 300A of Pt, and 2000 A of Au were 
e-beam evaporated and lifted off. This process 
should be compatible with recessed-gate devices. 
If the same etchant is used for gate-recessing, then 
sidewall-recessing could be performed simultane¬ 
ously. 

To study sidewall-leakage, we fabricated special- 
purpose heterojunction diodes with an active 
Schottky area of 10,000 /im^. Gate-metal/mesa- 
sidewall overlaps were created by etching grooves 
through the active diode during mesa formation, 
and then depositing gate metal on top. The ohmic 
contact surrounds the gate region. Twelve diodes 
were fabricated in each die, with sidewall-overlap 
lengths, U. of 0, 200, 400, and 600 am, running in 
each of [ 011 ], [ 001 ] and [ 011 ] crystallo¬ 

graphic directions. Figure 6 shows a photograph 
of the diode test structure with various sidewall- 
overlap lengths along one selected crystallographic 
direction. 



Figure 6. Photograph of sidewall-leakage test diodes with 10,000 /rm^ area and sidewall overlap lengths of (left to 
right) 600, 400, 200, and 0 ^m. 


*0 T.P.E. Broekaert and C.G. Fonstad, 'AlAs Etch-Stop Layers for InGaAsAs/lnP Heterostructure Devices and Cir¬ 
cuits," /fff Trans. Electron Dev., forthcoming (1992). 
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Figure 7. SEM photograph and explanatory sketch of the heterostructure sidewall, showing the isolation of the 
gate metal from the channel. 


Figure 7 shows an SEM photograph and a sketch 
of the finished heterostructure with 420 
A ino.B3Gao.47As channel processed with the other 
device samples. The presence of a sidewall cavity 
is clearly revealed. This confirms the successful 
action of the selective etchant. The overhang 
formed by the InAIAs insulator prevents the gate 
from contacting the channel edge. 

Figure 8 shows the l-V characteristics of typical 
heterojunction diodes from the non-sidewall- 
recessed portion of the wafer with grooves along 
the [011] direction. Both forward and reverse 
diode currents increase with sidewall-overlap 
length, as expected from the increase in contact 
area. This unmistakably proves the existence of a 
significant sidewall-leakage path. The figure 
shows that sidewall-leakage plays a major role 
from forward-bias till threshold ( - 0.82 V), when 
the channel is depleted. We also find (not shown) 
that sidewall-leakagejlepends on crystallographic 
orientation as l[011] > l[001] > l[011]. For 
typical diodes with 200 pm, I increases by 
38% from [011 ] to [011 ] at V = -2 V. 

In the sidewall-recessed portion of the wafer, both 
sidewall-length and orientation dependence of the 
diode current have disappeared. This is shown in 
figure 9, a plot of the l-V characteristics of twelve 


diodes (one of each stripe and orientation) from 
the sidewall isolated portion of the wafer. All 
reverse characteristics are tightly clustered and 
randomly distributed around the baseline area 
leakage characteristic observed for zero overlap 
length in figure 8. The absence of orientation and 
overlap length dependence confirms that sidewall- 
leakage has been completely eliminated for all 
crystallographic orientations. 



• 1.5 -w -oa 0 oa IP 

VOLTAGE (V) 


Figure 8. l-V characteristics of test diodes without 
sidewall recessing,_as a function of sidewall overlap 
length along the [011] direction on the wafer surface. 
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Figure 9. I-V characteristics of test diodes with 
sidewall recessing, for different orientations and 
sidewall overlap lengths. 



GATE-VOLTAGE (V) 


Figure 10. Gate-diode characteristics of HFETs, with 
and without sideyyall isolation, for « 1 and mesa 
sidewall along [011 ]. 


On the same wafer, we have also fabricated HFETs 
with nominal gate length, Lg, and width, Wg, of 1 
nm and 30 pm respectively. Figure 10 shows the 
impact of sidewall isolation on the gate character¬ 
istics. Sidewall isolation results in lower gate- 
leakage current, a larger gate turn-on voltage, and 
an increased breakdown voltage. The complete 
elimination of sidewall-leakage in HFETs was con¬ 
firmed by the disappearance in orientation depend¬ 
ence of the subthreshold current (not shown) at 
several gate-lengths. 

in conclusion, a simple self-aligned technique for 
eliminating mesa-sidewall gate-leakage has been 
developed. This technique uses selective etching 
to etch the exposed part of the InGaAs channel 
into the mesa-sidewall, creating a cavity to isolate 
the channel from the gate. Measurements on spe¬ 
cially designed diodes have shown complete elimi¬ 
nation of sidewall-leakage. This process should be 
useful for any kind of heterostructure in the 
InAlAs/InGaAs system. 


2.2.2 Doubly Strained 
InAIAs/n^-InGaAs HFET 

In this section, we pull together all the knowledge 
developed over the last three years on 
InAIAs/n-^-InGaAs MIDFETs and present an ori¬ 
ginal device design that features for the first time: 
( 1 ) complete elimination of sidewall leakage, ( 2 ) 
an optimized channel doping level, and (3) high 
electron confinement by using a strained 
(AlAs-rich) insulator and a quantized and strained 
(InAs-rich) channel. The resulting device, though 
moderate in L«, displays unprecedented power 
handling capabilities. 

The MBE grown HFET (figure 11) on S.l. InP 
consists of (bottom to top), a 1000A lno. 52 Alo. 4 «As 
buffer, a 75A lno. 53 Gao. 47 As subchannel, a 100 
A n'^-InxGai-xAs channel, a 300A lno. 41 Alo. 69 As 
strained insulator, and a 50A lno 53 Gao. 47 As cap. 
The insulator thickness and mole fraction were 
chosen to avoid misfit dislocations while providing 
maximum conduction band discontinuity. The 
subchannel thickness was chosen to. introduce 
electron quantization and increase Eg in the 
channel without significantly degrading current 
driving capability. The impact of channel doping 
was examined by growing three wafers with lattice 
matched (x = 0.53) channels nominally doped to 
4, 6 , and 8 x 10’* cm-*. From power consider¬ 
ations, as indicated below, an optimum doping 
level of 6 X 10’* cm-* was selected. This is the 
doping at which the InAs-rich channel (x=0.65) 
device was grown. 
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Figure 11. Cross-section of doubiy-strained device 
structure. 


Devices (see sketch in figure 11) were processed 
by mesa etching, sidewall isolation etching. RTA 
alloying of lifted-off Ge/Au/Ni ohmic contacts, 
and lifting-off Ti/Pt/Au gates and pads. Sidewall 
isolation was done by selectively etching the 
InGaAs channel edge into the sidewall to create a 
cavity. Gate-metal running up the mesa sidewall 
does not enter the cavity, and therefore does not 
contact the channel edge. Measurements are 
reported for 1.9 fim gate length (Lg) HFETs. All 
DC measurements are averaged over several 
devices. 

In our x=0.53 channel devices, higher channel 
doping results in higher lo(mix), but lower 
Vb lD{mM) for dopings of 4, 6, and 8 x 10’* cm-* 
is respectively 117, 220, and 371 mA/mm. Vs is 
21.1, 15.6, and 5.3 V, giving an lo(mn) x Vb 
product of 2.5, 3.4, and 2.0 W/mm respectively, f, 
and fmtx also increase with doping, ft is 9.2, 11.3, 
and 12.8 GHz, and fmtx is 34, 42, and 68 GHz, 
respectively. 



Figure 12. t-V characteristics of x=0.6S HFET. 


For our doubly-strained InAs-rich (x=0.65) 
channel device, we found 

OmtPMk) = 201 mS/mm, and Id(imi) = 276 mA/mm 
(figure 12). ft was 14.9 GHz and fm„ 101 GHz 
(figure 13). This gives an effective electron 
velocity of 1.8 x 10^ cm/sec for our 1.9 fim Lg, 
comparable to values obtained in MODFETs of 
similar Lg ‘s. The high imu is a consequence of 
low gd (0.9 mS/mm at peak ft) as this device is 
now able to enter accumulation at Vg, ^ 0 V. The 
voltage gain at peak ft is 150. The average Vb was 
12.8 V (figure 14). This results in an average 
Id X Vb product of 3.5 W/mm. Devices capable of 
handling as much as Id x Vb = 4 W/mm were 
obtained. These values represent an improvement 
over similar gate length MODFETs of 2 to 3 times. 
High electron confinement and reduced gate 
leakage in these HFETs are instrumental in 
obtaining these values. 

In conclusion, in the pursuit of high-power 
InAIAs/InGaAs HFETs, we have combined for the 
first time in a single device, an optimized heavily- 
doped channel, an AlAs enriched insulator, an 
InAs enriched channel, channel quantization, and 
mesa-sidewall leakage elimination. We have 
thereby been able to fabricate a high-voltage, 
high-current device with excellent microwave 
characteristics. Scaled down versions of this 
device are very promising for high-power micro- 
wave and photonics applications. 
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Figure 13. H 21 and MAG vs. frequency for x=0.6S. 



Figure 14. Reverse Ig vs. Vg, characteristics for typical 
HFETs showing that our improvements have mitigated 
the degradation in Ve for x=0.65 (Wg = 30pm). 


2.2.3 Publications and Meeting 
Papers 

Bahl, S.R.. and J.A. del Alamo. 

“InAIAs/n+’InGaAs HFETs: Impact of Dislo¬ 
cations and Channel Quantization.' Paper pre¬ 
sented at the Workshop on Compound 
Semiconductors Materials and Devices 
(WOCSEMMAD). Ft. Lauderdale, Florida, Feb¬ 
ruary 18-20,1991. 

Bahl, S.R., W.J. Azzam, and J.A. del Alamo. 'Ori¬ 
entation Dependence of Mismatched 
ln,Ali_*As/lno 5368047 AS HFETs.' J. Crystal 
Growth 111 (1 -4): 479-483 (1991). 

Bahl, S.R., W.J. Azzam, and J.A. del Alamo. 
‘ Strai ned -1 nsu lator I n^Ali -.As/n-^-1 no.e 3 Gao .47 As 
Heterostructure Field-Effect Transistors.' IEEE 
Trans. Electron Dev. 38(9): 1986-1992 (1991). 

Bahl, S.R., and J.A. del Alamo. 'Elimination of 
Mesa-Sidewall Gate Leakage in InAIAs/InGaAs 
HFETs by Selective Sidewall Recessing.' Pro¬ 
ceedings of the 18th International Symposium 
on Gallium Arsenide and Belated Compounds, 
Seattle, Washington, September 9-12,1991. 
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3.1 Project Description 

Sponsor 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
Contract DAAL03-92<C-0001 

Project Staff 

Professor John M. Graybeal, Professor Henry I. 
Smith. Or. Bernard S. Meyerson, George E. 
Rittenhouse 

In this program we seek to examine the supercon¬ 
ducting Josephson behavior in short-channel elec¬ 
tronic structures where interference effects can 
take place. Additionally, our hybrid superconduc¬ 
ting-semiconducting structure represents the first 
attempt at Josephson coupling via resonant 
tunneling. A significant consequence of this 
approach is that the energy scale for modulating 
the Josephson coupling is set by the quantum 
confinement levels and not the intrinsic semicon¬ 
ducting gap voltage. 

Our device uses a thin vertical Si membrane 
(< lOOOA) as a degenerate semi-conducting weak 
link between the superconducting source and 
drain electrodes which lie on either side. This 
structure is in many ways an electrical analogue of 
a Fabry-Perot interferometer. The membrane 
defines a vertical quantum well in which the 
charge carriers can interfere with each other (either 
constructively or destructively), thus creating a 
resonant tunneling structure between the source 


and drain for both the Cooper pairs (i.e., the 
Josephson channel) and the normal quasiparticle 
(i.e., electron) channel. A separately applied gate 
voltage biases the levels within the semicon¬ 
ducting well, and ideally will modulate the 
Josephson current by varying the wavelength of 
the tunneling particles. Previous efforts on such 
hybrid Josephson devices have not utilized either 
resonant tunneling or quantum confinement and 
their lack of gain arises from the significant mis¬ 
match between the semiconducting and supercon¬ 
ducting gap energy scales.’ A key motivation for 
our devices is that the energy scale is set by 
quantum confinement and can be tuned sepa¬ 
rately from either of the gap energies. For Si well 
widths in the range of 500-1OOOA, the energy 
levels separation in the well can be more compa¬ 
rable to the superconducting gap energy, i.e., 
many times smaller than the semiconducting gap. 

Recent progress has now validated our basic 
device geometry in uniformly doped low-mobility 
Si:P. We are now progressing to fabricate the 
devices on doped high-mobility epitaxial layers 
grown on Si wafers at IBM. Higher mobilities will 
increase the characteristic lengths for the carriers 
in the semiconducting well, thereby better ena¬ 
bling internal interferences. 

A schematic of the device geometry is shown in 
figure 1. More details of the fabrication process 
and progress will be published elsewhere.* The 
primary aspects of the fabrication are described 
here. The originating wafer is <110> Si. 
phosphorous-do|^ via ion implantation to yield 


1 IBM Corporation, Thomas J. Watson Research Center, Yorktown Heights, New York. 

2 A.W. Kleinsasser, T.N. Jackson, 6.D. Pettit, H. Schmid, J.M. Woodall, and D.P. Kern, IEEE Tnn$. M»g. 
MAG-25(2): 1274 (1989), and references within. 

* K. Early, G.E. Rittenhouse, J.M. Graybeal, and Henry I. Smith, *Sub-100 and sub 10-nm-thick Membranes Aniso- 
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Figure 1. Schematic of the planar device structure. 
The vertical membrane consists of two etch pits sepa¬ 
rated by a length L ^ 10OOA. The typical height for the 
vertical membrane is %3000A. 

1x10'’cm~3 carriers at low temperature. A 
300A SIsNa layer is used as an etchant mask, with 
the membrane defined by aljgning a mask stripe 
along either the [112] or [112] crystalline axis. 
The membrane fabrication is done via wet etching, 
using a potassium hydroxide/isopropyl alcohol 
solution which is well known for its highly aniso¬ 
tropic etching of Si. Because of the wafers' crys¬ 
tallographic orientation, the resultant terminating 
(111) planes are perpendicular to the wafer 
surface and form the walls of the vertical mem¬ 
brane structure. Such processing yields 
atomically-smooth surfaces which possess a sig¬ 
nificantly higher degree of parallelism and a signif¬ 
icantly lower surface defect density than could be 
obtained by reactive-ion-etching (RIE) or other 
energetic processes. The niobium counterelec¬ 
trodes are subsequently deposited, and source- 


drain isolation is obtained by planarization 
followed by a chlorine RIE etch. 

The electrical results of such devices fall into two 
categories. Those with wall thicknesses (denoted 
L in figure 1) greater than roughly lOOOA dis¬ 
played no Josephson currents, and showed 
current-voltage (l-V) characteristics consistent 
with superconducting-insulator-normal (SIN) 
metal behavior. The thin 'insulating' layer is pre¬ 
sumably a Schottky barrier at the Nb/Si interface, 
whereas the lack of a Josephson current reflects 
the fact that the device thickness L is significantly 
greater than the proximity-induced supercon¬ 
ducting coherence length {n in the Si 
({n« 65-100A at 4.2K for these low-mobility 
devices). However, devices with membrane thick¬ 
nesses comparable to 600A displayed finite 
Josephson currents, and fine l-V behavior con¬ 
sistent with SNS behavior. The l-Vs are shown as 
a function of temperature for one device in figure 
2. The product of the critical current and the 
high-bias device resistance URn was as high as 0.4 
meV at low temperatures, consistent with the best 
previous results on Nb/Si/Nb devices. Careful 
examination of the device characteristics plus 
imaging with a high-resolution STEM ruled out 
the existence of any Nb shorts over the top of the 
membrane. Additionally, modulation of the critical 
current versus applied magnetic flux was entirely 
consistent with the device geometry. 

The next step for this work is to fabricate these 
devices using high-mobility Si. Note that conven¬ 
tional uniformly doped Si cannot provide the nec¬ 
essary characteristic lengths as its mobility is too 
low. Therefore, we are turning to chemical-vapor 
deposited epitaxial layers of Si/Si;B grown on 
(110) oriented wafers which can have mobilities 
approaching 10* cm*/V-sec. This work is presently 
underway. 


trcpically Etched in (110) Silicon,' submitted for publication; G.E. Rittenhouse, Henry I. Smith, and J.M. Graybeal, 
'A Silicon Weak-Link Josephson Junction Using a Vertical Membrane Structure,* forthcoming. 
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Figure 2. Current-voltage characteristics of one device versus temperature. Note presence of a finite Josephson 


current. 
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4.1 Facility for the Gas Source 
Epitaxy of Compound 
Semiconductors 

Sponsors 

3M Company Faculty Development Grant 
AT&T Research Foundation 
Special Purpose Grant 
Charles S. Draper Laboratories 
Contract DL-H-418484 
Defense Advanced Research Projects Agency 
Subcontract 216-25013 and 542383 
Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
Contract DAAL03-92-C-0001 
National Science Foundation 
Grant ECS 88-46919 
Grant ECS 89-05909 

The new laboratory for the growth of both ll-VI 
and lll-V compound semiconductor heterostruc¬ 
tures, using gas source epitaxy techniques, has 
received all hardware necessary to complete the 
multichamber system. The epitaxy system consists 
of six interconnected ultrahigh vacuum chambers: 
two introduction chambers, a transfer chamber, a 
ll-VI-dedicated reactor, a Ill-V-dedicated reactor, 
and an analytical/metalization chamber. The 
analytical/metalization chamber is now attached to 
the transfer chamber; figure 1 shows a photograph 


of the new chamber. In situ analytical tools which 
are available include Auger electron spectroscopy, 
reflection high energy electron diffraction 
(RHEED), additional ports for laser illumination of 
the sample, and a three-source electron beam 
evaporator for metals. The ultrahigh vacuum inter¬ 
connection to both of the growth systems allows 
surface analysis of the many ll-VI, lll-V, and 
ll-VI/lll-V heterostructures which are envisioned, 
immediately following fabrication, without concern 
of contamination from atmospheric gases. The 
analysis chamber is expected to provide additional 
insight into the effect of both electron and photon 
beams on the surface of films grown using 
metalorganic molecular beam epitaxy. Metal 
evaporators will enable us to probe the initial 
stages of ohmic contact formation specifically with 
regard to the difficulties experienced in forming an 
ohmic contact to p-type ZnSe. The presence or 
absence of mixed interfacial layers in multilayered 
ll-VI/lll-V heterostructures may also be detected 
(in some instances) in the anal^ical chamber. 

The final addition was the connection of a gas 
source molecular beam epitaxy (GSMBE) reactor 
to a small buffer chamber attached to the transfer 
chamber. The new reactor, shown in figure 2, was 
purchased from Instruments S.A., Riber Division, 
and is now in the final stages of installation. The 
gas source reactor has solid elemental sources of 
Ga, In, Al, Si, and Be, and gaseous hydride 
sources of arsenic and phosphorus. The purpose 


' Physics Department, University of Puerto Rico, Rio Piedras, Puerto Rico. 
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Figure 1. Photograph of the analytical/metalization chamber. Connection to the ultrahigh vacuum transfer 
chamber is yet to be completed. 


of the lll-V chamber is two-fold: (1) layers con¬ 
taining P and As will be crucial for utilizing lattice- 
matched buffer layers for the ll-VI-based quantum 
well structures which are under investigation, and 
(2) Ill-V-based optical and electronic devices will 
be explored including quantum-effect electronic 
devices (in collaboration with Professors Dimitri A. 

Antoniadis and Henry I. Smith), and channel drop¬ 
ping filters and waveguide structures (in collab¬ 
oration with Professor Hermann A. Haus). In 
addition, the lit-V chamber will be employed for 
fabrication of ll-VI/lll-V multilayered heterostruc¬ 
tures which contain lll-V layers having substantial 
thickness. The II-VI reactor, which has gaseous 
sources for In, Ga, and As, will only be used for 
the growth of very thin lll-V layers on the order of 
hundreds of angstroms or less. 

Optical sources operating in the blue portion of 
the visible light spectrum have many potential 
commercial and military applications in, for 
example, optical information storage systems. 


4.2 Metalorganic Molecular 
Beam Epitaxy (MOMBE) of 
ZnSe 

Sponsors 

Charles S. Draper Laboratories 
Contract DL-H-418484 
Defense Advanced Research Projects Agency 
Subcontract 542383 and 5300716-07 
Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
Contract DAAL03-92-C-0001 
National Science Foundation 
Grant ECS 88-46919 
U.S. Navy - Office of Naval Research 
Contract N00014-88-K-0564 
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Figure 2. Photograph of the new Ill-V-dedicated GSMBE (which is attached to the small loadlock/buffer 
chamber). Final connection to the hemispherical ultrahigh vacuum transfer chamber can be seen. 


medical imaging systems, electroluminescent flat alloy buffer layers such as InGaP and GaAsP. The 
panel displays, and color laser printers. attached lll-V gas source molecular beam epitaxy 

reactor will be the source of any Ill-V-based 
Although there are a variety of material systems, lattice-matched buffer layer. Although the lattice 

such as SiC, frequency doubled Ill-V-based lasers, mismatch difficulty is circumvented here, we are 

and GaN, which can potentially offer a source of 3(50 fabricating 'epilayer-on-epilayer" structures, 

blue light, each of these systems has some limita- These epilayer-on-epilayer structures are necessary 

tions relative to integration and/or fabrication. The jq provide the minimum possible number of 

ll-VI family of semiconductors (specifically the Zn defects created at the heterovalent interface and 

chalcogenides) offers a solution to the problem of ultimately within the active region of the ll-VI 

integration, and impressive recent gains have been optical device. In addition, a key advantage 

achieved in fabrication technology. ll-VI-based offered by gas source epitaxy techniques is the 

quantum well structures exhibiting excellent costal ability to control, via mass flow controllers, the 

quality are routinely fabricated on technologically f|yx of species which exhibit very high vapor pres- 

important lll-V substrates, suggesting that future sures as a function of temperature. Sulfur and 

integration needs will be straightforvyard. phosphorus containing materials are excellent can- 

However, some substantial problems still remain to didates for fabrication utilizing gas source epitaxy 

be solved. techniques with impressive successes reported. 

In order to meet the requirements for room temper- j^e objectives of the research described here 

ature, blue optical sources, the ll-VI materials will revolve around the growth of various ll-VI-based 

require greater concentrations of sulfur to increase heterostructures using gaseous source epitaxy 

the bandgap. However, the alloying of ZnSe with including metalorganic molecular beam epitaxy 

sulfur results in a contraction of the lattice param- (MOMBE), chemical beam epitaxy (CBE), and gas 

eter and, therefore, an increase in the lattice mis- source molecular beam epitaxy (GSMBE). The 

match with binary compounds such as GaAs. This various techniques utilize combinations of gaseous 

problem is eliminated by utilizing lattice-matched 
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sources such as metalorganics and hydrides, as 
well as solid elemental sources. Thus far we have 
explored the MOMBE of ZnSe using diethylzinc 
and diethylselenium. 

A large number of films have been grown with the 
details of the growth parameters described below. 
ZnSe epilayers were normally grown on carbon- 
controlled semi-insulating (001) GaAs wafers, 
although n-type, p-type and 2°-off n-type GaAs 
were also used. Typical growth temperatures 
ranged from 250°C to with the majority of 

the film grov»rths performed near 320'’C. The Vl/ll 
flow ratio was varied from 0.2 to 5.0, with 
diethylzinc (DEZn) and diethylselenium (DESe) 
gas flows ranging from 0.5 to 2.5 seem. Reactor 
pressure during growth was typically 1x10"^ Torr. 
In nearly all cases the DEZn cracker was set at 
50°C, and the DESe was thermally decomposed at 
800°C. Cracking experiments with the quadrupole 
mass spectrometer (QMS) indicated that DESe 
decomposed at 800°C, while DEZn thermally 
dissociated at temperatures above 300°C. 

Preliminary microstructural and optical character¬ 
ization of the ZnSe suggests that the thin films are 
of good quality. A typical RHEED pattern 
observed following growth shows the presence of 
narrow streaks and readily apparent Kikuchi bands 
suggesting the existence of single crystalline films 
with smooth surface morphology. Secondary ion 
mass spectroscopy (SIMS) results indicate very 
low levels of carbon incorporation which is near 
the detectability limit of our instrument (VG-Fisons 
Instruments Model IX70S SIMSLAB). X-ray dif¬ 
fraction rocking curve analysis confirms the 
pseudomorphic nature of thin 10OOA films of ZnSe 
on GaAs. Typical full width at half maximum 
measurements for the (400) reflection is around 
220 arc seconds. Optical characterization from 
photoluminescence measurements indicates the 
presence of bandedge excitonic features which are 
shifted to higher energy reflecting the presence of 
strain, in addition to features originating from deep 
level defects. We speculate that the presence of 
these particular deep level features is due to 
stoichiometry such as Zn vacancies, and we are 
therefore investigating growth with various flux 
ratios. 

The growth rate using DEZn and DESe has con¬ 
tinued to be very low. Although the actual species 
at the surface front are not known, we expect to 
find ethyl species of Zn and Se, in addition to 
dimers of Se from the thermal decomposition of 
the metalorganic. The growth temperatures inves¬ 


tigated and the cracking data suggest that there is 
sufficient thermal energy to decompose the Zn 
molecules. (The RHEED patterns exhibit recon¬ 
structions of a Zn-stabilized surface.) In addition, 
the nearly constant growth rate observed for a 
wide variety of growth conditions (flow ratio and 
rate, substrate temperature, and with and without 
hydrogen) suggests that the growth is surface 
kinetics limited. Similar to metalorganic vapor 
phase epitaxy (MOVPE), although lacking the gas 
phase reactions, the MOMBE growth technique is 
expected to exhibit complex surface chemical 
reactions. We speculate that at the necessary low 
growth temperatures employed here, lattice sites 
for incorporation of Zn and Se are being blocked 
by alkyl radicals. The alkyl radicals are a copious 
byproduct of the surface pyrolysis of DEZn and 
thermal decomposition of the DESe. At higher 
growth temperatures (^425‘’C), where one might 
expect the alkyl radicals to desorb, competing 
desorption of Zn and Se limits the growth rate. 
Previous work^ reported using these metalorganic 
sources show growth rates approximately an order 
of magnitude greater under similar growth condi¬ 
tions. Differences which may exist include (1) 
pressure during growth, (2) use of hydrogen 
carrier gas, or (3) the actual design of the gas 
crackers. 

4.3 Photo-Assisted MOMBE of 
Wide Bandgap ll-Vl Compound 
Semiconductors 

Sponsors 

Defense Advanced Research Projects Agency 
Subcontract 530-0716-07 
Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
Contract DAAL03-92-C-0001 
National Science Foundation 
Subcontract DMR 90-07890 

Photo-assisted heteroepitaxial growth of ZnSe on 
GaAs by MOMBE is currently under intense inves¬ 
tigation. MOMBE has several potential advan¬ 
tages when compared to other epitaxial growth 
techniques such as metalorganic vapor phase 
epitaxy (MOVPE) and molecular beam epitaxy 
(MBE). The advantages include: (1) The flux 
ratio of gaseous metalorganics can be precisely 
controlled using mass flow controllers; (2) in-situ 
analysis techniques such as reflection high-energy 
electron diffraction (RHEED) and quadrupole mass 


2 H. Ando, A. Taike, R. Kimura, M. Konagai, and K. Takahashi, 'Metalorganic Molecular Beam Epitaxial Growth of 
ZnSe Using Diethylzinc and Diethylselenide,' Japanese J. Appi Phys. 25(4): L279-L281 (1986). 
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spectroscopy (QMS) are available; and (3) surface 
chemistry allows for fine tuning of the growth 
conditions. Photo-assisted MOMBE adds an extra 
degree of versatility because both selective deposi¬ 
tion of material and modification of the surface 
chemistry are possible. 

Photo-assisted growth was performed by selec¬ 
tively illuminating a portion of the wafer's surface. 
Radiation from a 5 Watt Ar ion laser (Spectra- 
Physics, Model 2025) was introduced at a 45* 
angle to the substrate normal. The laser spot on 
the wafer was typically an ellipse with major and 
minor axes of 1.4 and 1.0 cm. The effect of wave¬ 
length on growth has been studied by using laser 
radiation ranging from 780 nm to 350 nm. To 
produce the 780 nm line, a Tiisapphire crystal was 
pumped with all lines of the argon ion laser. The 
power density was kept low at approximately 
180mW/cm^ thus allowing thermal heating 
effects to be neglected. 

As mentioned in Section 4.2, under the growth 
conditions described here, the growth rate during 
MOMBE of ZnSe has consistently remained very 
low. The growth rate without laser illumination is 
typically 250 to 375A/hour. As a function of cali¬ 
brated growth temperature, the unilluminated 
growth rate shows a remarkable temperature 
insensitivity over the temperature range of 
250-425*C. This is quite unexpected since 
MOMBE has typically demonstrated a strong tem¬ 
perature dependence.3 Laser illumination, however, 
has been found to ameliorate the low growth rate. 
For the portion of the film grown with photo¬ 
assistance, a significant temperature dependence is 
observed which is more typical of MOMBE 
growth. 

Figure 3 shows the ratio of illuminated to unillumi¬ 
nated growth rates (enhancement ratio) as a func¬ 
tion of laser wavelength. At a growth temperature 
of 320°C the bandgap energy of ZnSe corresponds 
to a photon energy of ~520 nm.* The 780 nm 
radiation was found to produce no growth rate 
enhancement, while the higher energy photons 
produced significant growth rate enhancement. 
Furthermore, the growth rate enhancement was 
found to be directly proportional to the laser inten¬ 
sity (see insert of figure 3). When the TEMoo 
mode of the argon ion laser was used, the 



Figure 3. Growth rate enhancement ratio (ER) as a 
function of laser wavelength. The DESe and DEZn 
flow rates are both 0.5 seem, and the substrate temper¬ 
ature is ~320*C (arrow indicates ZnSe bandgap). The 
insert shows the growth rate enhancement ratio of the 
illuminated region (457.9 nm laser line) as a function of 
spot radius (mm). The dotted curve is a Gaussian fit 
to the data suggesting a direct correlation between 
growth rate enhancement and power density. 

enhancement ratio (solid circles) across the illumi¬ 
nated growth region was found to correlate well 
with the Gaussian intensity profile of the laser. 
The largest enhancement in the growth rate has 
been a factor of 9, although we have not yet 
attempted to optimize the various growth para¬ 
meters to maximize the growth rate. 

Thus far our data suggest that the creation of 
electron/hole pairs are required for the growth rate 
to be enhanced during photo-assisted growth. 
Photons having energies above the bandgap of 
ZnSe (at the growth temperature) resulted in 
growth rate enhancement, while photons with 
longer wavelengths do not increase the growth 
rate. Similar dependencies have been observed by 
Yoshikawa et al* for growth of ZnSe by photo- 
assisted MOVPE. For the metalorganic molecules 
utilized here, gas phase absorption^ is expected to 
occur at wavelengths much shorter (248 nm for 
DEZn) than investigated here. Although a more 


) H. Ando, A. Taike, R. Kimura, M. Konagai, and K. Takahashi, 'Metalorganic Molecular Beam Epitaxial Growth of 
ZnSe Using Diethylzinc and Diethylselenide,' Japanese J. Appl. Phys. 25(4); L279-L281 (1986). 

* A. Yoshikawa, T. Okamoto, T. Fujimoto, K. Onoue, S. Yamaga, and H. Kasai, 'Ar-lon Laser-Assisted MOVPE of 
ZnSe Using DMZn and DMSe as Reactants,' Japanese J. Appl. Phys. 29(2): L225-L228 (1990). 

* H. Hou, Z. Zhang, U. Ray, and M. Vernon, 'A Crossed-Laser Beam Study of the Photodissociation Dynamics of 
Zn(C 2 H 5)2 and (Zn(C 2 H 6 ) 2)2 at 248 nm and 193 nm,'J. Cham. Phys. 92(3); 1728-1746 (1990). 
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detailed investigation is required, we speculate that 
a potential mechanism is that a photo-generated 
carrier (we suspect holes) eliminates an ethyl- 
metal bond at the surface leading to a higher 
growth rate. Preliminary data (not shown here) 
provide evidence that there is a dependence of the 
growth rate enhancement on the gas flow ratio; 
this suggests that the growth rate enhancement is 
not due to liberation of surface alkyls which are 
blocking lattice sites by photo-generated carriers. 
The photo-assisted growth experiments currently 
underway are expected to aid our understanding of 
the dominant surface chemistry and growth mech¬ 
anisms. By adjusting conditions such that growth 


will only occur in regions irradiated by the laser, 
selective area epitaxy of ZnSe will be possible. 


4.4 Publications 

Coronado, C.A., E. Ho, L.A. Kolodziejski, and C.A. 
Huber. 'Photo-Assisted Molecular Beam 
Epitaxy of ZnSe.' Submitted to Appl. Phys. 
Lett. 

Kolodziejski, L.A. 'Chemical Beam Epitaxy of II-VI 
Semiconductors.' Invited paper presented at 
the American Vacuum Society Meeting, Seattle, 
Washington, November, 1991. 
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5.1 Electromigration and 
Microstructure 

Electric-current-induced atomic self-diffusion, or 
electromigration, causes open and short circuit 
interconnect failures in integrated circuits. Electro¬ 
migration occurs primarily by diffusion along grain 
boundaries, and interconnect failure initiates at 
grain boundaries, so that interconnect reliability is 
strongly dependent on the grain sizes and orien¬ 
tations in polycrystalline interconnects. The 
relationships between these microstructural char¬ 
acteristics and interconnect reliability have become 
more complex as interconnect dimensions (widths 
and thicknesses) have become comparable to 
grain sizes. We have studied the effects of micro¬ 
structure on the reliability of submicron and larger 
interconnects as well as vias and contacts. We 
have done this by using a variety of techniques for 
controlling grain sizes and grain size distributions 
in thin films. We are also developing models for 
microstructural evolution in order to develop accu¬ 


rate models for electromigration in submicron 
interconnects. 


5.1.1 Control of Microstructures in 
Aluminum Films 

Sponsor 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 

Project Staff 

Hai P. Longworth, Professor Carl V. Thompson 

In previous years, we reported the use of 
precipitate-induced abnormal grain growth to 
produce grains as large as 1 mm in 1 ^tm thick 
films. In the past year, we completed an extensive 
study in which we have established the require¬ 
ments for precipitate-induced abnormal grain 
growth in thin films and have demonstrated similar 


^ Presently at the University of Michigan, Ann Arbor, Michigan. 

* Ural Polytechnical Institute, Sverdlovsk, Russia. 

* Dartmouth College, Hanover, New Hampshire. 

* Hitachi Limited, Tokyo, Japan. 

■ IBM Corporation, Yorktown Heights, New York. 

* Formerly of Dartmouth College, currently with Digital Equipment Corporation, Hudson, Massachusetts. 
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phenomena in Al-Cu-Cr, Al-Ag-Cr, Al-Mn-Cr, and 
Al-W alloys. We have found that to optimize the 
precipitate-induced abnormal grain growth 
process, the alloy additions should initially be 
present as pure layers in the middle of the Ai film. 
In alloys containing Cr, the Cr serves to aid in the 
nucleation of precipitates, leading to localization of 
the precipitates in the middle of the film, and 
leading to a larger number of smaller, more closely 
spaced precipitates. Through controlled abnormal 
grain growth, the microstructure of thin films can 
be controllably varied over broad ranges. 


5.1.2 Reliability of Vias 

Sponsor 

Semiconductor Research Corporation 
Project Staff 

Harold Kahn, Professor Carl V. Thompson 

Over the past year, we have developed a via test 
structure that allows statistical characterization of 
via reliability by stressing multiple vias in parallel. 
We have used this process to characterize the reli¬ 
ability of vias with W plugs and have shown that 
under normal test conditions, joule heating of the 
plugs contributes to accelerated open failures of AI 
interconnects. We are developing simple thermal 
models to evaluate the effect of Joule heating at 
service conditions. We have also demonstrated 
that the grain size in the AI interconnect greatly 
affects via reliabilities and also the statistics of 
failure. Large grains lead to bimodal failure dis¬ 
tributions. We are developing tools for statistical 
analysis of via reliability. 


5.1.3 Modeling of Microstructures 
in Interconnects 

Sponsors 

Semiconductor Research Corporation 
National Science Foundation 


Project Staff 

Derek T. Walton, Dr. Harold J. Frost, Professor 
Carl V. Thompson 

Over the years, we have developed and tested 
models for crystal nucleation and growth to form 
continuous polycrystalline films. We have also 
modeled grain growth in continuous polycry¬ 
stalline films. In the past year we have modifi^ 
these models to investigate the evolution of the 
microstructures of polycrystalline interconnects to 
"bamboo* structures, as illustrated in figure 1. We 
find that the approach to a fully bamboo structure 
is exponential in time and that the rate of transfor¬ 
mation is inversely proportional to the square of 
the strip width. When the simulation is extended 
to model grain boundary pinning due to grooving 
at grain boundary-free surface intersections, we 
find that there exists a maximum strip width-to- 
thickness ratio (approximately 3) beyond which 
the transformation to the bamboo structure does 
not proceed to completion. In earlier work, we 
have shown that reliability of nearly-bamboo and 
fully-bamboo lines is greatly different. Our mod¬ 
eling work is now allowing us to identify the con¬ 
ditions under which it is reasonable to expect to 
obtain only fully-bamboo lines. 
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Figure 1. Evolution of the grain structure within a 
strip. The transformation to a ‘bamboo* structure is 
dominated by the growth of bamboo sections into 
regions of grain clusters. The rate of this transforma¬ 
tion depends strongly on the strip width. 
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6.1.4 Electromigration in Bicrystala 

Sponsors 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
Semiconductor Research Corporation 

Project Staff 

Hai P. Longworth, Young-Chang Joo. Brett D. 
Knowiton, Professor Carl V. Thompson 

In previous work we have shown that individual 
grain boundaries greatly affect the reliability of an 
interconnect line. In order to understand the reli¬ 
ability of interconnects, it is therefore important to 
understand the reliability of individual grain 
boundaries. To do this we have developed a new 
technique that allows us to characterize the reli¬ 
ability of populations of lines with single, identical 
grain boundaries with controlled types, locations, 
and orientations. The technique is schematically 
shown in figure 2. We start by making NaCI 
bicrystals of the desired type, on which we grow 
epitaxial bicrystal Al films which are transferred 
and sintered to oxidized silicon wafers. We then 
use conventional lithographic techniques to 
pattern multiline test structures with the same con¬ 
trolled Al grain boundary running through all the 
lines. We can then investigate, as an example, the 
effect of the boundary orientation on the reliability 
of the line, as illustrated in figure 3. 

We have shown that (1) failure times of lines con¬ 
taining single identical boundaries are lognormally 
distributed, (2) the medium time to failure 
depends more strongly on the boundary orien¬ 
tations than the types of grain boundaries, (3) the 
deviation in the time to failure has a large compo¬ 
nent not dependent on microstructure, and (4) 
both interfacial diffusion and grain boundary dif¬ 
fusion appear to contribute to failure in bicrystal 
lines and probably in bamboo and near-bamboo 
lines as well. 



Figure 2. The fabrication process of bicrystal test pat¬ 
terns. (a) NaCI bicrystal substrate made from two NaCI 
single crystals; (b) Al bicrystal thin film epitaxially 
deposited on the NaCI substrate; (c) Al film transferred 
to an oxidized Si wafer; and (d) patterning. 
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Figure 3. Failure distributions of £13[100] bicrystal 
lines as a function of 8, the angle between the current 
direction and the grain boundary plane. 
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1.1 Project Description 

Sponsor 

Joint Services Electronics Program 
Contract DAAL03-92-C-0001 

We are investigating equilibrium properties of 
quantum dots (isolated metallic systems with sizes 
smaller than the typical length of an electron 
dephasing). Until recently, we concentrated on 
such properties of an individual dot as capacitance 
and electric and magnetic susceptibility. The most 
dramatic manifestation of the quantum nature of 
the magnetic polarizability takes place for a mul> 
tiply connected geometry, e.g., for a ring, and is 
known as a persistent current. This current is a 
metallic system response to the magnetic flux 
through the hole. It can exist even when there is 
no field in the material, and it does not dissipate. 

Due to the impurity or surface scattering of the 
conductivity electrons, there is certain dispersion 
of the properties in the ensemble of macro- 
scopically equivalent samples. Therefore, we must 
consider the sample to sample (mesoscopic) fluc¬ 
tuations as well as ensemble averaged properties 
such as persistent current. 

Earlier^ we developed a theory on the persistent 
current neglecting the interaction between the 
electrons. This problem involves two energy 
scales: (1) the mean energy spacing between the 
exact one electron states in a random potential A, 
and (2) the inverse typical time of an electron dif¬ 
fusion around the sample h/t - Ec also known as 
the Thouless shift. The dimensionless ratio EJA is 


the sample conductance in natural units e^/h. We 
found that the mean root square of the mesoscopic 
fluctuations of the persistent current is of the order 
of Ece/h, while the averaged persistent current 
although nonzero, is much smaller and determined 
by A. 

The persistent current was recently observed^ and 
found to be orders of magnitude larger than the 
theoretical predictions in the absence of the inter¬ 
action between electrons. Also, superconduc¬ 
tive-like fluctuations cannot explain why the effect 
is so large. On the other hand, the magnitude of 
the observed persistent current is more or less con¬ 
sistent with the charging energy e^/R, where R is 
the radius of a ring. This was one of the moti¬ 
vations for us to study dielectric properties of 
mesoscopic systems. 

We have developed a technique to investigate the 
local charge fluctuations and local electric fields 
within a mesoscopic system. We have found that 
due to quantum interference, there are long range 
correlations in these local fluctuations (similar to 
the well-known Friedel oscillations). These corre¬ 
lations determine the additional energy caused by 
the Coulomb interaction between the fluctuations, 
and, therefore, their contribution to the thermo¬ 
dynamics. 

We have calculated the distribution of electric 
fields outside a disordered, neutral quantum dot as 
well as the quantum correction to its polarizability. 
Although both of the effects are small, they are 
very important ideologically from the point of view 
of the scaling theory of metal-insulator transition, 
which we believe are observable. 


' B.L. Altshuler, Y. Gefen, and Y. Imry, Phys. Rev. Lett. 66:88 (1991); B.L. Altshuler and B.Z. Spivak, Zh. Eksp Teo. 
92:607 (1987). 

2 L.P. Levy, G. Doland, J. Dunsmuir, and H. Bouchiat, Phys. Rev. Lett. 64: 2074 (1990); V. Chandrasekhar, R.A. 
Webb, M.Y. Brady, M.D. Ketchen, W.J. Gallagher, and A. Kleinsasser, Phys. Rev. Lett. 67; 3578 (1991). 
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These small effects are caused by the screening of 
the Coulomb interaction. Due to this screening, 
additional energy is determined only by the charge 
distribution in a layer near the sample surface of a 
width of the screening length. As a result, inter¬ 
action between the charge fluctuations is not suffi¬ 
cient to account for the persistent current. On the 
other hand, there is another magnetic flux depen¬ 
dent contribution to the Coulomb energy: the 
interaction between the charge fluctuations and 
the confining potential in a quantum dot (e.g., a 
finite work function of a metallic grain is deter¬ 
mined by a dipole layer on the surface). The 
energy of the local charge fluctuations in the field 


of this layer apparently depends on the magnetic 
field and explains the observed persistent current 
This interaction can have a number of important 
and observable consequences such as substantial 
mesoscopic fluctuations of the work function and 
its dependence on the external parameters such as 
low magnetic field. 

Publications 

Altshuler, B.L. 'Thermodynamics of Mesoscopic 
Systems.* Paper presented at the International 
Symposium on Nanostructures and Mesoscopic 
Systems, Santa Fe, New Mexico, May, 1991. 
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2.1 Project Description 

Sponsors 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
Contract DAAL03-92-C-0001 
National Science Foundation 
Grant ECS 88-13250 

Several years ago, we discovered, completely 
unexpectedly, new behavior in very small Si 
MOSFETs. Whereas conventional transistors turn 
from the off state to the on state only once as the 
gate voltage increases, these small transistors turn 
on and off periodically, in some cases as many as 
one hundred times. We have since learned how to 
control the period of these oscillations in GaAs 
field effect devices and have shown that a nm-size 
transistor turns on and off once for every electron 
added to it. Because of this, we call such a device 
a Single Electron Transistor. 

Having established that each cycle corresponds to 
the addition of one electron, the period in gate 
voltage provides a very precise measure of the 
capacitance of our devices. Because it is very 
small, '^lO'^'F, such transistors may eventually 
find application in ultra-sensitive charge detection. 
Furthermore, capacitive coupling may make pos¬ 
sible arrays of single-electron transistors with 
unique properties. Unfortunately, the single elec¬ 
tron behavior has only been observed, so far, 
below about 1 K. Whether this behavior can be 


observed at higher temperatures in the future 
depends on the details of the mechanism that 
causes it. Our objective is, therefore, to better 
understand that mechanism. 

Our initial discovery was made^ with inversion 
layers in Si which were about 30 nm wide and 
several mm long. These were dual gate MOSFETs 
with a narrow slot, <^600 nm wide, in the lower 
gate, fashioned with X-ray lithography. This lower 
gate was made of refractory metals allowing a high 
temperature anneal after all unconventional lithog¬ 
raphy was complete. In the end, these devices 
were the narrowest transistors ever made, but 
whereas most previous narrow devices had low 
mobility, these had mobilities comparable to wide 
devices. The successful fabrication of these ultra¬ 
narrow devices was the result of a close collab¬ 
oration between Professor Henry I. Smith, 
Professor Dimitri A. Antoniadis, and Professor 
Marc A. Kastner and their students. 

Soon after these devices were first fabricated, it 
was discovered that their conductance oscillated 
periodically as a function of the density of elec¬ 
trons in the conducting channel, proportional to 
the gate voltage, Vg. To prove that the period of 
the oscillations corresponded to the addition of a 
single electron, new devices had to be made. 

The first step was to produce devices close in their 
operation to a MOSFET but in GaAs where the 
influence of interface charges is much weaker.’ 
This was accomplished through a remarkably suc- 


' IBM Thomas J. Watson Research Canter, Yorktown Heights, New York. 

’ J.H.F. Scott-Thomas, S.B. Field, M.A Kastner, H.l. Smith, and D.A. Antoniadis, ‘Conductance Oscillations Peri¬ 
odic in the Density of a One-Dimensional Electron Gas,‘ PAys. Pev. Latt. 62: 583 (1989); S.B. Field, M.A. Kastner, 
U. Meirav, J.H.F. Scott-Thomas, D.A Antoniadis, H.l. Smith, and S.J. Wind, ‘Conductance Oscillations Periodic 
in the Density of One-Dimensional Electron 6ases,‘ Phys. Pav. B 42; 3523 (1990). 

3 U. Meirav, M.A. Kastner, M. Heiblum, and S.J. Wind, ‘One-Dimensional Electron Gas in GaAs: Periodic Conduc¬ 
tance Oscillations as a Function of Density,' Phys. Rtv. B (Rapid Comm.) 40: 5871 (1989). 
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Figure 1. Conductance in units of the quantum of conductance e^h versus voltage on the n-^GaAs substrate, the 
'back gate.’ On the logarithmic plot it is clear that the transistor has an on-to-off conductance ratio that exceeds 10^ 
for the lowest gate voltages. The measurements were made at 30 mK with a source-drain voltage of only 2.5 pV. 


cessful collaboration with M. Heiblum of IBM. 
The new structure was created by growing AIGaAs 
on a conductive substrate, followed by a layer of 
undoped GaAs. The electron density at the 
GaAs/AIGaAs interface, inverted from the usual 
configuration, was varied by applying a voltage to 
the substrate. Heiblum was able to grow such 
inverted heterojunctions with mobilities of 
~ 500,000cm2/V-s. 

Next the conducting channel is defined by depos¬ 
iting a metal gate on top of the GaAs, and a gap is 
patterned in the gate by electron beam lithography. 
Because a depletion region is created under the 
gate, the electrons accumulate only under the gap. 
in order to emulate the effect of two charg^ 
impurities, two constrictions are patterned in the 
gap. The electron beam lithography was done at 
IBM with the help of Shalom Wind. 


Devices of this kind worked exactly as we hoped 
they would.* Periodic oscillations of the conduc¬ 
tance were seen again, but now they were con¬ 
trolled: The period was the same for different 
devices with the same spacing between con¬ 
strictions and was larger when the spacing was 
shorter, consistent with the idea that the same 
number of electrons is added for one period in all 
devices. Calculations of the capacitance of such 
devices are only consistent with each period corre¬ 
sponding to one electron added per oscillation. 

This is an amazing result: Figure 1 shows that the 
conductance consists of periodic, narrow, well- 
separated resonances. The conductance rises and 
falls by a factor '^1000 with a variation of gate 
voltage that corresponds to only >^1 /3 of an elec¬ 
tron. 


* U. Meirsv, M.A. Kastner, and S.J. Wind, 'Single Electron Charging and Periodic Conductance Resonances in 
GaAs Nanostructures,* Phys. P»v. L»tt. 65: 771 (1990). 
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We now understand quite well why the Single 
Electron Transistor works the way it does. The 
conductance results from resonant tunneling, and 
Professor Patrick Lee and his collaborators have 
shown that the states arise from single particle 
quantum states (the Fabry-Perot states induced by 
the two barriers, the same as expected for any 
'quantum dot'), but that the energies of the levels 
are dominated by the Coulomb interaction 
between the electrons in the isolated segment. 

Experiments in the quantum Hall regime* at first 
seemed to indicate that the nature of the quantum 
states was very simple. By measuring the position 
in gate voltage of a series of conductance reso¬ 
nances, we were able to measure the energies for 
adding a series of electrons as a function of mag¬ 
netic field. Figure 2 shows (top panel) how the 
peak positions oscillated with magnetic field. 
Subtracting a constant Vg difference from each pair 
of cunres we emulated the effect of subtracting a 
constant Coulomb energy. We thereby obtained 
the level spectrum shown in the lower panel. 

The level spectrum looks remarkably similar to that 
for non-interacting electrons in a magnetic field 
and a confining electrostatic potential. We 
thought, therefore, that the simple Coulomb 
blockade model, with a constant Coulomb inter¬ 
action between any two electrons in the isolated 
region, was adequate. However, more careful 
analysis of data like that in figure 2 shows that this 
simple picture is inadequate. In particular, the dif¬ 
ference in energies for electron spin parallel and 
antiparallel to the magnetic field are much too 
large in the experiment. 

Recently, a density functional model for the iso¬ 
lated region has been developed in collaboration 
with Professor Lee's group. The good agreement 
between the model and the experimental results 
shows that electrons in different spatial and even 
spin states have very different Coulomb inter¬ 
actions. This is a major step in deepening our 
understanding of the physics of single electron 
transistors. 


Publications 

McEuen, P.L., E.B. Foxman, U. Meirav, M.A. 
Kastner, Y. Meir, N.S. Wingreen, and S.J. Wind. 
'Transport Spectroscopy of a Coulomb Island 
in the Quantum Hall Regime.' Phys. Rev. Lett. 
66:1926 (1991). 
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Figure 2. (a) The gate voltages corresponding to a 
series of consecutive conductance peaks like those in 
figure 1, as a function of magnetic field, (b) Result of 
subtracting a constant voltage from each successive 
peak in (a). This gives the energy for adding the next 
electron as a function of field with the average 
Coulomb interaction subtracted. The arrows follow a 
particular level as it moves through successive peaks. 


• P.L. McEuen, E.B. Foxman. U. Meirav, M.A. Kastner, Y. Meir, N.S. Wingreen and S.J. Wind, 'Transport Spectro¬ 
scopy of a Coulomb Island in the Quantum Hall Regime,' Phys. Rev. Lett. 66:1926 (1991). 
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3.1 Project Description 

Sponsor 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
Contract DAAL03-92-C-0001 

We continue the study of electron transport 
through a quantum dot. It is now well established 
both experimentally' and theoretically^ that 
Coulomb blockade plays a key role in this problem 
and is responsible for the conductance peaks 
which have been observed to be periodic as a 
function of gate voltage. Each of these conduc* 
tance peaks occurs at a gate voltage where the 
energy of N and N+1 electron state are degenerate 
and each peak corresponds to the addition of a 
single electron to the quantum dot. In the past 
year, the general formalism developed by Meir et 
al.> has been extended and applied to the case 
where a strong magnetic field is applied to the 
system. The field is strong enough so that only 
the lowest two Landau levels are occupied, i.e., 
the system would be in the integer quantum Hall 
regime in the bulk. The advantage of this system 
is that the single particle states are well known to 
be edge states, and the energy of these states are 
well understood as a function of magnetic field. 
Thus, the successful application of the theory to 
explain the experimental data^ provides strong 
support for the general theory^ of Coulomb 
blockade in resonant tunneling. 

Experimentally it was found that in the presence of 
8 strong magnetic field, the location in gate 


voltage of each conductance peak exhibits small 
periodic fluctuation. At the same time, the magni¬ 
tude of each conductance peak goes through deep 
minima. This can be explained based on a model 
where the low lying states are the edge states of 
the first and second Landau levels. For a given 
conductance peak, the number of electrons on the 
quantum fluctuates between N and N+1. 
However, as the magnetic field is increased, the N 
electrons are distributed between the two Landau 
levels in a different way, and electrons are trans¬ 
ferred from the second Landau level to the first. 
The transfer occurs when the edge state of the first 
Landau level crosses the edge state of the second 
Landau level. Thus, as a function as of a magnetic 
field, we are watching the transfer of electrons 
from the second to the first Landau level, one elec¬ 
tron at a time. Furthermore, since the dependence 
of the energy of the edge state on magnetic field is 
understood, the periodic fluctuation of the reso¬ 
nant energy with magnetic field can be under¬ 
stood. The dips in the magnitude of the 
conductance peak can also be explained as being 
due to the fact that the edge state of the second 
Landau level is more tightly bound, and its cou¬ 
pling to the leads is much smaller so that a dip 
occurs every time the second Landau level moves 
through the chemical potential of the leads. Based 
on this picture, quantitative calculation of the tem¬ 
perature dependence of the peak position was 
made and successfully compared with the exper¬ 
iment. Indeed, we now have sufficient confidence 
in the model that the problem can be turned 
around, and we can use the experiment as a spec¬ 
troscopy of the single particle level, i.e., we can 


' U. Meirsv, M.A. Kastner, and S.J. Wind, P/iys. Rev. Lett 65: 771 (1990). 

2 Y. Meir, N.S. Wingraen, and P.A. Lee, 'Transport through a Strongly Interacting Electron System: Theory of Peri¬ 
odic Conductance Oscillations,' Rhys. Rev. Lett. 66: 3048 (1991). 

* P.L. McEuen, E.B. Foxman, U. Meirav, M.A. Kastner, Y. Meir, N.S. Wingreen, and S.J. Wind, “Transpon Spectro¬ 
scopy of a Coulomb Island in the Quantum Hall Regime,* P/iys. Rev. Lett. 66:1926 (1991). 
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extract the dependence of the edge state energy 

level as a function of magnetic field. 

Publications 

Feng, S., and P.A. Lee. ‘Mesoscopic Conductors 
and Correlations in Laser Speckle Patterns.* 
Sc/. 251:633 (1991). 

Kinaret, J., and P.A. Lee. ‘Conductance of a Dis¬ 
ordered Narrow Wire in a Strong Magnetic 
Field.‘ Phys. Rev. B. Forthcoming. 

Lee, P.A. ‘Small Structures: A Laboratory for 
Studying Strongly Interacting Systems,* Pto- 


ceedings of the Taniguchi Symposium on 
Mesoscopic Physics 1991. Shima, Japan. 
Forthcoming. 

McEuen, P.L., E.B. Foxman, U. Meirav, M.A. 
Kastner, Y. Meir, N.S. Wingreen, and S.J. Wind. 
‘Transport Spectroscopy of a Coulomb Island 
in the Quantum Hall Regime.‘ Phys. Rev. Lett. 
66:1926 (1991). 

Meir, Y., N. Wingreen, and P.A. Lee. ‘Transport 
Through a Strongly Interacting Electron 
System: Theory of Periodic Conductance Oscil¬ 
lations.* Phys. Rev. Lett. 66: 3048 (1991). 
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4.1 Submicron Structures 
Laboratory 

The Submicron Structures Laboratory at MIT 
develops techniques for fabricating surface struc¬ 
tures with feature sizes in the range from nano¬ 
meters to micrometers and uses these structures in 
a variety of research projects. These projects, 
which are described briefly below, fall into four 
major categories: (1) development of submicron 
and nanometer fabrication technology: (2) nano¬ 
meter and quantum-effect electronics; (3) crystal¬ 
line films on non-lattice-matching substrates; and 
(4) periodic structures for x-ray optics, spectro¬ 
scopy and atomic interferometry. 


4.2 Microfabrication at 
Linewidths of 100 nm and 
Below 

Sponsors 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
Contract DAAL03-92-C-0001 
National Science Foundation 
Grant ECS 90-16437 

Project Staff 

Martin Burkhardt, James M. Carter, William Chu, 
Kathleen R. Early, Reza A. Ghanbari, Scott D. 
Hector, Yao-Ching Ku, Alberto M. Moel, George E. 
Rittenhouse, Dr. Mark L. Schattenburg, Professor 
Henry I. Smith, Vincent V. Wong, Anthony Yen 

A variety of techniques for fabricating structures 
with characteristic dimensions of 100 nm (0.1 pm) 


^ University of Cairo, Giza, Egypt, and IBM Corporation, Thomas J. Watson Research Center, Yorktown Heights, 
New York. 


* Outokumpu Electronics, Espoo, Finland. 

3 Princeton University, Princeton, New Jersey. 
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and below are investigated. These include: x-ray 
nanolithography, holographic lithography, achro¬ 
matic holographic lithography, electron-beam 
lithography, focused-ion-beam lithography, 

reactive-ion etching, electroplating, and liftoff. 
Development of such techniques is essential if we 
are to explore the rich field of research applications 
in the deep-submicron and nanometer domains. 
X-ray nanolithography is of special interest 
because it can provide high throughput and broad 
process latitude at linewidths of 100 nm and 
below. As discussed in section 4.3, we have 
established a reliable mask technology, based on 1 
/im-thick SiNx membranes, as shown schematically 
in figure 1. Mask blanks are patterned by several 
techniques including electron-beam lithography, 
focused-ion-beam lithography, holographic lithog¬ 
raphy, and transfer from other x-ray masks. 

For high-quality x-ray lithography, the mask-to- 
sample gap, G, must be precisely controlled. The 
mesa-rim architecture of figure 1 yields optically 
flat membranes enabling us to control gap down 
to about 2 microns by means of spacer studs. For 
100 nm features a gap of 10 microns is suitable. 
However, for 50 nm features a gap of 2.5 microns 
is required. Figure 2 shows the replication of a 
planar-resonant-tunneiing field-effect transistor 
pattern having 50 nm minimum features. The gap 
in this case was 2.7 microns. This gap is signif¬ 
icantly larger than predicted by simple Fresnel dif¬ 
fraction formulae which neglect diffraction inside 
the mask ab^r'rber. Taking this into account, using 
vector elec*'7magnetic calculations and 


rasK for Soft X-Rav LUhoqraDhv 
(A.S t.3 nm) 



Figure 1. Structure of the new standard, 31 mm 
diameter, mesa-rim, x-ray mask suitable for lithography 
from 1000 to 20 nm feature sires. Membranes are flat 
to < 0.25 /im. 


Proximity X-rav Nanollthoaraphy 



Electrode 

Pattern 

PRESTFET* 


• 50 nm lines and spaces 

• k s 1.32 nm G = 2.72 pm 

G-aVIF/X 
as 1.44 

• Implication: if W s 0.1 pm, X. s 1 nm 

Gs 10 to 15 pm! 

* PRESTFET -►Planar Resonant-Tunneling 
Field-Effect Transistor 


Figure 2. Replication, using the Cut x-ray (1.3 nm) 
and a mask-substrate gap of 2.7 microns, of 50 nm line 
and space pattern. 


accounting for the finite source size, we have 
developed a simulation that predicts image con¬ 
trast as a function of feature size and gap. In order 
to increase the allowable gap for features of 70 nm 
and below, we are investigating the feasibility of 
near-field in-line x-ray holography. 

Using the figure 1 mask architecture, patterns are 
readily transferred onto standard Si wafers. 
However, in research on GaAs or other lll-V 
devices, substrates are often much smaller than the 
31 mm-diameter mask membrane. We have devel¬ 
oped an apparatus that allows us to bring the 
membrane into contact with such substrates 
during exposure, without risk of mask damage. 
Resolution down to 10-20 nm is feasible under 
such exposure conditions. 

Alignment of x-ray masks to patterns pre-existing 
on a substrate is currently done with an optical 
imaging system equipped with a video camera. 
We are investigating electronically enhanced moir4 
imaging methods in order to achieve sub-10 nm 
alignment precision. In this approach the gap is 
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set at <^4 microns and the mask translated 
piezoelectrically. 

The achromatic holographic lithography system in 
conjunction with an anti-reflection coating (ARC) 
specially formulated for the 193 nm wavelength is 
now producing high quality images, as shown in 
figure 3. The ARC consists of PMMA and a bis- 
azide, 4, 4'-diazidodiphenyl sulfone. 


Achromatic Holographic 
Lithography 



100nm-Period Grating 
in PMMA on top of ARC 


Figure 3. A 100 nm-period grating exposed in PMMA 
on top of an anti-reflection coating (ARC) by achro¬ 
matic holographic lithography using the ArF laser at 
193 nm. 

4.3 Improved Mask Technology 
for X-Ray Lithography 

Sponsors 

National Science Foundation 
Grant ECS 90-16437 
Semiconductor Research Corporation 
Contract 90-SP-080 
U.S. Navy - Naval Research Laboratory 
Contract N00014-90-K-2018 


Project Staff 

Juan Ferrera, Yao-Ching Ku, Mark K. Mondol, 
Lee-Peng Ng, Or. Mark L. Schattenburg, Professor 
Henry I. Smith, Flora S. Tsai, Sabah Yunus 

The x-ray mask architecture shown in figure 1 is 
compatible with both our research in the sub-100 
nm domain and the commercial fabrication of 
submicron integrated electronics. For the latter 
application especially, distortion in the mask must 
be eliminated. Distortion can arise from stress in 
the absorber, which is usually gold (Au) or 
tungsten (W). We have developed means of sput¬ 
tering W so that the stress is below 5 E 7 dynes/ 
cm^ and the distortion induced by stress is below 
the level of detectability in a Linnik interferometer. 
The corresponding in-plane pattern distortion is 
calculated to be less than 1 nm. We are currently 
developing an interferometer to directly measure 
in-plane distortion. 

Membranes of SiNx are preferred at the present 
time because they can be made by conventional 
LPCVD methods and are extremely rugged and 
resistant to breakage. (SiNx membranes only 1 
micron thick can even be used as vacuum 
windows, 20 mm in diameter.) Membranes are 
characterized by their optical transmission, reso¬ 
nant frequency and bulge under pressure. From 
these measurements we extract index of refraction, 
optical attenuation, thickness, stress, and Young’s 
modulus. We have also made radiation hardness 
measurements on membranes by subjecting them 
to the equivalent of one million exposures on the 
University of Wisconsin synchrotron. Effects of x 
radiation on membrane stress are extremely small, 
and we believe this can be eliminated by reducing 
the oxygen content of the SiNx. To this end a 
special purpose LPCVD reactor for low-stress SiNx 
will be set up within the Integrated Circuit Labora¬ 
tory at MIT. 

We are collaborating with the Naval Research Lab¬ 
oratory (NRL) in using their JEOL JBX5D11 
e-beam lithography system to pattern x-ray masks. 
The patterns are generated using the CAD tools at 
MIT, converted to an acceptable binary interme¬ 
diate format, transferred to NRL via ARPA Internet, 
converted there to the JEOL format, and written 
onto mask blanks that we send in advance to NRL 
(by express mail). Once the masks are written, 
they are shipped back to MIT for further pro¬ 
cessing. This collaboration has been fruitful for a 
number of projects including short-channel 
MOSFETs (section 4.4), multiple-parallel one¬ 
dimensional conductors (section 4.6), and the 
electron waveguide device (section 4.7). 
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4.4 Study of Electron Transport 
in Si MOSFETs with 
Deep-Submicron Channel 
Lengths 

Sponsors 
IBM Corporation 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
Contract DAAL03-92-C-0001 
Semiconductor Research Corporation 
Contract 90-SP-080 

Project Staff 

Professor Dimitri A. Antoniadis, Professor James 
Chung, Dr. Hao Fang, Hang Hu, Professor Henry I. 
Smith, Lisa T-F. Su 

We have continued using x-ray lithography to fab¬ 
ricate both N and P channel MOSFET devices with 
effective channel length down to 80 nm. As 
channel lengths decrease below about 150 nm, 
electron velocity overshoot, and thus an increase 
of device transconductance has been observed 
both at room temperature and liquid nitrogen tem¬ 
perature. It appears that the necessary conditions 
for this phenomenon are high carrier mobility close 
to the Si/SiOj interface, and low parasitic resist¬ 
ance. The former is achieved in both our 
NMOSFET and PMOSFET devices by utilizing e 
super steep retrograde doping of the channel. In 
the NMOSFET case, indium, by virtue of its 
heavier mass and desirable doping profile in 
silicon, has been used to give retrograde threshold 
control doping with low surface concentration and 
yet high peak concentration. In the PMOSFET 
case, arsenic has been used to achieve the same 
goals. 

Record saturation transconductance of over 
IS/mm in deep submicron NMOSFET with 
channel length down to 100 nm was obtained 
with the indium doping technique. This under¬ 
scores the achievement of increased surface 
mobility with the steeper retrograde channel 
doping. A transconductance of 0.23 S/mm was 
achieved in deep-submicron PMOSFET. 

As shown in figure 4 we have developed a tech¬ 
nology for self-aligned silicide N/PMOSFET device 
fabrication for the purpose of significantly 
reducing parasitic resistance. We have used cobalt 
deposition on the exposed silicon of source/drain 
and the exposed polysilicon of the gate electrode 
with a subsequent two-step rapid thermal 
annealing to form cobalt silicide (CoSi 2 ). Thin 
oxide spacers around the gate electrode have been 
used to prevent shorts between gate and 
source/drain. These process improvements were 



Figure 4. Configuration of self-aligned, 100 nm- 
channel-length, P-MOSFET. 


tested first with conventional lithography where 
the submicron gate was achieved by resist erosion 
in an Oj plasma after resist exposure and develop¬ 
ment. The first ultrashort channel self-aligned 
silicide MOSFET devices with channel length 
down to 0.12 /im for PMOS and 0.15 ^m for 
NMOS and parasitic resistance down to 400 ohms 
per micron width have been successfully fabri¬ 
cated. 

More recently, SiNx-membrane x-ray mask tech¬ 
nology has allowed us to use x-ray lithography for 
the definition of the gates. X-ray masks were pat¬ 
terned by electron-beam lithography at the Naval 
Research Laboratory, as shown in figure 5. Fabri¬ 
cation of both NMOSFET and PMOSFET devices 
with the new lithography technology and new 
cobalt silicide shallow junction source/drain tech¬ 
nology should give us 100 nm channel length 
CMOS circuits compatible with commercial mass 
production. 

4.5 Studies of Coulomb 
Charging in Uitrasmall 
Semiconductor Devices 

Sponsor 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
Contract DAAL03-92-C-0001 
U.S. Air Force - Office of Scientific Research 
Grant F49620-92-J-0064 

Project Staff 

Professor Dimitri A. Antoniadis, Professor Marc A. 
Kastner, Arvind Kumar, Professor Terry P. Orlando, 
Professor Henry I. Smith 
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X-ray Mask for O-I^m MOSFET 
[Written by E-Beam Lithography at NRL] 

Gate Pattern 


Gate Lengths0.2 pm 
Gate Width s 25 pm 

'igure 5. Scanning electron micrograph of x-ray mask 
for the too nm channel-length MOSFET. The e-beam 
patterning was done at the Naval Research Laboratory, 
the plating and x-ray lithography was done at MIT. 

in previous studies on very narrow silicon 
MOSFETs fabricated using x-ray lithography a 
fundamentally different mechanism to modulate 
the low-temperature electrical conduction was 
found. When the devices were cooled to the milii- 
Kelvin range, the electrical conductance exhibited 
a series of sharp oscillations which were almost 
perfectly periodic in the voltage of the gate used to 
control the electron density. This periodic modu¬ 
lation of the conductance was thought to be 
caused by presence of two dominant impurities in 
the electron channel which resulted in an isolated 
segment of the electron gas. In order to test the 
two-impurity hypothesis, the structure illustrated in 
figure 6(a) was fabricated using electron beam 
lithography and liftoff. The structure, fabricated on 
an inverted GaAs/AIGaAs heterostructure, consists 
of two Schottky gates used to squeeze the elec¬ 
tron gas laterally, while the n-*- substrate serves as 
a bottom gate to independently modulate the 
channel electron density. The two sets of con¬ 
strictions were intended to emulate the potential 
due to the two impurities in the MOSFET device. 


CoulomD Itlsns structurt utilizing an 
‘inverud* conriguratton of MSE grown layara 


Elactrodo 6ata with 



I haavlly dopad n-typa 



Gate volioge (mV from fhreshald) 


Figure 6. (a) Device structure used to study Coulomb 
blockade effects. A negative voltage on the two 
Schottky gates on the top depletes the electron gas at 
the lower GaAs/AIGaAs interface below, forming a 
narrow electron channel interrupted by two potential 
barriers at the constrictions. The n-i- GaAs substrate 
serves as a bottom gate to modulate the channel elec¬ 
tron density, (b) Periodic oscillations in the conduc- 
tence as a function of bottom gate voltage for the 
device shown in (a). 


while the mobility was sufficiently high to reduce 
the probability of unwanted accidental impurities. 

When cooled to milli-Kelvin temperatures, conduc¬ 
tance oscillations periodic in the bottom gate 
voltage were again observed, as shown in figure 
6(b). Moreover, the voltage separation between 
successive conductance peaks was found to 
depend on the separation between the two poten¬ 
tial barriers. To better understand these structures, 
we performed semi-classical calculations of poten¬ 
tial profiles and electron densities for the structure 
in figure 6(a). Equipotentiai contours near the 
AIGaAs/GaAs interface are shown in figure 7 for a 
typicai bias point with the shaded regions indi¬ 
cating areas occupied by inversion electrons. 
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Figure 7. Equipotential contours for the device in 
figure 6(a) near the lower GaAs/AIGaAs interface at a 
typical bias point. The shaded regions indicate areas 
occupied by electrons. The two sets of constrictions 
result in an isolated electron segment, or 'island,' 
through which transport can be thermally activated. 

We calculated the integrated electron charge in the 
isolated electron segment, or 'island,' between the 
sets of constrictions as a function of bottom gate 
voltage and found that it is approximately linear in 
bottom gate voltage, and that the voltage required 
to add one electron to the island is very nearly the 
spacing between successive conductance peaks. 
Thus, each conductance peak corresponds to the 
addition of a single electron to the island. This 
modulation of the conductance by changing the 
charge by a single electron can be understood in 
terms of a theory known as Coulomb blockage. 
According to the Coulomb blockade theory, a con¬ 
ductance valley corresponds to the electrostatic 
energy of the island being minimized by an integer 
number of electrons, so that addition of another 
electron costs some finite charging energy and 
transport is thermally activated. A conductance 
peak, on the other hand, corresponds to a state of 


the island where the electrostatic energy is mini¬ 
mized by a half-integer number of electrons, so 
that there is no energy barrier to adding another 
electron. These two states repeat cyclically in 
bottom gate voltage, resulting in the observed 
periodic modulation of the conductance. 

We are presently enhancing our understanding of 
Coulomb blockade effects both experimentally and 
theoretically. We first plan to reproduce the effect 
in a conventional modulation-doped GaAs/AIGaAs 
heterostructure, which is much more widely avail¬ 
able than the inverted heterostructure. Secondly, 
we are working to understand all the capacitance 
components to the island in addition to the bottom 
gate capacitance described above. This is impor¬ 
tant because the observability temperature for 
Coulomb effects depends inversely on the total 
capacitance of the island, so that reducing the 
island capacitance is key to observing these effects 
at temperatures above the milli-Kelvin regime. 
Finally, we plan to study experimentally transport 
through more than one dot as well as the coupling 
of dots in parallel. 

Because the modulation of the conductance is so 
strong. Coulomb blockade is among the most 
promising of novel effects in ultrasmall structures 
for long-term applications. 


4.6 Study of 

Quasi-One-Dimensional Wires 
and Superlattice Formation in 
GaAs/AIGaAs Modulation 
Doped Field-Effect Transistors 

Sponsor 

U.S. Air Force - Office of Scientific Research 
Grant F49620-92-J-0064 

Project Staff 

Professor Dimitri A. Antoniadis, Martin Burkhardt, 
William Chu, Professor JesOs A. del Alamo, Reza 
A. Ghanbari, Professor Terry P. Orlando, Professor 
Henry I. Smith 

We have continued to use the modulation-doped 
field-effect transistor (MODFET) to explore elec¬ 
tron back diffraction effects and reduced dimen¬ 
sionality effects in the GaAs/AIGaAs system. 
Because such devices have critical dimensions on 
the order of 100 nm, the use of x-ray lithography 
for fine feature definition has been central to our 
strategy. Recent results in our laboratory indicate 
that there may be a significant advantage in using 
x-ray lithography over direct-write e-beam lithog¬ 
raphy to fabricate quantum-effect electronic 
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Figure 11. Schematic of an electron waveguide 
device defined by the split-gate configuration. The 
shaded region represents the electron concentration at 
the AIGaAs/GaAs interface. 



Figure 12. Top view of a parent x-ray mask showing 
the pattern for a 0.1 pm-long waveguide device. 



Figure 13. Conductance of a L > 0.75 pm electron 
waveguide device as a function of split-gate voltage. 


4.8 Arrays of 

Field-Effect-Induced Quantum 
Dots 

Sponsors 

Joint Services Electronics Program 
Contract OAAL03-89-C-001 
Contract DAAL03-92-C-001 
U.S. Air Force - Office of Scientific Research 
Grant F49620-92-J-0064 

Project Staff 

Professor Dimitri A. Antoniadis, Martin Burkhardt, 
Reza A. Ghanbari, Professor Terry P. Orlando, Pro¬ 
fessor Henry I. Smith, Professor M. Shayegan,* 
Professor Daniel Tsui,* Yang Zhao* 

A metal grid on a modulation-doped AIGaAs/GaAs 
substrate (depicted in figure 14a) produces a two- 
dimensional periodic potential modulation at the 
AIGaAs/GaAs interface via the Schottky effect. If 
a gate electrode is attached to the grid, the poten¬ 
tial can be further modified with an external 
voltage source. By changing the gate voltage from 
positive to negative values, the potential seen by 
the electrons located at the AIGaAs/GaAs interface 
can be varied from uniform (in which case the 
electrons behave as a 2-D electron gas), to weakly 
coupled zero-D quantum wells (figure 14b), to 
isolated zero-D quantum dots (figure 14c). We 
have made such structures with spatial periods of 
200 nm in both orthogonal directions and covering 


* Princeton University, Princeton, New Jersey. 
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Figure 14. (a) Metal grid gate on a modulation- 

doped AIGaAs/GaAs substrate; (b) Depiction of poten¬ 
tial seen by electrons at the AIGaAs/GaAs interface for 
weakly coupled quantum dots; (c) Potential for the 
case of isolated quantum dots. 


areas of several square millimeters. The isolated 
quantum dots and the attendant zero-dimensional 
electronic sub-bands were examined in collab¬ 
oration with D. Tsui at Princeton University using 
far-infrared (FIR) cyclotron resonance. Transitions 
between the discrete energy levels in the quantum 
dots were observed as a function of magnetic filed. 
Results were in agreement with a theoretical 
model. 

Currently, we are continuing our study using 
extremely high quality samples prepared by M. 
Shayegan's group at Princeton. With mobilities 
typically greater than 10* cm^/Vsec, the resolution 
of the experiments should improve dramatically. 


4.9 Planar-Resonant-Tunneling 
Field-Effect Transistors 
(PRESTFET) 

Sponsor 

U.S. Air Force - Office of Scientific Research 
Grant F49620-92-J-0064 

Project Staff 

Professor Dimitri A. Antoniadis, Mike T. Chou, 
William Chu, Professor Henry I. Smith 

Previously, we reported on the performance of a 
planar-resonant-tunneling field-effect transistor 
(PRESTFET) depicted in figure 15, in which the 
gate electrodes were 60 nm long and separated by 
60 nm. Clear evidence of resonant tunneling 
through the bound states in the well between 
electrodes was observed, as shown in figure 
15(b). 

In order to reduce the electrode separation while 
retaining a large process latitude, we have chosen 
to pursue a new technology for making the 
PRESTFET. In collaboration with S. Rishton of 
IBM, a high-performance e-beam nanolithography 
system was used to write PRESTFET patterns on 
SiN| x-ray mask membranes, 1 nm thick. Reduced 
backscattering from the thin membrane has per¬ 
mitted us to achieve 50 nm lines and spaces. We 
hope to get even finer features in the future. The 
written masks are processed at MIT, replicated 
using the Cul x-ray (1.3 nm), and Schottky 
electrodes formed by liftoff. Masks are aligned to 
GaAs substrate using a specially built alignment 
and exposure station. 


4.10 Fabrication of 
Distributed-Feedback Lasers 
and Channel-Dropping Filters 

Sponsor 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
Contract DAAL03-92-C-0001 

Project Staff 

James M. Carter, Woo-Young Choi, Jay N. 
Damask, Professor Clifton G. Fonstad, Jr., Pro¬ 
fessor Herman Haus, Professor Leslie A. 
Kolodziejski, Professor Henry I. Smith, Vincent V. 
Wong 

We are developing techniques to fabricate large 
area, spatially-coherent gratings for both quarter- 
wave-shifted distributed-feedback (DFB) lasers 
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Figure 15. (a) Layout of a 4-terminal double-barrier 
pianar-resonant-tunneling field-effect transistor 
(PRESTFET). (b) Plot of source-drain current versus 
gate voltage for a PRESTFET with 60 nm well width. 

and channel-dropping filters (CDFs) on InP-based 
materials. The quality of these gratings is the most 
important factor determining the performance of 
these devices. Typical grating periods will be 
~200 nm, and grating lengths will range from 
^500 {im to several millimeters. Much of the 
nanofabrication technology developed in our lab 
over the years is directly applicable to the fabri¬ 
cation of these devices. SiNx x-ray masks will be 
patterned using both holographic lithography and 
e-beam lithography, as described below. X-ray 


lithography will then be used to define the device 
patterns on the InP-based material. 

Through the use of holographic lithography, large- 
area, spatially-coherent grating patterns can be 
defined. Fine control of the grating period is also 
possible. However, the incorporation of the 
quarter-wave shift into the grating structure is not 
easily accomplished by holographic means. We 
are pursuing a technique that takes advantage of a 
unique aspect of x-ray lithography: the phase of a 
segment of a grating can be shifted by pi-radians 
by transferring it onto the back surface of the mask 
membrane. Grating structures containing single 
and multiple quarter-wave shifts, which are neces¬ 
sary for both DFB lasers and COFs, will be fabri¬ 
cated using this technique. 

The incorporation of quarter-wave shifts into a 
grating structure is straightforward with e-beam 
lithography. However, precise control of the 
grating period, and large-area, spatially-coherent 
grating patterns are not easily obtained with 
e-beam lithography. For precise control of the 
grating period, a continuously-variable field size is 
necessary. To achieve spatial coherence in the 
written gratings one must eliminate the stitching 
errors, which occur when writing adjacent fields. 
We are pursuing a technique based on a global- 
fiducial grid scheme which potentially can solve 
both of these problems. Once developed, this 
technique will be incorporated into the e-beam 
lithography process to define the appropriate DFB 
laser and CDF grating patterns onto a SiN. x-ray 
mask. 

4.11 Novel Superconducting 
Tunneling Structures 

Sponsor 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
Contract DAAL03-92-C-0001 

Project Staff 

Professor John M. Graybeal, Dr. Bernard S. 
Meyerson,* George E. Ritten house. Professor 
Henry I. Smith 

In this program we seek to examine the behavior 
of short-channel hybrid Josephson devices, whose 
geometry is designed to display quantum elec¬ 
tronic interference behavior. These devices repre¬ 
sent the first attempt to observe Josephson 


B IBM Corporation. 
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coupling via resonant tunneling. The characteristic 
energy scale for such a device would be set by 
quantum confinement, and consequently may 
provide a new avenue for obtaining finite gains in 
Josephson three-terminal devices. 

These novel structures consist of an ultrathin ver¬ 
tical Si membrane (the quantum well) sandwiched 
between two superconducting Nb counterelec¬ 
trodes. These structures are fabricated using a 
combination of x-ray and optical lithography, 
anisotropic wet chemical etching, planarization 
and reactive dry etching techniques. The basic 
geometry of these new structures has now been 
fully validated using low-mobility uniformly doped 
Si wafers, where we find Josephson behavior con¬ 
sistent with the best published work on 'conven¬ 
tional' Nb/Si/Nb junctions. However, no resonant 
tunneling effects were expected or observed in 
these devices due to their low mobility. Present 
work now centers on the fabrication of similar 
structures in high-mobility modulation-doped epi¬ 
taxial layers grown on Si, where we expect reso¬ 
nant tunneling effects to appear. 

4.12 Submicrometer-Period 
Transmission Gratings for 
X-Ray and Atom-Beam 
Spectroscopy and 
Interferometry 

Sponsor 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
Contract DAAL03-92-C-0001 

Project Staff 

James M. Carter, Daniel B. Olster, Dr. Mark L. 
Schattenburg, Professor Henry I. Smith 

Transmission gratings with periods of 100-1000 
nm are finding increasing utility in applications 
such as x-ray, vacuum-ultraviolet, and atom- beam 
spectroscopy and interferometry. Over 20 labora¬ 
tories around the world depend on MIT-supplied 
gratings in their work. For x-ray and UV spectro¬ 
scopy, gratings are made of gold and have periods 
of 100-1000 nm, and thicknesses ranging from 
100-1000 nm. They are most commonly used for 
spectroscopy of the x-ray emission from high-tem¬ 
perature plasmas. Transmission gratings are sup¬ 
ported on thin (1 jum) polyimide membranes or 
made self supporting ('free standing') by the 
addition of crossing struts (mesh). (For short 
x-ray wavelengths membrane support is desired, 
while for the tong wavelengths a mesh support is 
preferred in order to increase efficiency.) Fabri¬ 


cation is generally performed by holographic 
lithography, x-ray lithography and electroplating. 
Progress in this area tends to focus on improving 
the yield and flexibility of the fabrication proce¬ 
dures. 

Another application is the diffraction of neutral 
sodium beams (de Broglie-wavelength ~17 pm) 
by mesh-supported gratings. Professor Pritchard's 
group at MIT has clearly demonstrated atomic dif¬ 
fraction and interference. Because good spatial 
coherence (low distortion) of the grating is critical 
to ensure measurable interference of the beams, 
efforts are concentrated on the use of holographic 
lithography and the reactive-ion etching of free¬ 
standing gratings in low stress and high stiffness 
materials such as silicon nitride. Results are 
shown in figure 16. 

4.13 High-Dispersion, High 
Efficiency Transmission 
Gratings for Astrophysicaf 
X-Ray Spectroscopy 

Sponsor 

National Aeronautics and Space Administration 
Contract NAS8-36748 

Project Staff 

Professor Claude R. Canizares, Robert C. Fleming, 
Dr. Mark L. Schattenburg, Professor Henry I. Smith 

This work involves a collaboration between the 
Center for Space Research and the Submicron 
Structures Laboratory (SSL), providing trans¬ 
mission gratings for the Advanced X-ray 
Astrophysics Facility (AXAF) x-ray telescope, cur¬ 
rently scheduled for launch in 1998. Many hun¬ 
dreds of low-distortion, large area transmission 
gratings of 200 nm period (gold) and 600 nm 
period (silver) are required. These will provide 
high resolution x-ray spectroscopy of astrophysical 
sources in the 100 eV to 10 keV band. 

Because of the requirements of low distortion, 
high yield, and manufacturability, a fabrication 
procedure involving the replication of x-ray masks 
has been selected. Masks are made of high¬ 
stiffness silicon nitride membranes to eliminate 
distortion, and patterned using a process involving 
holographic lithography, reactive-ion etching, and 
electroplating. The masks are then replicated 
using soft x-rays (1-1.5 nm) and the resulting pat¬ 
terns electroplated with gold or silver. An etching 
step then yields membrane-supported gratings 
suitable for space use. Flight prototype gratings 
have been fabricated and continue to undergo 
space-worthiness tests. Progress in this area tends 
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Extensive numerical analysis of EGG using mean 
field coarsening theory has been performed. We 
have examined model interface energy functions 
(interface energy versus orientation) and shown 
which of the specific aspects of the interface 
energy function are important in promoting 
bimodal secondary or epitaxial grain growth, i.e., 
grain growth which is strongly orientation selec¬ 
tive. 

4.16 GaAs Epitaxy on 
Sawtooth-patterned Si 

Sponsor 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
Contract DAAL03-92-C-0001 
Spire Corporation 

Project Staff 

Professor Henry I. Smith. Kenneth Yee. Or. Khalid 
Ismail. Nasser Karam> 

The growth of GaAs on Si offers the possibility of 
combining high-speed and optoelectronic GaAs 
devices with Si integrated-circuit technology. Ori¬ 
ented gratings of 200 nm period are fabricated in 
Si 3 N 4 on (100) Si substrates. Anistrotropic 
etching in KOH is then used to produce sawtooth- 
profile gratings in the Si. These then serve as sub¬ 
strates for GaAs growth by MOCVD. The 
dislocation density in the grown GaAs films is 
orders of magnitude lower than the density in films 
formed on planar Si substrates. 
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1.1 Additive Pulse Modelocking 
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Farzana I. Khatri, John D. Moores, Stephen Y. Hon 

Additive Pulse Modelocking (APM) is a technique 
of mode locking that simulates a fast saturable 
absorber via nonlinear interferometric action using 
a Kerr medium. Often this medium is a fiber 
placed in an auxiliary cavity, but the Kerr effect 
may be produced in the components internal to 
the main cavity of the laser. The Kerr Lens Mode¬ 


locking pioneered by Coherent produces ^uiv- 
alent saturable absorber action by self focusing. A 
general theory has been developed for both kinds 
of systems.’ Its main results are: 

1. An analytic treatment of the modelocking 
process in the presence of self-phase modu¬ 
lation and group velocity dispersion; 

2. Criteria for the stability of the modelocked 
pulses; 

3. The condition for self starting of the pulses 
and the explanation of a power threshold for 
mode locking.^ 

The theory is helpful in pointing the way to the 
generation of shorter pulses with lower fluctu¬ 
ations. 

In conjunction with this theory, an analysis of 
noise in modelocked systems hes been carried 
out.’ This analysis provides an explanation of the 
observed spectra of the detector current illumi- 


’ H.A. Haus, J.G. Fujimoto, and E.P. Ippen, 'Structures for Additive Pulse Modelocking,* J. Opt. Soc. Am. B 8: 
2068-2076 (1991). 

2 H.A Haus and E.P. Ippen, 'Self Starting of Passively Modelocked Lasers,' Opt Lett 16:1331-1333 (1991). 

’ S. Hon, K. Bergman, A. Mecozzi, and H.A. Haus, 'Noise Spectra of Modelocked Laser Pulses,' submitted for pres¬ 
entation St CLEO '92, Anaheim, CA (1992). 
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nated by modelocked pulses observed by many 
researchers. Stephen Hon has performed extensive 
measurements of three different laser systems 
(actively modelocked NdiYAG, soliton laser, APM 
TirSapphire laser). We have shown theoretically 
that the current spectrum reveals not only the 
amplitude and timing jitter, as heretofore assumed, 
but also gives information on the correlation 
between the two. These experiments have been 
an important guide to the theory. Conversely, the 
newly developed theory is now being applied to 
extract full information from the detector current 
spectra. The fluctuations that are not observable 
from the current spectrum of direct detection are 
the phase fluctuations (that give rise to the well- 
known Schawlow-Townes linewidth) and the 
carrier frequency fluctuations. 

The analysis has also been extended to allow for 
the inclusion of slow, as well as fast, absorber 
effects. Such a combination may be used to 
advantage in lasers where pulse initiation is a 
problem. There is a trade-off in the range of gain 
and dispersion over which steady-state stability 
can be achieved. 


1.2 Ultrashort Pulse Fiber Laser 

Sponsors 

Charles S. Draper Laboratories 
Contract DL-H-441629 
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Project Staff 
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Ippen, Janice M. Huxley, Lynn E. Nelson, Kohichi 
R. Tamura 

High bit-rate optical communication networks and 
ultrafast optical signal processing systems will 
require compact, chirp-free picosecond pulse 
sources. Erbium-doped fiber lasers are prime can¬ 
didates because of their compatibility with trans¬ 
mission fibers, their operation at 1.5 nm. and their 
amenability to diode pumping. The nonlinear 
refractive index effects in these fibers also make 
them attractive lasers for additive pulse mode¬ 
locking (APM). We have therefore constructed an 
erbium fiber laser with intent to study how APM 


may be implemented in an all-fiber device and to 
optimize ultrashort pulse generation. 

In a fiber, the nonlinear interferometric action nec¬ 
essary for APM is most easily achieved through 
intensity dependent polarization rotation. Thus our 
laser uses fiber polarizers and polarization control¬ 
lers in a fiber ring with an erbium amplifying 
section. Fiber couplers to the ring facilitate 
pumping (with a Ti;sapphire laser for the time 
being) and output coupling. Short pulse operation 
has already been obtained and studied under 
several operating conditions. Pulse spectra indi¬ 
cate durations as short as 300 fs. Autocorrelation 
measurements and variations with pumping level 
indicate that the output consists of strings of 
pulses quantized by soliton effects. A computer 
model has been developed that extends our APM 
theory to cases like this in which pulse shaping 
effects per round trip are large and cannot be 
treated perturbationally. This model is now being 
used to map regions of stability and single pulse 
operation. 

1.3 Long Distance Fiber 
Communications 

Sponsor 

MIT Lincoln Laboratory 
Project Staff 

Professor Hermann A. Haus, Dr. Antonio Mecozzi, 
Yinchieh Lai, John D. Moores 

In 1986, J.P. Gordon and Professor Hermann A. 
Haus published a paper* pointing out that 
transatlantic transmission using solitons would 
suffer from timing jitter of the pulses caused by 
frequency jitter induced by the spontaneous emis¬ 
sion of the amplifier stages. This effect, now 
called the Gordon-Haus effect, has had an impor¬ 
tant influence on the planning of the next 
transatlantic cable of ATErT which will be laid in 
1995. This cable will not use solitons, but will 
operate at (or near) the zero dispersion wave¬ 
length of the fiber. It is designed for 5 Gbit trans¬ 
mission and so that future Wavelength Division 
Multiplexing will not be permitted, at least not 
when using similar channels shifted in frequency. 
The reason for this design is that the dynamic 
range of the input pulse stream of the soliton 
version of the system is limited by two effects. On 
one hand, the pulses have to be energetic enough 


* J.P. Gordon and H.A. Haus, 'Random Walk of Coherently Amplified Solitons in Optical Fiber Transmission,* Opt. 
ierMI; 665-667 (1986). 
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Figure 2. 


by keeping its injection current relatively constant. 
Such structures have achieved wide tuning ranges, 
but not in a continuous manner. 

We are pursuing an alternate approach that relies 
on the ability of controlling the injection currents 
into two quantum wells, one on top of the other, 
by contacts properly fabricated by ion implanta¬ 
tion. The two quantum wells are different and 
have transitions at different frequencies. The 
respective gain peaks are displaced in frequency, 
and, along with them, the index characteristics. 
By controlling the currents injected into the two 
quantum wells separately, one can control the gain 
and index independently. Such structures would 
be ideal for tuning of monolythic resonators at 
constant gain. 


1.7 Gain Dynamics in 
Semiconductor Amplifiers 

Sponsors 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
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National Science Foundation 
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U.S. Air Force - Office of Scientific Research 
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Project Staff 
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Lenz, Dr. Kazunori Naganuma 

The potential application of semiconductor ampli¬ 
fiers to wideband communication networks and 
integrated photonic circuits motivates our studies 
of these devices. We use femtosecond optical 
pulses to probe propagation effects in amplifier 
structures, characterize optical nonlinearities in the 
gain, and observe fundamental ultrafast carrier 
dynamics in the semiconductor active layers. 

During the past year we have performed the first 
femtosecond time-domain measurements of group 
velocity dispersion (GVD) in diode lasers at 
1.5 /im.® The effects of GVD are two-fold: (1) it 
can distort pulses being transmitted and amplified, 
and (2) it limits the temporal resolution, or fre¬ 
quency bandwidth, of nonlinear optical inter¬ 
actions. Our measurements, in both bulk and 
multiple-quantum-well (MOW) samples yield 
values for >ld*n,fl/d2* In the range — 0.63 to 
0.95 ^m~'' and indicate that intrinsic material 
dispersion is the dominant factor. In addition, we 
discovered differences in group velocity between 
TE and TM modes on the order of 200 fs/mm. 
These must be accounted for in studies of aniso¬ 
tropic nonlinearities. 

We have also demonstrated a new technique for 
femtosecond pump-probe measurement that 
makes it possible to use parallel-polarized pump 
and probe beams even in strictly collinear geom¬ 
etries.'*® The method makes use of acousto-optic 
frequency shifting and relies on heterodyne 
detection to distinguish between pump and probe 
pulses. In experiments with bulk (V-groove) 
diode amplifiers, we have demonstrated the obser- 


• K.L. Hall, G. Lenz, and E.P Ippen, 'Femtosecond Time Domain Measurements of Group Velocity Dispersion in 
Diode Lasers at 1.5 f/m,' J. Lightwave Tech., forthcoming. May 1992. 

*0 K.L. Hall, G. Lenz, E.P. Ippen, and G. Raybon, 'Heterodyne Pump-Probe Technique for Time-Domain Studies of 
Optical Non linearities in Waveguides,' submitted to Opt. Lett. 
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vation of differences due to nonlinear as well as 
linear anisotropy. Preliminary results have been 
obtained on strained-layer quantum-well devices 
which have highly anisotropic gain. 

1.8 Uitrafast Optical Kerr 
Effect in Active Waveguides 

Sponsors 

National Center for Integrated Photonics 
National Science Foundation 
Grant EET 88-15834 

U.S. Air Force - Office of Scientific Research 
Contract F49620-91 -C-0091 

Project Staff 

Professor Erich P. Ippen, David J. Dougherty. 
Charles T. Hultgren 

Active optical waveguides offer interesting possi¬ 
bilities for performing nonlinear functions in pho¬ 
tonic circuits. By varying the injected carrier 
density one can change not only the linear compo¬ 
nents of gain and index of refraction but also the 
nonlinear ones. A particularly interesting regime of 
operation is that near the point of nonlinear trans¬ 
parency. At that point, an optical pulse traveling 
through the waveguide neither creates carriers by 
absorption nor removes carriers by stimulated 
emission. There is no long-lived change. Yet the 
optical wavelength is above bandgap and the 
density of carriers in the medium is high. 

In previous studies of gain dynamics in AIGaAs 
lasers near this point, we discovered a large tran¬ 
sient gain nonlinearity which we identified with 
nonequilibrium carrier heating. The initial purpose 
of the work described here was to investigate the 
index of refraction dynamic associated with this 
gain nonlinearity. We have done so.” using a 
time-division interferometric technique originally 
demonstrated by Professor Fujimoto's group in 
passive guides. With this technique we have 
observed the predicted index dynamic in 
channeled-substrate planar (CSP) AIGaAs diode 
lasers (Hitachi HLP 1400) biased just below 
threshold. The dynamic is identified by its approx¬ 


imately 1.7 ps recovery time and positive sign that 
is consistent with a theoretical calculation for non¬ 
equilibrium heating. To our surprise, however, 
these experiments have also uncovered a faster 
and even larger index change. We attribute this 
new component to an above-band optical Stark 
shift.’* 

Our experiments were performed with orthogonally 
polarized pump and probe, so that what we 
measure is 021 . The ultrafast Stark-like component 
has a negative sign, as do below-bandgap values, 
but an enhanced value of about 
- 5 X 10"’*cm*/W. We believe that the resonant 
enhancement in the active region has contributions 
from filled states below the laser photon energy as 
well as from filled states, above. A theoretical 
description of this possibility has now been put 
forward by Sheik-Bahai.’* Such an effect may 
have application to ultrafast optical switching. It 
can also provide a new means for studying femto¬ 
second carrier relaxation phenomena in semicon¬ 
ductors. 

1.9 Coherent Phonons in 
Electronic Materials 

Sponsors 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
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U.S. Air Force - Office of Scientific Research 
Contract F49620-91 -C-OOSI 

Project Staff 

Professor Erich P. Ippen. Tak K. Cheng, Siegfried 
B. Fleischer 

With optical pulses shorter than 100 fs it has 
become possible to excite and observe in the time 
domain the coherent oscillation of optical phonons 
with frequencies up to 10 THz (333 cm~’). In col¬ 
laboration with Professor Dresselhaus' group we 
have been investigating the application of this 
possibility to electronic materials. Recently, we 
have reported that excitation of such phonons is 


C.T. Hultgren and E.P. Ippen, 'Ultrafast Refractive Index Dynamics in AIGaAs Diode Laser Amplifiers,' Appl. Phys. 
Lett. 59:635 (1991). 

C.T. Hultgren and E.P. Ippen, 'Ultrafast Optical Kerr Effect in Active Semiconductor Waveguides,' Proceedings of 
the OS.A. Topical Meeting on Nonlinear Guided-Wave Phenomena, paper TuB2, Cambridge, England, September 
1991. 

r* M. Sheik-Bahai and E.W. Van Stryland, 'Theory of Ultrafast Nonlinear Refraction in Semiconductor Lasers,' Pro¬ 
ceedings of the O.SA. Annual Meeting, San Jose, California, November 1991, paperTuB4, post-deadline paper. 


86 RLE Progress Report Number 134 






Chapter 1. Optics and Quantum Electronics 


particularly efficient in certain semimetals and 
semiconductors (e.g., Bi, Sb, Te and Ti 203 ).'* 

For our experiments we used pulses from a 
dispersion-balanced colliding-pulse-modelocked 
(CPM) dye laser, focused to a 2 pm spot. 
Reflectivity oscillations were observed in these 
materials with amplitudes of several percent. That 
is several orders of magnitude larger than can be 
expected from impulsive stimulated Raman scat¬ 
tering. In each material only one oscillation fre¬ 
quency was observed, the one corresponding to 
the Aig mode; and the phase of the amplitude was 
cosinusoidal with respect to excitation at t=0. Ail 
of these facts are inconsistent with a stimulated 
Raman mechanism. To explain them, we have 
proposed a new theory, dubbed displacive excita¬ 
tion of coherent phonons (DECP).’® It attributes 
behavior to a rapid displacement of the ion equi¬ 
librium coordinates by the electronic excitation. 

The temperature dependence of this effect in Ti 203 
is particularly interesting. In the range 300 K-570 
K dramatic changes take place in Ti 203 with 
respect to both its lattice and electronic properties, 
turning it from a semiconductor to a metal. In our 
femtosecond experiments, we have observed tem¬ 
perature dependent changes in the optical phonon 
frequency and electronic relaxation rate that are 
consistent with this transition. On a dynamic basis 
this may mean that the bandgap and conductivity 
of the material are being modulated at the 6 THz 
phonon frequency. Work in progress will test this 
hypothesis. Higher excitation levels may also 
make it possible to observe enharmonic lattice 
dynamics. In recent preliminary experiments with 
an amplified CPM system we have measured 
refractive index modulations at the phonon fre¬ 
quency of as much as 15 percent. 


1.10 Femtosecond Studies of 
Superconductivity 

Sponsors 

Joint Sen/ices Electronics Program 
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Femtosecond optical experiments that monitor the 
rate of hot electron cooling in materials provide an 
estimate of the strength of the electron-phonon 
interaction. If one is dealing with BCS-like super¬ 
conductors, and one knows the electron-phonon 
spectral density, femtosecond data on the cooling 
rate can be related directly to Tc . We have been 
successful in making such determinations in a 
range of metallic superconductors and in some 
high-Te materials.''® Recently, we have performed 
the first femtosecond pump-probe experiments on 
Ceo and two of its intercalants. 

Most of the properties of Ceo and its intercalants 
can be explained by treating them as rolled up 
sheets of graphite. For this reason it is surprising 
that KaCeo has a superconducting transition tem¬ 
perature (18 K) more than an order of magnitude 
higher than that of the other graphite intercalation 
compounds. Determining the reason for this dif¬ 
ference may help develop a better understanding 
of the mechanisms leading to high-Tc. 

For our experiments, thin films (lOOOA) of 
Ceo, KsCeo. and KeCeo were deposited on quartz 
substrates and encapsulated to prevent contam¬ 
ination. Femtosecond dynamics were induced and 
observed, in reflectivity and transmission at 
X = 630 nm, using pulses from a colliding-pulse- 
modelocked (CPM) dye laser. We believe we are 
monitoring carrier cooling rates and have deduced 
time constants of 9 ps, 0.2 ps and 0.7 ps for the 


T.K. Cheng, J. Vidal, H.J. Zeiger, G. Dresselhaus, M.S. Dresselhaus, and E.P. Ippen, ‘Mechanism for Displacive 
Excitation of Coherent Phonons in Sb, Bi, Te, and Ti 203 ,* Appl. Phys. Lett. 59:1923 (1991). 

'® H.J. Zeiger, J. Vidal, T.K. Cheng, E.P. Ippen, G. Dresselhaus, and M.S. Dresselhaus, ‘Theory for Displacive Excita¬ 
tion of Coherent Phonons,' Phys. Rev. B 45: 768 (1992). 

S. D. Brorson, A. Kazeroonian, J.S. Moodera, D.W. Face, T.K. Cheng, E.P. Ippen, M.S. Dresselhaus, and G. 
Dresselhaus, 'Femtosecond Room-temperature Measurement of the Electron-phonon Coupling Constant 2 in 
Metallic Superconductors," Phys. Rev. Lett. 64: 2172 (1990); S.D. Brorson, A. Kazeroonian, D.W. Face, G.L. Doll, 

T. K. Cheng, M.S. Dresselhaus, G. Dresselhaus, E.P. Ippen, T. Venkatesan, X.D. Wu, and A. Inam. ‘Femtosecond 
Thermomodulation Study of high-Tc ,' So/id. State Commun. 74:1305 (1990). 
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three samptes, respectively. As might be expected, 
the superconducting sample has the shortest time 
constant. But the absolute value of the rate, in the 
context of the material's phonon spectral density, 
is surprisingly low for a Tc of 18 K. One possible 
explanation is that electron-phonon coupling is 
specific to certain types of phonons. This will be 
investigated by experiments with different material 
compositions as well as by more detailed studies 
of these particular intercalants. 

1.11 Femtosecond Pulse 
Generation in Solid State 
Lasers 

Sponsors 
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1.11.1 Introduction 

Ultrashort pulse laser technology is of primary 
importance for the development of high speed 
optical communication and measurement tech¬ 
niques and for studies of ultrafast phenomena. 
The principle objectives of our research program 
are the investigation and development of new 
ultrashort laser pulse generation and measurement 
techniques. Special emphasis is placed on studies 


of solid state ultrashort pulse laser technology. 
Solid state lasers represent a third generation tech¬ 
nology that has several advantages over existing 
dye-laser technology. Solid state lasers feature 
lower cost, greater compactness, higher power, 
and broader tunability. The large gain bandwidths 
of solid state laser media can support pulse dura¬ 
tions in the femtosecond range with wavelength 
tunability. These features make solid state laser 
technology a versatile tool for engineering and 
research applications in ultrafast phenomena. 

One of the most promising solid state laser mate¬ 
rials that has emerged to date is TiiA^Os-'^ The 
TiiAljOs laser operates at room temperature, has a 
high energy storage suitable for high power opera¬ 
tion, and a broad gain bandwidth with an 
extremely wide tuning range from 670 to 1000 nm. 
These features make it an ideal system for femto¬ 
second pulse generation and spectroscopy. 

Recently, we have developed and demonstrated a 
new ultrashort pulse generation technique known 
as Additive Pulse Modelocking (APM).'* Self 
starting, passive modelocking was demonstrated in 
Ti:Al 203 and pulse durations of 200 fs were 
achieved.'* An APM laser utilizes an external 
cavity containing a nonlinear element such as an 
optical fiber to produce a pulse shortening mech¬ 
anism similar to saturable absorber action. A 
saturable absorber is a nonlinear element that exhi¬ 
bits greater transmittivity or reflectivity for higher 
intensities. This nonlinearity produces pulse short¬ 
ening by enhancing the peak of the pulse com¬ 
pared to the wings. 

Our research in modelocked Ti;Al 203 represented 
the first demonstration of self starting passive 
Additive Pulse Modelocking. The self-starting 
APM system represented a significant advance¬ 
ment because it eliminated the need for a high fre¬ 
quency loss modulation element and thus reduced 
system complexity and cost. We have extended 
the self-starting APM technique to modelock 
diode pumped laser systems and have generated 
1.7 ps pulses in a Nd:YAG system and 2 ps pulses 
in a Nd:YLF system.*® 


P.F. Moulton, "Spectroscopic and Laser Characteristics of Ti:Al 203 ," J. Opt. Soc. Am. 9 3:125 (1986). 

'• E.P. Ippen, 'Additive Pulse Modelocking," J. Opt. Soc. Am. B 6; 1736-1745 (1989); J. Goodberlet, J. Wang, J.G. 
Fujimoto, and P.A. Schulz, "Femtosecond Passive Modelocked Ti;Al 203 Laser with a Nonlinear External Cavity," 
Opt. Lett. 14:1125(1989). 

J. Goodberlet, J. Wang, J.G. Fujimoto, and P.A. Schulz, "Femtosecond Passive Modelocked Ti:Al 203 Laser with a 
Nonlinear External Cavity," Opt. Lett. 14; 1125 (1989). 

w J. Goodberlet, J. Jacobson, J.G. Fujimoto, P.A. Schulz, and T.Y. Fan, 'Self-starting Additive-pulse Mode-locked 
Diode-pumped Nd:YAG Laser," Opt. Lett. 15: 504-506 (1990). 
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1.11.2 Femtosecond Generation 
Using Nonlinear Intracavity 
Elements 

Studies by other investigators have recently dem¬ 
onstrated that a single cavity Ti:Al 203 system can 
be modelocked without an active modulator, 
saturable absorber, or nonlinear external cavity.^' 
This phenomena was termed self-modelocking. 
Although these investigations did not elucidate the 
short pulse generation mechanisms, they demon¬ 
strated that a solid state laser could be mode- 
locked solely by intracavity optical nonlinearities. 

Solid state femtosecond lasers such as Ti:Al 203 
can achieve peak intracavity intensities of 0.5 MW 
when they are modelocked. These high peak 
intensities make possible a wide range of new 
modelocking techniques which utilize intracavity 
nonlinearities. Saturable absorber action based on 
all-optical solid state refractive index nonlinearities 
is of particular interest because of its extremely fast 
response and recovery times. Our group has 
pursued this approach further and demonstrated 
two different types of passive modelocking 
schemes in Ti:Al 203 using all solid state nonlinear 
intracavity modulators.^^ These modelocking tech¬ 
niques permit the generation of femtosecond pulse 
durations and can be extended to a wide range of 
solid state lasers at different wavelengths. In addi¬ 
tion, they represent a significant reduction in cost 
and system complexity over previous modelocktng 
techniques. 

One approach for modelocking uses a nonlinear 
intracavity element which functions by nonlinear 
polarization rotation. This device is called a Kerr 
Polarization Modulator (KPM)^^ and consists of a 
quarter-wave plate, a focusing lens, an SF11 glass 
plate, a second quarter-wave plate, and a curved 


end mirror. The KPM works by generating an 
elliptical polarization state which undergoes a 
nonlinear rotation due to the intensity dependent 
refractive index of the SF11 glass plate. The linear 
and nonlinear transmission functions of the KPM 
can be controlled by setting the retardations of the 
quarter-wave plates. Thus the KPM can be made 
to have an increasing transmission with increasing 
intensity and function as a saturable absorber. We 
have demonstrated this technique for modelocking 
Ti:Al 203 and achieved pulse durations of 230 fs. 

Another approach for modelocking relies on the 
use of intracavity self focusing. This technique 
has been termed Kerr Lens Modelocking (KLM),^* 
and pulse durations in the 100 fs range have been 
demonstrated in Ti:Al 203 . KLM produces a fast 
saturable-like action by using self focusing in the 
laser rod combined with intracavity aperaturing. 
Research in our group has emphasized developing 
both a fundamental understanding of the under¬ 
lying nonlinear processes as well as an approach 
for engineering these systems. We have intro¬ 
duced a significant simplification in the design of 
KLM lasers by developing a modular nonlinear 
intracavity element called a micro-dot mirror.^s 

The end mirror of the laser is substituted with a 
focusing lens and a mirror with small diameter 
micro-dots. The mirror is antireflection (AR) 
coated on the first surface (i.e., the surface that 
laser beam first encounters) and high reflection 
(HR) coated on the second surface. The second 
surface is patterned into an array of HR spots 
(micro-dots) with varying diameter using micro¬ 
electronic techniques. Fast saturable absorber 
action is produced by nonlinear self focusing 
which produces a smaller spot size for higher 
intensities. When a pulse is reflected from the 
micro-dot mirror, the higher intensity peaks of the 
pulse are reflected more than the lower intensity 


21 O.E. Spence, P.N. Kean, and W. Sibbett, '60-fsec Pulse Generation from a Self-mode-locked Ti.sapphire Laser,* 
Opt. Lett. 16: 42 (1991). 

22 G. Gabetta, 0. Huang, J. Jacobson, M. Ramaswamy, H.A. Haus, E.P. Ippen, and J.G. Fujimoto, ‘Femtosecond 
Pulse Generation in Ti.AlaOs Using Nonlinear Intracavity Elements,* CLEO 1991, Baltimore, Maryland, 1991; G. 
Gabetta, 0. Huang, J. Jacobson, M. Ramaswamy, H.A. Haus, E.P. Ippen, and J.G. Fujimoto, ‘Femtosecond Pulse 
Generation in Ti:Al 203 Using a Microdot Mirror Modelocker,* Opt. Lett. 16:1756 (1991). 

23 G. Gabetta, D. Huang, J. Jacobson, M. Ramaswamy, H.A. Haus, E.P. Ippen, and J.G. Fujimoto, ‘Femtosecond 
Pulse Generation in Ti:Al 203 Using Nonlinear Intracavity Elements,* CLEO 1991, Baltimore, Maryland, 1991. 

24 D.K. Negus, L. Spinelli, N. Goldblatt, and G. Feugnet, ‘Sub-100 Femtosecond Pulse Generation by Kerr Lens 
Modelocking in Ti:Al 203 ,' OSA Meeting on Advanced Solid State Lasers, Hilton Head, South Carolina, March 
18-20, 1991, postdeadline paper. 

26 G. Gabetta, D. Huang, J. Jacobson, M. Ramaswamy, H.A. Haus, E.P. Ippen, and J.G. Fujimoto, ‘Femtosecond 
Pulse Generation in Ti:Al 203 Using a Microdot Mirror Modelocker,* Opt. Lett. 16: 1756 (1991); D. Huang, M. 
Uiman, L. Acioli, H.A. Haus, and J.G. Fujimoto, ‘Self-focusing Induced Saturable Absorber Loss for Laser Mode 
Locking," Opt. Lett., forthcoming. 
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wings of the pulse. Thus the micro-dot mirror 
effectively acts as a fast saturable absorber. We 
have demonstrated the micro-dot mirror for mode¬ 
locking the Ti:Al 203 laser and achieved pulse dura¬ 
tions of 190 fs. 

The modelocking techniques described above have 
led to a better understanding of passive mode- 
locking mechanisms in solid state laser systems. 
In the case of self-focusing modelocking, we 
developed a new theoretical framework which is 
applicable to the majority of single cavity passively 
modelocked Ti:Al 203 systems.^^ Theoretical studies 
are important not only for understanding pulse 
shortening mechanisms but also for investigating 
factors which limit pulse duration. 

We have developed an analytical description of 
modelocking which includes fast saturable absorp¬ 
tion, self phase modulation, gain bandwidth lim¬ 
iting, and dispersion. A pulse in a modelocked 
laser can be mathematically described by a single 
equation similar to the nonlinear Schrodinger 
equation. We have studied the behavior of pulse 
characteristics such as pulse duration, chirp, band¬ 
width, and stability.^B Our primary findings show 
that there are two basic operating regimes where 
stable short pulse generation may be achieved. A 
laser may be constructed either with positive 
(normal) or negative (anomalous) dispersion. The 
pulse shaping mechanisms for these two regimes 
are different and a laser that has anomalous 
dispersion will in general produce dramatically 
shorter pulses than one with normal dispersion. 
These results suggest design criteria for a wide 
range of modelocked solid state lasers. 


1.11.3 Reduction of Pulse Duration 
by Dispersion Compensation 

One of the key topics of investigation in ultrashort 
pulse generation is the development of approaches 
for generating the shortest possible pulse dura¬ 
tions. The fundamental limit to pulse duration is 
determined by the energy-time uncertainty prin¬ 
ciple which implies that the pulse duration is 
inversely related to the laser gain bandwidth. The 
gain bandwidth in TirA^Os is in excess of 100-200 
nm27 and thus pulse durations of a few femtosec¬ 
onds are theoretically possible. Pulse compression 
techniques using amplified femtosecond pulses 
have been demonstrated to achieve pulses as short 
as 6 fs.’B However, pulses in the 10 fs range have 
not as yet been generated directly from a laser. 
The broad bandwidths associated with these pulse 
durations allow the use of frequency synthesis 
techniques to control pulse shape or generate mul¬ 
tiple pulse, multiple wavelength pulse trains.^* 
Thus, the development of new sources with pulse 
durations in the 10 fs range would have a major 
impact on the study of ultrafast phenomena. 

One area of investigation in our program is the 
exploration of the factors which limit the short 
pulse duration that can be generated. The theore¬ 
tical considerations discussed above are valid only 
in systems where the pulses are long enough that 
dispersive terms higher than those of second order 
can be neglected. Previous pulse compression 
studies have demonstrated that higher order 
dispersion presents an important limitation to pulse 
duration as pulses become very short and their 
corresponding frequency bandwidths become 
broad.s® 

In order to explore and extend the limits of short 
pulse generation, we have developed and demon¬ 
strated a solid state laser with independently 
adjustable second and third order dispersion. 
Compensation of dispersion is achieved using a 
pair of intracavity prisms which are set to produce 
a net negative second order dispersion in the 


2B H.A. Haus, J.G. Fujimoto, and E.P. Ippen, 'Structures for Additive Pulse Modelocking,' J. Opt. Soc. Am. B 8: 
2068 (1991). 

27 p.F. Moulton, 'Spectroscopic and Laser Characteristics of TiiAbOs,* J. Opt. Soc. Am. 0 3:125 (1986). 

2S R.L. Fork, C.H. Brito Cruz, P.C. Becker, and C.V. Shank, 'Compression of Optical Pulses to Six Femtoseconds by 
Using Cubic Phase Compensation,' Opt. Lett. 12; 483 (1987). 

M A.M. Weiner, J.P. Heritage, and E.M. Kirschner, 'High-resolution Femtosecond Pulse Shaping,’ J. Opt. Soc. Am. 
B 5:1563 (1988). 

W.J. Tomlinson and W.H. Knox, 'Limits of Fiber-grating Optical Pulse Compression,' J. Opt Soc. Am. B 4:1404 
(1987): C.H. Brito Cruz, P.C. Becker, R.L. Fork, and C.V. Shank, 'Phase Correction of Femtosecond Optical Pulses 
Using a Combination of Prisms and Gratings," Opt Lett. 13:123 (1988). 
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cavity. Intracavity prisms are a standard technique 
for dispersion compensation which have been 
applied to a number of solid state and dye lasers. 
An undesirable consequence of this prism arrange¬ 
ment, however, is that it also introduces a signif¬ 
icant negative third order dispersive term which 
produces pulse distortion effects.^ In order to 
explicitly compensate for third order dispersion, we 
have developed an intracavity thin film element 
which has positive third order dispersion but zero 
second order dispersion. Working in collaboration 
with investigators at CVI Incorporated, we 
designed and fabricated a thin film Gires-Tourneau 
Intederometer (GTI) on a mirror substrate. The 
amount of third order compensation introduced in 
the laser cavity can be varied by adjustment of the 
number of times the pulse bounces off each GTI 
surface during each round trip. This is the first 
laser that has directly adjustable second and third 
order dispersion. Using this technique we have 
generated pulse durations as short as 28 fs directly 
from a modelocked TrAhOs laser. These pulses 
are the shortest pulses which have been directly 
generated from a Ti:Al 203 laser without using 
pulse compression techniques.^^ 

1.12 Studies of Ultrafast 
Phenomena in Optoelectronic 
Materials 

Sponsors 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
Contract DAAL03-92-C-0001 
National Science Foundation 
Grant ECS 85-52701 

U.S. Air Force - Office of Scientific Research 
Contract F49620-88-C-0089 
U.S. Navy - Office of Naval Research 
Grant N00014-91-J-1956 

Project Staff 

Dr. Lucio H. Acioli, Morrison Ulman, David Huang, 
Professor Erich P. Ippen, Professor James G. 
Fujimoto 


1.12.1 Four Wave Mixing and 
Temporal Encoding in Barium 
Titanate 

The idea of optically encoding information in 
photosensitive materials is relevant to applications 
in optical signal processing, image processing, and 
optical logic. Holography is an example of a tech¬ 
nique in which three-dimensional spatial informa¬ 
tion can be recorded in a suitable matrix such as 
film or photorefractive material. In addition to 
spatial information, temporal information may also 
be encoded. The search for systems to record 
information in the temporal domain have usually 
been based on the concept of encoding informa¬ 
tion using inhomogeneously broadening transi¬ 
tions in atomic systems. The techniques most 
often applied for temporal storage have been vari¬ 
ations of photon echo like schemes. Working in 
collaboration with investigators from Tufts Univer¬ 
sity, we have performed the first studies of femto¬ 
second wave-mixing in BaTiOa and demonstrated 
the possibility of using volume holography for 
recording temporal information.^ 

Photorefractive materials such as BaTiOa, or SBN, 
have been demonstrated for a wide range of 
optical phase conjugation and signal processing 
experiments using continuous wave lasers.^^ These 
materials exhibit extremely large nonlinearities, and 
studies can be performed without the need for 
high laser intensities. However, the response and 
relaxation times of the photorefractive materials are 
extremely slow, and thus only recently has the use 
of such materials under pulsed laser excitation 
been considered. In order to circumvent the slow 
material response, our technique records temporal 
information in the spatial domain using four wave 
mixing and modulated volume photorefractive 
gratings. 

Experiments were performed using a colliding 
pulse modelocked, CPM, dye laser source which 
generated 40 fs pulse durations. Temporal infor¬ 
mation was encoded using a signal and a refer¬ 
ence beam. When the signal and reference beams 
interact in the photorefractive crystal, they form a 


3^ J.M. Jacobson, A.G. Jacobson, K. Naganuma, and J.Q. Fujimoto, 'Generation of 28 fs Pulses from a TiAlgOs 
Laser Using Second and Third Order Intracavity Dispersion Compensation,' paper to be presented at CLEO 1992, 
Anaheim, California. 

*2 L.H. Acioli, M. Ulman, E.P. Ippen, J.G. Fujimoto, H. Kong, B.S. Chen, and M. Cronin-Golomb, 'Femtosecond 
Temporal Encoding in Barium Titanate,' Opt. Lett. 16:1984 (1991). 

33 M. Cronin-Golomb, B. Fischer, J.O. White, and A. Yariv. 'Theory and Applications of Four-wave Mixing in Pho¬ 
torefractive Media," /£££ J. Quant Electron. QE-20:12 (1984). 


91 


Chapter 1. Optics and Quantum Electronics 


standing wave interference pattern where they 
geometrically overlap. This in turn induces a pho- 
torefractive diffraction grating. If the signal team 
contains a femtosecond temporal waveform, and 
the reference beam consists of a single femto¬ 
second pulse, a spatially varying interference 
pattern will be formed according to the geomet¬ 
rical overlap of the signal and the reference optical 
waveforms. The temporal information on the 
signal beam is recorded as a volume grating in the 
BaTiOa crystal. The photorefractive gratings have a 
long persistence time and can be used to store the 
temporal information. The femtosecond waveform 
can subsequently be reconstructed by diffracting a 
femtosecond reference pulse from the stored 
volume grating.3^ This process is analogous to 
image reconstruction in holography. 

Different phase conjugation geometries were 
examined including the ring resonator and the two 
beam coupling geometry. A surprising discovery 
was that temporal signals are influenced only by 
material dispersion effects and that pulse durations 
of 40 fs could be preserved. We have determined 
that the factors that influence the temporal resol¬ 
ution of the temporal encoding and reconstruction 
process and developed a simple geometrical model 
to predict resolution and information storage 
capacity. Our experiments are closely related to 
femtosecond holography which uses holographic 
recording to store transient femtosecond images,^^ 
except that our results demonstrate an extension of 
this technique to both the storage and recon¬ 
struction of temporal information. Variations of 
our approach using acousto-optic modulators or 
other programmable volume diffraction devices 
could make possible the generation of program¬ 
mable optical pulses at THz repetition rates. 


1.12.2 Self-Focusing and Nonlinear 
Propagation 

Nonlinear optical media with Kerr-like refractive 
index nonlinearities are relevant to a wide range of 
applications including all-optical switching, optical 
modulation, short pulse generation, and short 
pulse propagation in optical fibers. Self-focusing 
is a classic nonlinear optical process which has 
been the subject of investigation for many years. 
This process has been of interest in previous 
studies because it poses a limitation to high power 
laser generation. Recently, however, self-focusing 
has become technologically important for applica¬ 
tions which involve high speed all-optical modu¬ 
lation and short pulse generation. A new 
technique for modelocking solid state lasers, called 
Kerr lens modelocking has been developed which 
uses self focusing to achieve fast saturable 
absorber action.^B Pulse durations of less than 100 
fs have been achieved in the Ti:Al 203 laser. Self- 
focusing also plays a major role in other short 
pulse generation techniques such as high power 
femtosecond pulse compression and continuum 
generation.*^ 

Self-focusing is produced by the intensity depen¬ 
dent nonlinear index of refraction of the material 
and the spatial intensity variation of the laser 
beam. The transverse spatial intensity variation of 
the beam induces a graded index lens in the non¬ 
linear material which, in turn, produces a 
geometric focusing of the beam. Several 
researchers have studied this problem in the past, 
and, with the advent of high intensity femtosecond 
lasers, a description of the temporal as well as the 
spatial behavior of self-focusing has become of 
interest. There is no theory yet that fully describes 
the problem of nonlinear femtosecond pulse prop¬ 
agation in bulk nonlinear media. 

Our investigations of self-focusing are directed 
toward applications in laser modelocking and short 
pulse generation. Self-focusing has recently been 
used to achieve fast saturable absorber mode- 


M L.H. Acioli, M. Ulman, E.P. Ippen, J.G. Fujimoto, H. Kong, 6.S. Chen, and M. Cronin-Golomb, ‘Femtosecond 
Temporal Encoding in Barium Titanate,' Opt. Lett. 16:1984 (1991). 

30 J.A. Valdmanis, H. Chen, E.N. Leith, Y. Chen, and J L. Lopez, 'Three Dimensional Imaging with Femtosecond 
Optical Pulses,' CLEO Technice/ Digest paper CTUA1,1990, p. 64. 

30 D.K. Negus, L. Spinelli, N. Goldblatt, and G. Feugnet 'Sub-100 Femtosecond Pulse Generation by Kerr Lens 
Modelocking in Ti:Al 203 ,' OSA Meeting on Advanced Solid State Lasers. Hilton Head, South Carolina, March 
18-20,1991, postdeadline paper. 

37 c. Holland and P.B. Corkum, 'Compression of High-power Optical Pulses,' J. Opt. Soc. Am. B S: 641 (1988); 
D.H. Reitze, A.M. Weiner, and O.E. Leaird, 'High-power Femtosecond Optical Pulse Compression by Using 
Spatial Solitons,' Opt. Lett. 16: 1409 (1991); R.L. Fork, C.V. Shank, C. Hirlimann, R. Yen, and W.J. Tomlinson, 
‘Femtosecond White-light Continuum Pulses,’ Opt. Lett. 8:1 (1983). 
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locking in Ti;Al 203 . This process has been termed 
Kerr Lens Modelocking and can form the basis for 
ultrashort pulse generation in a wide range of solid 
state lasers.^ As discussed previously, we have 
developed a modular device called a micro-dot 
mirror which uses self-focusing to modelock the 
Ti:Al 203 laser.^B In addition, self-focusing in bulk 
materials can also be used for short pulse gener¬ 
ation and nonlinear measurement. High power 
pulse compression has been demonstrated using 
self phase modulation in bulk nonlinear materials.^ 
it has also been shown that self-focusing plays a 
fundamental role in the process of continuum gen¬ 
eration with femtosecond pulses and can be used 
to generate spatial optical solitons for pulse com¬ 
pression.^' An experimental technique, based on 
the self-focusing effect and referred to as 
‘z-scan,’42 has been introduced as a sensitive 
probe of the magnitude of the nonlinear refractive 
index. 

In order to investigate the use of self-focusing for 
modelocking and ultrashort pulse propagation, we 
have developed a new approach to describe the 
effects of self-focusing in thick (compared to the 
confocal parameter of the beam) nonlinear 
media.*3 This method is based on a complex 
scaling operation applied to the q parameter of a 
gaussian beam. Self-focusing can be described in 
the perturbative limit using a compact formulation 
that is compatible with the ABCD matrix 
description of gaussian beam propagation. Exper¬ 
imental 'z-scan* measurements of self-focusing 
were performed using a femtosecond CPM (col¬ 


liding pulse modelocked) dye laser and high repe¬ 
tition rate copper vapor laser pumped amplifier.** 
Excellent agreement between the theoretical pred¬ 
ictions and the experimental results have been 
obtained for thick optical glasses (BK7 and fused 
silica). These theoretical and experimental studies 
provide a basis for designing and developing Kerr 
lens modelocking techniques in a wide range of 
solid state lasers. 

We have also studied the effects of self-focusing 
in continuum generation. Our investigations use 
variations of the *z-scan' technique*^ in order to 
quantify the continuum generation process by 
varying the nonlinear propagation parameters of 
the beam. Our current studies focus on describing 
energy transfer to wavelengths different from that 
of the input beam, the amplitude stability, and the 
divergence properties of the output beam. An 
increased understanding of these processes can 
form the basis for developing new short pulse 
generation techniques based on nonlinear propa¬ 
gation. 

1.12.3 Femtosecond Carrier 
Dynamics in Semiconductors 

The efficient design of high speed electronic and 
photonic devices requires a detailed understanding 
of transient carrier dynamics in semiconductor 
materiais. The dynamics of electrons and holes 
determine fundamental limitations to device speed. 
For example, intervalley scattering processes affect 


D.K. Negus, L. Spinelli, N. Goldblatt, and G. Feugnet, *Sub-100 Femtosecond Pulse Generation by Kerr Lens 
Modelocking in Ti;Al 203 .’ OSA Meeting on Advanced Solid State Lasers, Hilton Head, South Carolina, March 
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D.H. Reitze, A.M. Weiner, and D.E. Leaird, 'High-power Femtosecond Optical Pulse Compression by Using 
Spatial Solitons.' Opt. Lett. 16:1409 (1991). 
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the high field transport behavior of GaAs devices.^ 
Carrier dynamics have been studied with a number 
of techniques including luminescence spectro¬ 
scopy and femtosecond absorption saturation.*^ 
Our research program combines state-of-the-art 
experimental and theoretical approaches to 
develop an accurate model for carrier dynamics on 
an ultrashort timescale. The ultrahigh temporal 
resolution of femtosecond lasers permits direct 
measurement of electronic scattering processes 
that occur in less than 100 femtoseconds. Carrier 
scattering events may be examined both exper¬ 
imentally and by supercomputer simulation. As a 
result of these investigations, transient electronic 
behavior in a variety of important semiconductor 
materials and quantum structures may be pre¬ 
dicted. Recently, we have focused our efforts on 
intervalley scattering in AIGaAs as a paradigm for 
investigating femtosecond carrier dynamics. 

Our research program is a collaborative effort with 
Professor C.J. Stanton's solid state theory group at 
the University of Florida. The objective of the the¬ 
oretical component of the program is to develop 
both numerical simulation and analytical methods 
which establish a correspondence between exper¬ 
imental femtosecond transient absorption satu¬ 
ration results and theoretical models. The problem 
of investigating transient carrier dynamics in lll-V 
semiconductors is extremely complicated because 
of the large number of scattering processes. Thus, 
advanced theoretical techniques such as Ensemble 
Monte Carlo theory are necessary to accurately 
interpret experimental results. A full zone, 30 
band, empirical k • p procedure is used to deter¬ 
mine the relevant electronic bandstructure. A 
Monte Carlo simulation of 40,000 electrons and 
holes tracks the time development of the electron 
and hole distribution functions. Finally, the differ¬ 
ential transmission is calculated and compared 
directly to the results of femtosecond pump-probe 
experiments.** In this way, fundamental material 
parameters such as deformation potential con¬ 
stants and carrier-carrier scattering rates may be 


determined. An analytic discretization of the 
Boltzmann equation combined with a rate 
equation approach has also been developed to 
provide intuitive insight into the transient carrier 
relaxation process. This interdisciplinary approach 
has already shown that collisional broadening and 
femtosecond hole redistribution must be included 
in the theory in order to match the experimental 
results. We are now investigating the role of 
carrier diffusion and coulomb enhancement on the 
differential transmission data. 

In the past year, work at M'T has been directed 
toward femtosecond pump-probe absorption satu¬ 
ration studies of the GaM and AIGaAs semicon¬ 
ductor systems. Femtosecond intervalley 
scattering in AIGaAs is investigated using a 
tunable wavelength technique. Systematic vari¬ 
ation of the wavelength and spectral content of 
ultrashort laser pulses permits a selective study of 
different scattering channels, such as scattering to 
the L satellite valley. An amplified colliding pulse 
modelocked (CPM) ring dye laser system provides 
the source for our experiments. Femtosecond 
optical pulses from the CPM are amplified by a 
copper vapor laser pumped dye amplifier to 
microjoule energies. The amplified pulses are 
focused on a jet of ethylene glycol to generate a 
femtosecond white light continuum using non¬ 
linear self phase modulation. A Fourier synthesis 
scheme selects variable wavelengths and pulse 
durations for use in pump-probe experiments.** 
This approach has the advantage that it permits 
studies using wavelength tunable pulses with 
durations as short as 30 fs. The samples used in 
our studies are grown by molecular beam epitaxy 
at MIT Lincoln Laboratory. The mole fraction of Al 
in the ALGai.^As system is chosen to vary the 
bandstructure and isolate optical transitions of 
interest. Scattering to the L valley in AloiGaogAs 
is only allowed for photoexcited carriers with 
energy greater than 1.8 eV. The possibility of scat¬ 
tering by this channel affects the behavior of the 
transient in femtosecond absorption saturation. 


** M.A. Littlejohn, J.R. Hauser, T.H. Glisson, D.K. Ferry, and J.W. Harrison, ‘Alloy Scattering and High Field Trans¬ 
port in Ternary and Quaternary lll-V Semiconductors,' Solid State Electron. 21:107 (1978). 
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'Hot-electron Recombination at Neutral Acceptors in GaAs: A cw Probe of Femtosecond Intervalley Scattering,' 
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Thus, a wavelength tunable study discerns salient 
features of the data in both the time and energy 
dimensions. 

The results already obtained will be extended to 
new materials systems and to quantum confined 
structures. Device design in the ^ture will depend 
on thorough knowledge of the transient electronic 
properties of semiconductor materials. The combi¬ 
nation of advanced femtosecond laser technology 
with sophisticated theoretical methods provides a 
potent tool for discovery of fundamental informa¬ 
tion on ultrafast carrier dynamics. 

1.13 Time Domain Diagnostics 
of Waveguide Devices 
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1.13.1 Tima Domain Interferometry 

Investigations of nonlinear processes such as 
measurements of the nonlinear index of refraction 
in waveguide devices are key to the development 
of high speed all-optical switching devices. In 
particular, the magnitude of the nonlinear index is 
an important parameter for all-optical switching 
and modulation, since it determines the laser pulse 
energy which is necessary to achieve switching 
behavior. Time domain measurements permit both 
the magnitude as well as the dynamics of the non¬ 
linear index of refraction to be characterized. Pre¬ 
vious investigators have used a number of time or 
frequency domain techniques to measure nonlinear 
index including four wave mixing,^ nonlinear 
waveguide couplers,^^ nonlinear Fabry-Perots,^^ 
fringe shift interferometry,®® and Mach-Zehnder 
interferometry.®* In general, direct measurements of 
index nonlinearities are complicated because 
measurement techniques are sensitive to thermal or 
acoustical parasitic signals, have low sensitivity, or 
require deconvolution. 

Recently, we have developed a new technique, 
femtosecond time division interferometry (TDI), 
that permits high sensitivity, direct measurements 
of the nonlinear index 02 in waveguide devices 
with reduced interference from thermal or acous¬ 
tical parasitic signals.®® In this technique, pump 
and probe pulses are coupled into a waveguide 
device and the nonlinear phase shift which is 
induced on the probe by the pump is measured by 
interfering it with a reference pulse. TDI was pre¬ 
viously limited to measurements of 021 , i.e., index 
changes induced on a probe pulse by an orthog¬ 
onally polarized pump. However, new materials 
such as quantum wells, quantum wires, organic 
polymers as well as other anisotropic systems are 
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becoming increasingly important candidates for all 
optical switching devices. A complete character¬ 
ization of n 2 requires the measurement of n 2 |. 

During the past year, we have developed a new 
variation of our original time division interferom¬ 
etry, TOI, technique which permits the complete 
characterization of the nonlinear index.B< By using 
a novel phase modulation technique which breaks 
the degeneracy between the pump and probe and 
actively stabilizes the interferometer, the magnitude 
and femtosecond transient response of n 2 | can be 
measured. Experiments have been performed to 
measure n 2 i and n 2 | in an optical fiber to demon¬ 
strate this new measurement technique. High sen¬ 
sitivities have also been achieved and phase shifts 
as small as 5 milirad can easily be detected. Thus, 
TOI can be applied to studies of nonlinear index 
and all-optical switching in a wide range of wave¬ 
guide devices. 

1.13.2 Time Domain Diagnostics of 
Strained Layer Devices 

Strained layer quantum well devices represent one 
of the most active and technologically promising 
areas of current optoelectronics device research. 
Strained layer materials greatly extend the available 
wavelength range of optoelectronic devices by 
permitting the growth of non-lattice materials, in 
addition, strained layer materials represent a 
bandgap engineering approach to achieve superior 
device performance. By modifying the valence 
band structure, lower threshold current density, 
higher efficiency, smaller chirp, wider modulation 
bandwidth, and improved high temperature perfor¬ 
mance^^ can be achieved over conventional 
devices.®* High power, high efficiency, long life¬ 
time, and low threshold current density®* semicon¬ 


ductor lasers have been achieved in InGaAs 
strained layer devices. With their extended wave¬ 
length range and high output power, 
InGaAs/GaAs strained layer lasers are important for 
applications such as coherent optical communi¬ 
cations, pumping erbium doped fiber amplifier and 
solid state lasers, and high speed optical signal 
processing. 

With the development of new tunable femto¬ 
second laser source of wavelengths between 900 
and 1000 nm, time domain diagnostics can 
become a powerful tool for characterizing 
dispersion, gain dynamics, nonlinear index change, 
and transient nonlinearities. We have recently 
begun a program in collaboration with Groups 67 
and 83 at Lincoln Laboratory to develop frequency 
and time domain diagnostics of strained layer 
diode lasers. The objective of our program is 
twofold; To develop and demonstrate new time 
and frequency domain diagnostics of diode lasers, 
and to examine fundamental processes and device 
design issues which limit high power single fre¬ 
quency performance in strained layer diodes. 

Studies will involve InGaAs/GaAs graded-index 
separate-confinement heterostructure strained layer 
single quantum well laser diodes.®<> Femtosecond 
pump probe measurements of gain dynamics, non¬ 
linear index, and dispersion are being performed. 
Continous wave measurements such as linewidth, 
spectral character, self heterodyne four wave 
mixing, and injection locking are being performed 
at Lincoln Laboratory. These two approaches rep¬ 
resent complementary techniques which, taken 
together, can permit a more comprehensive char¬ 
acterization of device and materials properties. 
When completed, this study should establish new 
approaches for performing time domain diagnos¬ 
tics on a wide range of devices as well as provide 
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information relevant to optimizing high power 
device performance in the InGaAs/GaAs system. 

1.14 Laser Medicine and 
Surgery 

Sponsors 
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1.14.1 The Ultrashort Pulse Laser 
Scalpel 

In recent years, the use of short optical pulses has 
enhanced the degree of control and localizability 
of intraocular laser surgery.^' Laser induced optical 
breakdown is a nonlinear laser tissue interaction 
which permits the photodisruption or cutting of 
transparent structures within the eye without the 
need for intervening surgical incision. The objec¬ 
tive of our program is to develop an optimized 
ultrashort pulse laser scalpel which can be used to 
perform surgical incisions of intraocular structures 
even in close proximity to sensitive ocular struc¬ 
tures such as the retina and cornea. These studies 
are part of an onging collaboration between inves¬ 
tigators at MIT, the New England Eye Center of 
New England Medical Center Hospitals, and the 
Wellman Laboratories of Photomedicine at Mass¬ 
achusetts General Hospital. 

The physical processes which occur in laser 
induced breakdown include plasma formation. 


acoustic wave generation, and cavitation. These 
processes produce the laser tissue surgical effects 
as well as collateral damage. To date, the majority 
of clinical applications of laser induced optical 
breakdown have utilized nanosecond duration 
pulses in the millijoule energy range and single 
pulse exposures. We have studied and compared 
the mechanisms, scaling behavior, and tissue 
effects of single pulses in the nanosecond and 
picosecond ranges.*^ in general, nanosecond and 
picosecond optical breakdown results in compa¬ 
rable damage zones if the same amount of energy 
is deposited; however, since the threshold energy 
for breakdown is much lower for picosecond 
pulses, near-threshold picosecond pulses produce 
greatly reduced collateral damage zones. For 
example, we have demonstrated collateral damage 
ranges in a corneal endothelial cell model of only 
100 fitn with 40 picosecond duration pulses at 8 
microjoules pulse energy. We have also studied 
tissue effects (corneal excisions) into the femto¬ 
second domain.*^ Ultrashort pulses with high peak 
intensities can produce a plasma mediated ablation 
of transparent tissues such as the cornea. 
Picosecond and femtosecond pulse durations have 
been demonstrated to produce much smoother 
excision edges and less damage to the adjacent 
tissue than nanosecond pulses. 

We have demonstrated that clinically viable sur¬ 
gical incisions can be made using multiple pulse, 
high repetition rate picosecond pulses. Each pulse 
produces minimal collateral damage while multiple 
exposures produce a cumulative incision effect. In 
order to study the effect of repetition rate, we have 
developed a new NdiYAG laser system which 
produces 100 picosecond duration pulses with a 
maximum energy of 140 //J and a repetition rate 
variable from 3 to 1000 Hz. This laser has been 
used to demonstrate in vitro cutting of a 
monolayer of cultured fibroblast cells suspended 
over corneal endothelium in saline as a model for 
vitreous membrane surgery. Membrane cutting in 
the deep vitreous in close proximity to the retina 
remains one of the most challenging problems for 
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1.1 Introduction 
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Our objectives are to develop, using renormaliza¬ 
tion-group theory and other methods of statistical 
mechanics, microscopic theories of quantum spin 
and electronic systems. Our approach is partic¬ 
ularly suited to systems with fluctuations due to 
finite temperatures, impurities, surfaces or other 
geometric constraints. 

1.2 Renormalization-Group 
Approach to Electronic Systems 

Project Staff 

Alexis Falicov, Professor A. Nihat Barker 

High-Tc superconductivity, metallic magnetism, the 
metal-insulator transition, and heavy fermion 
behavior are all phenomena produced by the 
strong correlations of electrons in narrow energy 
bands. It is therefore important to study theore¬ 
tical models that incorporate the strong correlation 
effects of electrons. The tJ model is such a 
system. It is defined on a lattice with one 
spherically symmetric orbital per site by the Hamil¬ 
tonian 

HtJ - 

P [ ~t 0 (Cj^ Cj^ -h h.c.) + 

J I<ij>(SiSj - njOj/A)] P. 

where P is an operator that projects out ail doubly 
occupied sites, Cii and Ci« are creation and annihi¬ 
lation operators for an electron in a Wannier state 
at site i with z-component of spin o. ni is the elec¬ 
tron number operator that counts the total number 
of electrons at site i, and s, is the spin operator. 


This Hamiltonian can be interpreted in two ways: 
(1) In the case of small J, the system can be 
thought of as the large U limit of the single-band 
Hubbard model of electronic systems. The 
antiferromagnetic exchange comes as a result of a 
virtual process where one electron hops onto a 
singly occupied nearest-neighbor site and then 
hops back. The energy gain for such a process is 
of the order of t^/U since a doubly occupied site 
has energy U. (2) The system can also be thought 
of as an electronic system with a ‘super” exclusion 
principle where no two electrons (like or unlike 
spins) are allowed on the same site. The second 
interpretation does not put any restriction on the 
size of J. 

There are few rigorous results available on this 
system: (1) At half filling, the system reduces to a 
Heisenberg antiferromagnet. (2) At 2t = ± J, the 
model has been solved by the Bethe-ansatz tech¬ 
nique. (3) In one dimension and J = 0, the model 
has been solved also by the Bethe-ansatz tech¬ 
nique. Other attempts to study this model have 
focused on the ground state and first few excited 
states. We have decided, instead, to focus on the 
thermodynamic properties of this system. Our 
method, the renormalization-group approach, 
involves solving a statistical mechanics problem by 
a recursive elimination of the degrees of freedom. 
Since the Hamiltonian involves a regular lattice, 
this problem is well suited for the position-space 
renormalization-group method. The solution is 
obtained in an expanded space with Hamiltonian 

H - 

P[-tS<ij>,,(4cja + h.c.)-h 

I<y> (Jsjsj -I- Vnjnj) + /iljnj] P. 

Because of the non-commutativity of quantum 
operators, new techniques of renormaliza¬ 
tion-group theory had to be developed. 

We have obtained the renormalization-group flows 
for the effective coupling constants for one-, two-, 
and three-dimensional systems. The flows deter¬ 
mine the phase diagrams and all thermodynamic 
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properties. tn one dimension, we find, as 
expected, no finite-temperature phase transition. 
In two dimensions, we find a single finite- 
temperature critical point, as previous workers 
conjectured but were unable to derive. In three 
dimensions, our preliminary results indicate a rich 
finite-temperature phase diagram with novel 
phases, phases transition behaviors, and conducti¬ 
vity phenomena. Our results are also confirmed by 
our exact small-cluster calculations. 

1.3 Phase Diagrams of 
Semiconductor Alloys 
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William C. Hoston, Jr., Roland R. Netz, Professor 
A. Nihat Berker 

A study of ternary compounds on face-centered- 
cubic lattices has been started. The aim of the 
work is to elucidate the phase behavior of ternary 
and quaternary semiconductor alloys. These alloys 
have received recent attention for both funda¬ 
mental and technological reasons. It has been 
found that these alloys can exist in the zincblende, 
chalcopyrite, or possibly stannite structures, invol¬ 
ving two interpenetrating fee lattices on which up 
to four atomic species exist. One atomic species 
occupies one of the fee lattices while three other 
atomic species may (chalcopyrite, stannite) or may 
not (zincblende) order on the other fee lattice. At 
present, not much is known about the chalcopyrite 
structure seen in the experimentally obtained lll-V 
compounds. There is great interest in learning 
how stable this chalcopyrite phase is against other 
possible structures. 

In recent work by K.E. Newman and collaborators, 
the Blume-Emery-Griffiths model has been 
adopted for the study of the zincblende to 
chalcopyrite or stannite transitions. This model is 
a spin-1 Ising model with Hamiltonian 

H * J Sj s^ -h K £<jj> Sj Sj — AZjSj, Sj = 0, ± 1. 

The three spin values are each associated with a 
different species of atom. A, B, or C, which exist 
on one of the fee lattices. The other fee lattice is 
considered occupied by atomic species D. The 
systems under consideration have the composition 
[(AB)i.,C20D2. 

The model includes interactions between the A, B, 
and C atoms. The parameters J and K in the Ham¬ 
iltonian above are fixed as combinations of these 
interaction energies. They are chosen so as to give 
the chalcopyrite structure at low temperature 
(J < 0) and to control the phase transition 


between the chalcopyrite and the zincblende. A 
controls the relative densities of the species (AB) 
and C. 

As a preliminary study, we have completed the 
global phase diagram study of the Blume-Emery- 
Griffiths model. We found six new phase dia¬ 
grams, including novel multicritical topology and 
two new ordered phases. Thus, we determined 
that the phase diagram of this simple spin system 
includes nine distinct topologies and three ordered 
phases. These results were obtained by a global 
mean-field theory with four independent order 
parameters. 

It is important to note, however, that the choice of 
J < 0 makes the model on the fee lattice frus¬ 
trated, which requires analysis beyond mean-field 
theory. We have developed a new method, which 
we have called ‘hard-spin mean-field theory,* that 
incorporates the hard-spin condition of local 
degrees of freedom and thereby conserves frus¬ 
tration. We have tested this method on frustrated 
triangular and stacked triangular lattices with the 
spin-1/2 Ising model, obtaining excellent results. 
We are now applying the method to the Blume- 
Emery-Griffiths model. Our future studies of semi¬ 
conductor alloys will also include 
renormalization-group theory. 


1.4 Quantum Spin Systems 
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We are interested in calculating the thermo¬ 
dynamic properties of quantum mechanical 
systems at low temperatures, where the distinctive 
quantum phenomena are important. The first 
model we consider is the s - 1 /2 XXZ model, with 
Hamiltonian 

H = [J (8*8* -H srsj") -»- J, sfs*]. 

This model in two dimensions is of interest in 
high-Tc superconductors; it is also of general 
interest to develop our ability to do the statistical 
mechanics of quantum systems and to extract, 
thereby, the quantitative properties of their ele¬ 
mentary excitations. 

We map this d-dimensional quantum system 
quantitatively onto a (d+1 }-dimensional classical 
system with complicated interactions. For the 
ferromagnetic quantum system, we have been able 
to obtain magnon dispersion relations and 
magnon-magnon interactions. Eventually, we 
expect to apply our method to antiferromagnetic 
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quantum systems, where even the ground state 
presents a many-body problem. 
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2.1 Introduction 
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In this research program we use modern x-ray 
scattering techniques to study structures and 
phase transitions in thin films and on surfaces. We 
have two principal experimental facilities, one at 
MIT and the other at the National Synchrotron 
Light Source at Brookhaven National Laboratory. 
At MIT, we have four high-resolution computer- 
controlled x-ray spectrometers using high- 
intensity, rotating anode x-ray generators. The 
angular resolution can be made as fine as 1.8 
seconds of arc, which enables us to probe the 
development of order from distances of the order 
of the x-ray wavelength ~lA up to 30,000A. The 
sample temperature can be varied between 2 K 
and 500 K with a relative accuracy of 2x10~’ K. At 
the National Synchrotron Light Source, in collab¬ 
oration with IBM, we have three fully instrumented 
beam lines. Two of these beam lines allow us to 
make studies with photons varying in energy 
between 3 and 12 keV; the third has a fixed energy 
of 17 keV. These facilities make possible high- 
resolution scattering experiments with a flux of 
more than three orders of magnitude larger than 
that of a rotating anode x-ray generator, opening 
up a new generation of experiments. 

Several years ago, as part of this JSEP program, 
we built an x-ray-compatible, high vacuum single 
crystal apparatus. This enabled us to use 
synchrotron radiation to study the structures and 
transitions occurring at a single surface, and such 
experiments are now becoming routine. As a 
result of our recent research, we determined that a 
new chamber allowing access to a wider range of 
reciprocal space was required. In collaboration 
with Professor Simon G.J. Mochrie, we have 


designed, built and are now commissioning this 
second-generation x-ray surface facility. 

Our basic scientific objective is to understand the 
morphologies and microscopic structures of simple 
semiconductor and metal surfaces at high temper¬ 
atures. Possible phase changes include surface 
roughening, surface reconstruction, melting, 
amorphization and dilution. These phenomena are 
particularly interesting on stepped surfaces, where 
there may be an interplay between step structures, 
faceting, reconstruction, and roughening. 


2.2 Metal Surface Studies 

In 1987 we reported a detailed study of a possible 
roughening transition of the AgpiO) surface. 
Specifically, we discussed the detailed line shape 
of the Ag(110) surface peak for temperatures 
between 10O'C and 525*C. At we 

observed a crossover to a pure power law 
lineshape, indicative of a rough surface, for the 
near-surface diffraction. Above 450°C the profiles 
could be described by the power law singularity 
form (Q - with rj changing rapidly with 

temperature. Accordingly, we identified 450°C as 
the roughening temperature. 

Our more recent experiments and those by others 
on crystals with very small miscuts (^ 0.5°) 
suggest that the story may be more complicated. 
Specifically, the experiments suggest that below 
~450°C to 500°C, the Ag(110) surface may be 
facetted. In order to explore this further, we have 
begun a systematic set of experiments on vicinal 
Ag(110) crystals with miscuts in the (110) direc¬ 
tion of up to 4°. Our initial experiments on a 
crystal miscut by ~ 2.3° indeed show facetted 
behavior below '-500°C. Further experiments 
should elucidate the full nature of the roughening- 
facetting transition on this prototypical noble metal 
surface. 
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2.3 Semiconductor Surface 
Studies 

2.3.1 Stepped Si(111) Surfaces 

The behavior of vicinal surfaces (surfaces cut a 
few degrees away from low index facets) involves 
essential aspects of surface structures including 
equilibrium crystal shapes, roughening, faceting, 
and surface reconstruction. The stability of vicinal 
surfaces is determined by the shape of the equilib¬ 
rium crystal surface. By studying a series of vicinal 
surfaces with various surface normals, one can 
map out the equilibrium crystal shape. A vicinal 
surface, when it is stable, forms a structure con¬ 
sistent with low-index terraces separated by steps 
whose spacing is usually incommensurate with the 
bulk substrate lattice. As in the case of incom¬ 
mensurate domain wall problems, there can not be 
true long range order in the structure of the incom¬ 
mensurate steps. This is due to the overwhelming 
amplitude of the long wavelength thermal fluctu¬ 
ations in two-dimensional systems with contin¬ 
uous symmetry. As a result, an incommensurate 
stepped surface is rough in an exact sense. A 
vicinal surface becomes unstable when the free 
energy of a nearby low index facet is lowered 
drastically by some mechanisms such as recon¬ 
struction. One might expect that the lowering of 
free energy of a low-index facet would break a 
vicinal surface into a low-index facet and rough 
faces with a higher misorientation angle. 

We have carried out a detailed synchrotron x-ray 
scattering study of a vicinal Si(111) surface miso- 
riented by 3.5° toward the (112) direction. At 
temperatures higher than the ‘7 x 7-to-1 x 1" 
transition, we find that the surface is composed of 
steps with a mean separation of ^ 60A which are 
correlated beyond 3000A. Below the transition, 
the surface is split into the large (> 3000A) per¬ 
fectly flat ( 111 ) facets together with stepped 
regions whose misorientation angle increases as 
temperature is decreased. The transition is first 
order, indicating that it is induced by the 7x7 
reconstruction of the (111) facets. In the facetted 
phase, the facet angle a follows the law tan 
a ~ I Tc — TI *''' ^ with = 2.3i8;5. This disagrees 
somewhat with the theoretical predication A = 3/2 
for models in which the interaction between steps 
is entropic in character. 

At temperatures above 1110 K, we observed two 
step superlattice peaks bracketing the ( 10 ) posi¬ 
tion of the (111) surface. Figure 1 shows high 
resolution longitudinal and transverse scans 
through the first-order step peak at 1120 K. The 
simplest model leads us to believe that this peak 
results from the ordering of steps by the entropic 


repulsion. To test this idea quantitatively, we fit 
the high temperature peak profile to a 2D power 
law singularity of the form 

S(q) = (qS-»-Bq^)<-2+''i)/2 

convolved with the instrumental resolution func¬ 
tion. Figure 1 shows the optimal simultaneous fits 
to the longitudinal and transverse scans with 
h = 0.6 ± 0.1 and B = 15 ± 7. It is evident that 
the power law form fits the data very well. The 
large value of 171 is fully consistent with scattering 
from a rough surface whose height variance 
diverges logarithmically with distance. We note 
that as in the stripe domain commensurate- 
incommensurate problem the Pokrovsky-Talapov 
model leads to the prediction 171 = 2/9, indepen¬ 
dent of temperature. This is in clear disagreement 
with our results. More elaborate models of the 
interaction between the steps are needed to 
explain the large value of 171 . 

More recently, we have been studying the 
behavior of vicinal Si(111) at very high temper¬ 
atures. Remarkably, we find that at ^1300 K there 
is a second phase transition in the step structure 
which mimics a surface melting transition. We are 
currently exploring this new-found phenomenon. 
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Figure 1. Longitudinal and transverse scans through the first order Si(111) step surface peak. The dotted curves 
are the resolution. The solid curves are the result of fits to the power-law lineshape with i; 0.6 and B 15. 
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As we have shown previously, P 2 reacts with 
Si(100) with very close to unity probability 
(0.95) up to one monolayer of coverage, con¬ 
trary to popular belief. We have now also shown 
that the sites for fluorine adsorption under these 
conditions are the dangling bonds and that the 
adsorbed fluorine does not lift the reconstruction. 
This was accomplished using a new technique of 
He atom diffraction that we have recently imple¬ 
mented in our molecular beam-ultrahigh vacuum 
surface scattering apparatus. The experiment is 
carried out by directing a well-characterized, 
monoenergetic beam of He atoms at the surface 
and monitoring the angular distribution of the 
diffracted atoms. This arrangement adds a pow¬ 
erful new structural probe to our experimental 
capabilities. While the binding of fluorine to the 
dangling bonds along with the maintenance of the 
surface dimer bond is not surprising, our results 
show the first experimental evidence for it. No 
experimental information about this system had 
been available previously because the common 
structural probes, electron diffraction techniques, 
are not sensitive to halogens. 


Publications 

Ceyer, S.T. . D.J. Gladstone, M. McGonigal, and 
M.T. Schulberg. 'Molecular Beams: Probes of 
the Dynamics of Reactions on Surfaces.* In 
Physical Methods of Chemistry. 2nd ed. Eds. 
B.W. Rossiter, J.F. Hamilton, and R.C. 
Baetzold. New York: Wiley, 1991. Forth¬ 
coming. 
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We have previously shown that the reaction prob¬ 
ability of F 2 with a clean Si(IOO) surface decays 
from near unity to zero as the fluorine coverage 
increases to one monolayer. This lack of reactivity 
with the fluorinated Si surface is the source of the 
misconception that F 2 does not react with Si. The 
lack of reactivity with the fluorinated surface pre¬ 
cludes the build up of a sufficient layer of fluorine 
to produce the volatile etch product. SiF*. 
However, we have also previously shown that if 
the kinetic energy of the incident F 2 molecule is 
increased above a threshold value of 6 kcal/mol 
(0.25 eV), the dissociation probability of F 2 
with a fluorinated Si surface increases linearly with 
the normal component of kinetic energy. This is a 
result of a barrier to dissociation that is overcome 
by translational energy of the incident molecule. 
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The enhancement in the dissociation probability 
allows enough fluorine to be deposited to form the 
etch product, SiF 4 . This result establishes that Si 
can be etched with low energies using molecular 
beam techniques and without the use of plasmas. 
The low energies afforded by molecular beam 
techniques prevent the introduction of radiation 
damage or defects into the Si lattice, which are 
typical results of plasma etching. 

We have begun a new project aimed at deter¬ 
mining the origin of the barrier to dissociation of 
F 2 on the fluorinated surface. The hypothesis is 
that the barrier to F 2 dissociation is largely associ¬ 
ated with breaking the Si surface dimer bond. To 
test this hypothesis, we are using He atom diffrac¬ 
tion as a probe of the surface structure after expo¬ 
sure to the energetic F 2 . Preliminary measurements 
indicate that the disappearance of the dimer bond, 
as measured by the half order diffraction feature, 
correlates well with the increase in the dissociation 
probability as a function of the F 2 incident 
translational energy in the normal direction. This 
result may be the first determination of the phys¬ 
ical origin of a reaction barrier in any semicon¬ 
ductor system. 

3.3 New Mechanisms for 
Surface Processes 

Sponsor 

National Science Foundation 
Grant CHE 90-20623 
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We report the first detection and identification of 
vibrational modes of buried species by high resol¬ 
ution electron energy loss spectroscopy (HREELS) 
and the observation of a new mechanism for 
absorption of adsorbates, collision-induced 
absorption, as applied to the interaction of 
hydrogen with Ni(i 11). 

Subsurface or bulk sites of a Ni single crystal are 
populated with atomic hydrogen by exposure of 
the (111) face of the Ni crystal at 130 K to atomic 
hydrogen. Atomic hydrogen is generated by 
thermal dissociation of H 2 over a hot tungsten fila¬ 
ment positioned 0.25 inches from the crystal 
surface. This results in a flux of both atomic and 
molecular hydrogen to the front surface of the 
crystal. The H* ions and electrons generated by the 


filament are shown to play no role in the following 
observations. Both atomic H and H 2 result in pop¬ 
ulation of the threefold hollow surface sites up to 
a surface saturation coverage of 1 ML. The 
HR EEL spectrum at 1 ML is characterized by the 
Ni-H asymmetric stretch at 955 cm~' and the Ni-H 
symmetric stretch at 1170 cm-\ and the thermal 
desorption spectrum exhibits two maxima at 340 
and 380 K. Further exposure to atomic H results 
in the appearance of a new feature in the HREEL 
spectrum at 800 cm-' and in the thermal 
desorption spectrum at 180-220 K. The loss 
feature at 800 cm-', for which an impact scat¬ 
tering mechanism is operable, is identified as a 
Ni-H vibration of bulk H by the similarity of the 
dependence of its intensity on electron impact 
energy to the dependence of the electron inelastic 
mean free path on electron energy. At 3 eV impact 
energy, the intensity of the 800 cm-' loss feature 
is large but decreases monotonically to almost zero 
as the impact energy is increased beyond 15 eV. 
This kind of dependence of the intensity of a loss 
feature associated with a bulk vibrational mode is 
expected because the inelastic mean free path of 
an electron, which is 80 A at 3 eV, rapidly drops to 
1 A as the energy is increased above the plasmon 
frequency of the metal (15 eV), where plasmon 
creation is the predominant energy loss mechanism 
which shortens the mean free path. In contrast, 
the intensity of the two loss features associated 
with surface chemisorbed H exhibits almost no 
dependence on electron energy with the exception 
of a resonance at 7 eV. As much as an equivalent 
of 8 monolayers of hydrogen have been absorbed 
into the bulk. 

Subsurface or bulk sites of Ni can also be popu¬ 
lated via a newly observed mechanism, collision- 
induced absorption. A beam of monoenergetic Kr 
or Xe atoms, produced by the supersonic expan¬ 
sion of Kr or Xe seeded in He, is directed at a 
monolayer of atomic hydrogen chemisorbed on 
Ni(111). The impacts of the incident inert gas 
atoms are observed to induce the absorption of the 
chemisorbed hydrogen. The hydrogen thus 
absorbed is identified by the appearance of (1) a 
loss feature at 800 cm-' and (2) a feature at 
180-220 K in the thermal desorption spectrum, 
consistent with the study of atomic hydrogen 
absorption. The absorption probability scales line¬ 
arly with the energy of the incident Xe atom above 
a threshold energy of 2.5 eV. Collision-induced 
absorption may play a significant role in the mech¬ 
anism for hydrogen absorption in the high pressure 
environments of a hydrogen storage cell where the 
hydrogen saturated surface is continually subject 
to bombardment by a large flux of energetic parti¬ 
cles. 
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4.1 Introduction 
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Understanding the properties of surfaces of solids 
and the interactions of atoms and molecules with 
surfaces has been of extreme importance both 
from technological and academic points of view. 
The recent advent of ultrahigh vacuum technology 
has made microscopic studies of well- 
characterized surface systems possible. The way 
atoms move to reduce the energy of the surface, 
the number of layers of atoms involved in this 
reduction, the electronic and vibrational states that 
result from this movement, and the final symmetry 
of the surface layer are all of utmost importance in 
arriving at a fundamental and microscopic under¬ 
standing of the nature of clean surfaces, 
chemisorption processes, and the initial stages of 
interface formation. 

The theoretical problems associated with these 
systems are quite complex. However, we are cur¬ 
rently at the forefront of solving for the properties 
of real surface systems. In particular, we are con¬ 
tinuing toward our goal of calculating the total 
ground-state energy of a surface system from 'first 
principles,' so that we can provide accurate the¬ 
oretical predictions of surface geometries. Our 
efforts in this program have concentrated in the 
areas of surface growth, surface reconstruction 
geometries, structural phase transitions, and 
chemisorption. 


4.2 Heteroepitaxial Growth 

Epitaxial growth of dissimilar semiconductor mate¬ 
rials holds significant potential for technological 
applications and has been the subject of major 
international efforts in recent years. Nevertheless, 
relatively little theoretical work has been performed 
to understand the fundamental interactions gov¬ 
erning the initial stages of growth and the struc¬ 


ture of the first few monolayers in these systems. 
Of particular interest is the prototypical system 
involving growth of GaAs on Si(IOO) substrates. 
Since growth typically begins with an As over¬ 
layer, one first needs to understand this seemingly 
simple system. Experimentally, it is observed that 
the adsorption of As on a stepped Si(100) surface 
can drastically rearrange the distribution of steps of 
the original surface. These changes in surface 
morphology are controlled by different factors, 
including temperature, coverage of As, and the 
type and density of steps on the original surface. 
The structure of the the As-covered Si(100) 
surface [henceforth Si(100):As] determines many 
of the properties of the final epitaxial film. 

In particular, it has been observed that GaAs 
grown epitaxially on Si(IOO) can have two orien¬ 
tations, related by a 90 degree rotation, with 
respect to the Si (100). This puzzling observation, 
the so-called sublattice-orientation (or allocation) 
dilemma, has not been understood. Naively, this 
could be explained if the first epilayer above the Si 
substrate could be selected to be a Ga or As layer. 
However, bond strength and other theoretical 
arguments, as well as experimental evidence, show 
this selective wetting of the Si(IOO) surface by As 
or Ga does not occur. 

In this work we have performed theoretical calcu¬ 
lations which show that As adsorbed on a vicinal 
Si(IOO) surface with double-layer steps either can 
grow directly on top of the Si surface, or can rear¬ 
range the surface so that it replaces the original 
top Si layer. The former structure is metastable 
and results from growth at low temperatures. The 
latter has lower energy and occurs for growth at 
substrate temperatures where surface diffusion is 
activated. Replacing adsorption by substitution 
changes the sublattice of the diamond structure 
exposed at the surface and corresponds to a 90 
degree rotation of the crystal orientation about the 
surface normal. The resulting two configurations 
of the stepped Si(100):A$ surface are related by a 
90 degree rotation, but otherwise differ only in the 
type of steps present at the surface. Their differ¬ 
ence in energy, which we calculate from first prin¬ 
ciples, is related to the relaxation of stress at 
surface steps. This result shows that the rear- 
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rangement of steps upon As adsorption controls 
the orientation of an epitaxial GaAs film. 

Both Si(IOO) and Si(100):As surfaces have 2x1 
reconstructions where surface atoms (Si in the first 
case and As in the second) form dimers arranged 
in parallel rows. A simple picture for the 
adsorption of As on Si(IOO) is that the Si dimers 
break and As dimers are formed on top, perpendic¬ 
ular to the original Si dimers. This rotates the ori¬ 
entation of the surface dimers by 90”. A 
single-domain Si(100) 2x1 surface would result 
in a single-domain Si(100):As 1x2 surface. 
STM images show that the orientation of the As 
dimers depends on the substrate temperature Tt 
during deposition. For low initial temperatures 
(T, ^ 400°C) the surface reconstruction rotates 
from 2x1 to 1x2, as the simple picture 
describes. 

We begin by considering a single-domain Si(IOO) 
2x1 surface, as typically used in growth exper¬ 
iments. These are obtained by using double-layer- 
stepped vicinal surfaces misoriented towards the 
[011 ] direction. Here we will restrict ourselves to 
the case where both the starting Si(IOO) and the 
final Si(100):As surfaces are double-layer stepped 
and the As coverage is one complete monolayer. 
These are the conditions relevant to typical growth 
conditions. On the clean Si(IOO) surface the Si 
dimers are parallel to the step edges. Double-layer 
steps with dimers perpendicular to the edge have 
much higher energy and are not observed. In prin¬ 
ciple, there are two types of double-layer steps on 
the Si(100):As sjrface. Type A, where the As 
dimers are perpendicular to the step edge [Figure 
1(a)], and type B, where the As dimers are parallel 
to the edges [Figure 1 (b)]. The reconstruction of 
the surface is 1 x 2 with type-A steps, and 2x1 
with type-B steps. Figure 1 shows steps with a 
simple edge termination, where Si atoms are 
fourfold coordinated and As atoms are threefold 
coordinated with a doubly occupied lone-pair 
orbital. This is the same bonding configuration 
that passivates the flat Si(100):As surface. Other 
step reconstruction we considered, which are dis¬ 
cussed later in this chapter, were found to have 
higher energies. 

To determine the relative stability of the 1x2 and 
2x1 stepped Si(100):As surfaces, we compare 
the energies of surfaces with type-A and type-B 
steps, respectively. We calculate total energies 
from first principles within the framework of 
density-functional theory in the local-density 
approximation, using norm-conserving pseudo¬ 
potentials. 


(a) SI(100):As1 X2 
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Figure 1. Top and side views of double-layer steps on 
Si(100):As surfaces, (a) Type-A steps: As dimers per¬ 
pendicular to the step edge, (b) Type-B steps: As 
dimers parallel to the step edge (edge relaxation is 
illustrated in the top view, see text). Larger circles on 
top views represent atoms closer to the surface. The 
two_circle sizes on side views represent front and back 
[011 ] planes. 
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The calculated energy difference between type-A 
and type-B steps of Si(100):As is 

* (200 ± 20) meV/a, (1) 

where Xa and Xt are the energies of type*A and 
type-B steps, respectively, and aa>3.80A is the 
unit length along the step. This energy difference 
is large enough to change the surface from one 
step configuration to the other. Indeed, an energy 
difference of 110 meV/a between single- and 
double-layer steps on vicinal Si(IOO) surfaces (for 
a surface misorientation larger than ~ 3-4* ) drives 
the surface to have a majority of double-layer 
steps. It should be noted that the calculated step 
energies represent step-formation plus step-step- 
interaction energies, and the large size of these 
systems makes it prohibitive to increase the sepa¬ 
ration between steps and systematically separate 
these two contributions. However, it is assumed 
that the latter is approximately eliminated from the 
problem by taking energy differences between 
stepped surfaces and that this difference does not 
change significantly with larger step separations. 

An inspection of figure 1 shows no obvious elec¬ 
tronic origin for the large energy difference 
between type-A and type-B steps. In both cases 
the Si and As atoms have ideal bonding configura¬ 
tions, and there is no particular atom or bond 
whose relaxation or rehybridization might explain 
this result. We find that the energy difference is 
related to the relaxation of surface stress. The 
Si(100):As surface is under a large tensile stress, 
both parallel and perpendicular to the As dimers, 
that results from the tendency of the As atoms to 
from 90 degree bond angles and pull themselves 
up from the surface. This stress introduces a 
torque on surface steps, which allows a lateral 
contraction of the surface in the region near the 
steps, at the expense of introducing bulk strain. 

Figure 2 shows the response of type-A and type-B 
steps to the tensile surface stress. Type-A steps 
lead to a large lateral contraction. This contraction 
is ^ 3.5% for the first layer of Si atoms, most of it 
within two lattice constants from the step edge 
(see inset in figure 2), and decays to zero deeper 
into the bulk. The surface As dimers, perpendic¬ 
ular to the step edge, are ~ 2% shorter than their 
calculated value on a flat Si(100):As surface. The 
relaxation type-B steps is different, even though 
the strength the tensile stress is calculated to be 
approximately equal both parallel and perpendic¬ 
ular to the surface dimers [the large energy differ¬ 
ence in equation (1) is consistent with this]. 
Figure 2 shows that type-B steps undergo almost 
no contraction. Still, the relaxed coordinates 
suggest that the surface stress has been signif¬ 
icantly reduced. Although we have not directly 
calculated the stress remaining on the stepped 
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Figure 2. Layer-by-layer lateral contraction for a 
surface with type-A (solid line) and type-B (dashed 
line) steps as a function of penetrate into the bulk. 
Layer No. 1 denotes the first layer of Si atoms below 
the surface As dimers. Inset: Schematic representation 
of lateral contraction for a surface under tensile stress, 
defined with respect to the ideal terrace width shown 
by the vertical dashed lines. 


surface (such calculations are beyond current 
capabilities), the lack of strain in the subsurface of 
type-B steps provides evidence of reduced surface 
stress. On flat Si(100);As subsurface strain forces 
a displacement of atoms from their ideal bulk posi¬ 
tions. Atoms in third and fourth layers that are 
underneath surface dimers shift up, and subsurface 
atoms that are between dimer rows shift down. 
For type-A steps this distortion is further 
enhanced. On the other hand, for type-B steps the 
relaxed positions of the subsurface atoms are close 
to their ideal bulk value. This lack of strain in the 
subsurface for type-B steps is consistent with the 
reduction of surface stress. 

We now address the implications of these results 
to growth experiments. Our calculations show 
that the Si(100):As 1x2 surface with type-A 
steps is a metastable structure. This surface results 
from the growth of As directly on top of the initial 
Si surface. Experimentally, it is obtained by 
depositing As at low substrate temperatures, where 
the As overlayer caps the surface and freezes any 
surface mobility. However, if the temperature is 
high enough during As deposition to activate 
surface diffusion, the surface can then reach the 
lower-energy configuration with type-B steps. 
This explains the surprising experimental observa¬ 
tions. We note that this change in step type 
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requires a large rearrangement of the surface, but 
similar step redistributions requiring mass transport 
across the surface have been observed when the 
initial Si surface has single-layer steps. The trans¬ 
formation from type-A to type-B step may occur 
through other intermediate step configurations, 
including single-layer steps, and a complete con¬ 
version may not occur because of the large 
amount of mass transport required. This may 
further depend on the flux of As, and on whether 
As 2 or As 4 gas sources are used. Also, if surface 
diffusion is completely frozen or defects pin the 
edges of steps, type-A steps with a rebonded edge 
may occur on terraces that originally had an odd 
number of lattice sites, and both types of type-A 
steps may appear on the low-temperature surface. 
In principle, the surface with type-B steps ought to 
be reached from the surface with type-A steps 
upon annealing. However, As begins to leave the 
surface before this change occurs, and mixtures of 
steps with different heights are then observed. 

It should be noted that since the driving force for 
the rotation of the As surface dimers is the differ¬ 
ence in energy between type-A and type-B steps, 
as the surface misorientation becomes smaller and 
the density of steps decreases, the energy gain by 
step rearrangement diminishes. On nearly flat sur¬ 
faces only the 2 X 1 -»1 X 2 change in orientation 
of Si to As surface dimers should be observed. 

For the heteroepitaxial growth of GaAs on 
Si(100), our results predict the sequence of orien¬ 
tations of the Si and then As surface dimers in the 
growth process of 

Si(100)-*Si{100):As-GaAs/Si(100). For low 
initial substrate temperatures, the predicted 
sequence is 2x1-»1x2->4x2. and for high 
initial temperatures the sequence is 
2x1 -»2x1 ->2x4. The final orientation of the 
GaAs epilayer depends on whether the equilibrium 
Si(100):As surface is attained, illustrating the 
interplay between kinetics and energetics in epi¬ 
taxial crystal growth. This result explains the 
sublattice orientation dilemma and brings into a 
coherent picture experiments that seemed to be in 
contradiction. Also, type-A and type-B steps on 
Si(100)'.As may promote different nucleation 
channels for growth, leading to epitaxial films with 
different characteristics. 


4.3 Molecules 

in order to begin to address problems associated 
with molecules interacting with solid surfaces, it is 
important to test the applicability of the theoretical 
tools (developed for condensed matter systems) 
for isolated molecules. Moreover, if successful, 
these tools could be used to study large heteroge¬ 
neous nonperiodic systems such as biological mol¬ 


ecules that may consist of extended organic 
fragments, metal ions, and solvent water mole¬ 
cules. As a first step in this direction, aA initio 
theoretical calculations were performed on a large 
number of molecules with a variety of sizes and 
functional groups consisting of the first-row ele¬ 
ments. These include the simple molecules 
Hi, HiO. CH 4 and NH 3 that are constituents of the 
''biotic-soup' of the primordial atmosphere, neutral 
glycine which is the simplest amino acid, and 
formamide, which models the peptide-bond 
between amino acids in a protein. 

Table I lists the results of our theoretical calcula¬ 
tions (LDAP-PW) and compares them with both 
experimental values and quantum chemistry 
(MP2/6-31G*) Hartree Fock calculations. 

Overall agreement of LDAP-PW with experiment 
and MP2 calculations is excellent. Average devi¬ 
ations from experiment and MP2 are O.OSOA and 
0.034A for bond lengths and 0.5 degrees and 1.0 
degrees for bond angles. If the hydrogen-bonded 
system 3-hydroxy-acrylaldehyde is not included in 
the averages, the bond-length deviations become 
0.01 7A and 0.024A respectively. 

Upon closer examination, certain trends become 
apparent within the data contained in Table I. The 
most striking pattern is that a LDAP-PW calcula¬ 
tion of a bond distance between two heavy (non¬ 
hydrogen) atoms is almost always closer to 
experiment than the corresponding MP2 result. 
Conversely, bonds involving at least one hydrogen 
are almost always calculated more accurately using 
MP2. In addition, bond angles which include zero 
or one hydrogen atom are almost always more 
accurate using LDAP-PW, while those including 
two hydrogens are usually better reproduced by 
MP2. We believe that these effects are due to the 
use of the LDA rather that a more accurate 
exchange-correlation functional. The small size of 
hydrogen makes its charge density change over a 
short length scale. It appears likely that gradient 
corrections to LDA will improve these results. 
Another trend which we have observed is that 
long bonds are reproduced more poorly than short 
bonds by the LDAP-PW. This effect is also pro¬ 
bably due to the use of the LDA. As the atoms in 
a bond separate, they act more like free atoms, and 
the LDA does a poorer job on free atoms than on 
covalent materials. This effect is particularly 
striking in the 0-0 bond of H 2 O 2 and the 
hydrogen bond in 3-hydroxyacrylaldehyde. 

These results show that the LDAP-PW method can 
produce accurate structural parameters for mole¬ 
cules. The data suggest that the method which 
uses plane waves, pseudopotentials, and super¬ 
cells may be a useful alternative to standard at 
initio quantum mechanical methods for calcula- 
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Molecule Parameter 

Exp 

MP2 


LDAP-PW 

4^ 

Aq 

H2 

r(HH) 

0.742 

0.738 

0.004 

0.776 

0.034 

0.038 

02 

KOO) 

1.208 

1.242 

0.034 

1.183 

0.025 

0.059 

N2 

r(NN) 

1.098 

1.131 

0.033 

1.099 

0.001 

0.032 

P2 

r(PP) 

1.893 

1.936 

0.043 

1.875 

0.018 

0.061 

CO 

r(CO) 

1.128 

1.151 

0.023 

1.138 

0.010 

0.013 

OOz 

r(CO) 

1.162 

- 

- 

1.170 

0.008 

- 

HzO 

r(OH) 

0.9S8 

0.969 

0.011 

0.976 

0.018 

0.007 

<HOH) 

104.5 

104.0 

0.5 

104.0 

0.5 

0.0 

CH4 

r(CH) 

1.092 

1.090 

0.002 

1.102 

0.010 

0.012 

NH3 

r(NH) 

1.012 

1.017 

0.005 

1.028 

0.016 

0.011 

<(HNH) 

106.7 

106.3 

0.4 

105.7 

1.0 

0.6 

PH3 

r(PH) 

1.420 

1.415 

0.005 

1.440 

0.020 

0.025 

<(HPH) 

93.3 

94.6 

1.3 

91.6 

1.7 

3.0 

H2O2 

r<00) 

1.452 

1.467 

0.015 

1.414 

0.038 

0.053 

r(OH) 

0.965 

0.976 

0.011 

0.982 

0.017 

0.006 


<00H) 

100.0 

98.7 

1.3 

100.2 

0.2 

1.5 


(0(HOOH) 

119.1 

121.3 

2.2 

118.7 

0.4 

2.6 

N2H2 

r(NN) 

1.252 

1.267 

0.015 

1.233 

0.019 

0.034 

r(NH) 

1.028 

1.036 

0.008 

1.052 

0.024 

0.016 


<(NNH) 

106.9 

105.4 

1.5 

107.2 

0.3 

1.8 

H2C3O 

r(CO) 

1.208 

1.221 

0.013 

1.205 

0.003 

0.016 

r(CH) 

1.116 

1.104 

0.012 

1.135 

0.019 

0.031 


<(HCH) 

116.5 

115.6 

0.9 

116.4 

0.1 

0.8 

C2H2 

r{CC) 

1.203 

1.218 

0.015 

1.196 

0.007 

0.022 

r(CH) 

1.061 

1.066 

0.005 

1.058 

0.003 

0.008 

HNC 

r(NC) 

1.169 

1.187 

0.018 

1.178 

0.009 

0.009 


r(HN) 

0.994 

1.002 

0.008 

1.018 

0.024 

0.016 

HNO 

r(NO) 

1.212 

1.237 

0.025 

1.188 

0.024 

0.049 

r(HN) 

1.063 

1.058 

0.005 

1.093 

0.030 

0.035 


<(HNO) 

108.6 

107.3 

1.3 

109.2 

0.6 

1.9 

Formamide 

r(CO) 

1.193 

1.224 

0.031 

1.217 

0.024 

0.007 

(HCONH2) 

r(CN) 

1.376 

1.361 

0.015 

1.365 

0.011 

0.004 


<OCN) 

124.7 

124.8 

O.I 

124.9 

0.2 

0.1 

Glycine 

r(O0») 

- 

1.218 


1.204 


0.014 

(NHzCHzCQlO^H) 

r(C02) 

- 

1.358 


1.345 


0.013 

r(CC) 

- 

1.515 


1.535 


0.020 


r(NC) 

- 

1.451 


1.460 


0.009 


<01C02> 

- 

123.2 


123.7 


0.5 


<02CO 

- 

111.4 


111.9 


0.5 


<OlCC) 

- 

125.4 


124.4 


1.0 


<CCN) 

- 

114.9 


114.8 


0.1 

H2P04** 

r(PO) 

- 

1.510 


1.482 


0.028 

r{P-OH) 

- 

1.680 


1.640 


0.040 


<HO-P-OH) - 

101.1 


101.2 


0.1 


<(0-P-0H) 

• 

106.9 


107.1 


0.2 


<(0-P-0) 

• 

126.2 


125.2 


1.0 

3-hydroxy- 

r(O-H) 

0.969 

0.9‘'4 

0.025 

1.079 

0.110 

0.085 

acrylaldehyde 

(HOCHCHCHO) 

r(O...H) 

1.680 

l.fc 4 

0.014 

1.418 

0.262 

0,276 


Table I. A comparison of LDAP-PW geometry optimizations with experiment and with MP2. Al and A? are the 
deviations of the MP2 and LDAP-PW calculations from experiment, and Aq is the deviation of LDAP-PW from the 

V3tu6S- 
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tions on chemical systems. Future work will 
include more realistic density functionals and tests 
of energetics of chemical transformations. In addi¬ 
tion, research is continuing to develop more effi¬ 
cient energy minimization algorithms and superior 
pseudopotentials. These changes will improve the 
efficiency and accuracy of the method for the 
study of periodic and aperiodic chemical and 
biochemical systems. 
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5.1 Structure and Phase 
Behavior of the Si(113) Surface 

Sponsor 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
Contract DAAL03-92-C-0001 

Recently, our efforts under JSEP have been 
directed towards implementation of a design for an 
ultra-high-vacuum (UHV) chamber for surface 
x-ray scattering. The design and assembly of the 
UHV chamber and its associated diffractometer are 
now complete, and a first set of experiments using 
the new chamber at the MIT-IBM Beamline (X20) 
at the National Synchrotron Light Source has been 
carried out. 

Our initial program (in collaboration with RJ. 
Birgeneau) is to study the Si(113) surface.' 
Si(113) is interesting for a number of reasons, 
which involve both possible technological applica¬ 
tions and interest in the basic physics of surfaces. 
Firstly, we note that Si(113) is a stable facet with 
a surface energy comparable to that of the (111) 
or (001) orientation.^ Reports of superior epitaxial 
growth on Si(113)3 make the structure of this 
surface and its temperature dependence especially 
worthy of study. 


An unreconstructed (113) surface may be envis¬ 
aged as a sequence of (001) terraces separated by 
(111) steps. (The terraces are a single atomic row 
wide.) The corresponding surface unit cell con¬ 
tains two different types of surface atom: two-fold 
coordinated atoms deriving from the (001) terrace 
with two dangling bonds and three-fold coordi¬ 
nated atoms from the (111) step with a single 
dangling bond (see figure 1, which is taken from 
the article by Myles and Jacobi*). As noted by 
others before,^ this aspect of the (113) surface 
leads to the possibility of two different adsorption 
sites. This feature may alleviate the problem of 
twinning that occurs in the growth of lll-V semi¬ 
conductors on flat Si(OOI).^ Furthermore, in a lll-V 
material, the (113) surface is non-polar, in con¬ 
trast to the (001) and (111) surfaces, which 
makes it especially interesting for heteroepitaxy. 

The clean Si(113) surface reconstructs to form 
either a (1x3) or a (2x3) surface mesh. The (1x3) 
structure appears to be the equilibrium configura¬ 
tion,^ but the precise atomic geometry of the 
Si(113) (1x3) reconstruction has yet to be estab¬ 
lished experimentally Moreover, convincing cal¬ 
culations of the structure of Si(113) do not 
currently exist. At elevated temperatures (;^840K), 
the surface reverts from its (1x3) structure to a 
structure with true (1x1) symmetry of bulk (113) 


' B.F. Olshanetsky and V.l. Mashanov, Surf. Sci. Ill: 414 (1981); D.J. Chadi, Phys. Rev. B 29: 785 (1984). 

* J.M. Gibson, M.L. McDonald, and F.C. Unterwald, Phys. Rev. Lett. 55:1765 (1985). 

3 G.H. Olsen, T.J. Zamerowski, and F.Z. Hawrylo, J. Cryst. Growth 59: 654 (1982); S.L. Wright, M. Inada, and H. 
Kroemer, J. Vac. Sci. Techno). 21: 534 (1982). 

* U. Myles and K. Jacobi, Surf. Sci. 220: 353 (1989). 

B O.J. Chadi, Phys. Rev. B 29: 785 (1984); U. Myles and K. Jacobi, Surf. Sci. 220: 353 (1989); J. Knall, J.8. 
Pethica, J.D. Todd, and J.H. Wilson, Phys. Rev. Lett. 66:1733 (1991). 

* K. Fujikawa, J. Cryst. Growth 77: 509 (1986). 

t Y.N. Yang, E.D. Williams, R.L. Park, N.C. Barlett, and T.L. Einstein, Phys. Rev. Lett. 64: 2910 (1991). 
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Figure 1. 


pianesJ It has been suggested that the phase scales of & 4000A. (Other measurements establish 

transformation between these two surface config- that the existence of peaks associated with a (2x3) 

urations falls into a particular “universality class," structure is excluded at the 10"^ level.) The addi- 

that of the chiral, three state Potts model.^ Our tional profiles in figure 2 illustrate the behavior 

initial experiments, aimed at studying this aspect, observed on warming through the phase transfer- 

were carried out on a (1x3) reconstructed surface, mation. The commensurate 1/3-order rod 
produced by annealing at 1500 K and subse- becomes weaker, and at the same time a broad, 

quently cooling to below the (1x3)-to-(1x1) incommensurate peak may be seen emerging 

transformation temperature (^840K). The presence towards smaller wave vector (smaller K). Least- 

of commensurate ‘1/3-order' rods, located mean-squares fits (solid lines) reveal that the 

between the rods of scattering associated with width of the incommensurate peak is proportional 

bulk periodicities, indicates a (1x3) reconstruction to its displacement from the commensurate posi- 

of the crystal surface. A scan across one such tion, suggesting that the disordering indeed falls 

rod—the (4,2/3) rod (we employ the units of the within the universality class proposed, 

unreconstructed surface unit cell)—is shown for . 

several temperatures in figure 2 on a logarthmic Experimerits to investigate the detailed atomic 

intensity scale. The lowest profile shown corre- structure in the commensurate phase are planned 

spends to a temperature of 750 K. Evidently, the fo’’ near future. Another interesting possibility 
reconstruction is readily observable. In addition, that of surface roughening, which may occur at 
the width of this profile determines that the recon- temperatures. Experiments to investigate 

struction is translationaily well-ordered on length 1^'s aspect are also forthcoming. 


« Y.N. Yang, E D. Williams, R.L. Park, N.C. Barlett, and T.L. Einstein, Phys. Rev. Lett. 64: 2910 (1991); D.A. Huse 
and M E. Fisher, Phys. Rev. B 29: 239 (1984). 
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Figure 2. 
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2.2 Millimeter-Wave Frequency 
Measurement of the Rydberg 
Constant 
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Contract DAAL03-89-C-0001 
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The Rydberg constant, Roo. relates the wavelengths 
of the spectrum of atomic hydrogen to practical 
laboratory units. As such, R^o is the natural unit 
for measurements of numerous atomic properties. 
Experiments using laser spectroscopy have deter¬ 
mined Roo to nearly 1 part in 10’°.* Although R*. is 
the most accurately measured fundamental con¬ 
stant, high-precision experiments which depend 
on Roo as an auxiliary constant demand a more 
accurate determination. For example, new exper¬ 
iments to determine the Lamb shift in the ground 
state of atomic hydrogen will be limited by uncer¬ 
tainties in Roo- The accuracy of optical measure¬ 
ments has reached the limit of precision of 
wavelength metrology. Further progress in meas¬ 
uring Roo requires a frequency measurement, 
making use of the modern definition of length in 
terms of time intervals and the defined speed of 
light. 

We are attempting to advance the accuracy of Roo 
by directly measuring cRoo , the “Rydberg fre¬ 
quency.' By initially preparing highly excited 
“Rydberg" states of atomic hydrogen, we are able 
to measure millimeter-wave transitions to nearby 
states with the full precision of frequency 
metrology. 

The goal of our experiment is three-fold. First is 
the evaluation of cRoc • Second is the measure¬ 
ment of the ground state Lamb shift. Because our 
measurements involve high angular momentum 
states for which the Lamb shift is extremely small, 
a comparison of our results with optical measure¬ 
ments of transitions between low-lying states will 
yield an improved measurement of the Lamb shift. 


Third is the precise frequency calibration of the 
spectrum of hydrogen to provide an independent 
check on forthcoming optical frequency measure¬ 
ments based on laser spectroscopy. 

Our experiment is performed in an atomic beam 
configuration to provide a long interaction time 
and reduce Doppler and collisional perturbations. 
Atomic hydrogen is excited to the state (n-29, 
m=0) by two-photon stepwise absorption. The 
atoms are then transferred to the longer lived “cir¬ 
cular" state (n=29,m=28) by the method of 
crossed electric and magnetic fields.^ The atoms 
then enter a region of uniform fields in which the 
frequency of the transition {n=29.m=28) -» 
(n=30,m=29) is measured. The atoms interact 
with the millimeter-wave radiation at two locations 
in a Ramsey separated oscillatory fields geometry. 
The final state distribution of the atoms is meas¬ 
ured by a state-selective electric field ionization 
detector. The resonance signal is observed as a 
transfer of atoms from the n=29 state to the n=30 
state as the millimeter wave frequency is tuned 
across the transition. 

Figure 4 illustrates the main features of the appa¬ 
ratus. Atomic hydrogen or deuterium is produced 
by dissociation of H 2 or D 2 in a radio frequency 
discharge. The beam is cooled by collisions with 
the walls of a cryogenically cooled thermalizing 
channel in order to slow the beam and thereby 
increase the interaction time. The atoms are 
excited to the state (n=29,m=28) by two-photon 
stepwise excitation in the circular state production 
region. The development of the hydrogen beam 
and optical systems was described in the 1990 
RLE Progress Reporl Number 133. The magnetic 
field necessary to transfer the atoms to the circular 
state is provided by permanent magnets. The elec¬ 
tric field is produced by an arrangement of strip 
electrodes that causes the direction of the field to 
rotate. A detector in the circular state production 
region monitors the efficiency of the laser excita¬ 
tion and momentum transfer processes. 

After the atoms are prepared in the circular state, 
the beam enters the interaction region. Because 
Rydberg atoms interact strongly with external 
fields, accurate measurement of the energy level 
structure requires careful control of the interaction 
environment. Thermal radiation is reduced by 
cooling the interaction region to — 10 K by a 
liquid helium flow system. The ambient magnetic 
field is shielded by a double-wall high permeability 


* M.G. Boshier et al., Phys. Rev. A 40: 6169 (1989); P. Zhao et al.. Phys. Rev. A 39: 2888 (1989); F. Biraben et al., 
Rhys. Rev. Lett. 62: 621 (1989). 

■ 0. Delande and J.C. Gay, Europhys. Lett. 5: 303 (1988). 
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shield. The small electric field, necessary to define 
the quantization axis of the atoms, is applied with 
high uniformity by field plates with corrective strip 
electrodes along the sides. The millimeter waves 
intersect the atomic beam at two locations sepa¬ 
rated by 50 cm. The millimeter-wave optics were 
described in the 1990 RLE Progress Report 
Number 133. 

The state distribution of the atoms emerging from 
the interaction region is analyzed by an electric 
field ionization detector. The atoms enter a region 
of increasing electric field produced by a ramped 
plate held at constant potential. The atoms in the 
state (n=30), which ionize at lower field, are col¬ 
lected in the first detector, while atoms in the state 
(n=29), which ionize at higher field, are collected 
at the second detector. An important feature of 
the detector is that it provides time resolution. The 
atomic beam is pulsed, and the time-resolved 


signal from the detector allows different velocity 
classes to be observed independently. 

During this past year we have observed millimeter 
wave transitions in each of the two millimeter 
wave beams. In this fashion, we have observed 
the (n=29 -♦ n=30) transition in both hydrogen 
and deuterium, as well as the (n’=27 -» n=28) 
transition in hydrogen. A typical resonance signal 
is shown in figure 5. In addition, we have 
observed 'Rabi oscillations' in each of the milli¬ 
meter wave beams. Recently, we have begun 
experiments to study signals from interactions with 
both millimeter-wave beams in the Ramsey sepa¬ 
rated oscillatory fields geometry. Sample interfer¬ 
ence data is shown in figure 6. 

Improvements to increase the signal-to-noise ratio 
are underway, and we anticipate starting on high 
accuracy measurements in the coming year. 


ATOMIC BEAM APPARATUS 


Millimeter—wave interoction regions 
(Ramsey spectroscopy) 



Cold Hydrogen 
beam source 


Laser Excitation and > 

circular state production 

State selective detection 


50 cm 
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Figure 5. Typical signal for the (n=29,m=28) -» (n°30,m=29) transition in atomic deuterium. The width of the 
resonance line is close to the time-of-flight limit as the atoms cross the millimeter wave beam. 




Figure 6. First results for the a narrow tinewidth 'Ramsey' fringe for the transition n«29.m'*28 -» n>30,m«29) in 
deuterium. 
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2.3 Atom Interferometry 
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Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
Contract DAAL03-92-C-0001 
U.S. Army Research Office 
Contract DAAL03-89-K-0082 
U.S. Navy - Office of Naval Research 
Grant N00014-89-J-1207 
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John E. Berberian, Christopher R. Ekstrom, Troy D. 
Hammond, David W. Keith, Marc-Oliver Mewes, 
Dr. H. Joerg Schiedmeyer, Quentin A. Turchette, 
Professor David E. Pritchard 

During 1991, we demonstrated the first atom inter- 
ferometer.^ Using transmission gratings, which we 
fabricated,^ as optical elements for atom deBroglie 
waves, we constructed a three grating atom inter¬ 
ferometer which physically separates atom waves 
before recombining them. Our demonstration was 
closely followed by two demonstrations of atom 
interferometers which used laser light as the beam 
splitters.* 

Atom interferometers will make possible quali¬ 
tatively new types of experiments involving inertial 
effects, studies of atomic and molecular properties, 
and tests of fundamental physics, and may ulti¬ 
mately open the way to making ultra-small struc¬ 
tures using atom holograms. 

• The relatively large mass and low velocity of 
atoms makes atom interferometers especially 
sensitive to inertial effects such as rotation, 
acceleration, and gravity. Sagnac rotation has 
been observed in agreement with theoretical 
predictions,* and sensitivity to gravitational 
acceleration at the 3 x 10~* level has been 
demonstrated.’® Atom interferometers may 
become the best absolute accelerometers and 
gravimeters in the next few years. 


• Atom interferometers can be applied to a 
number of experiments in fundamental 
physics; tests of quantum mechanics such as 
the Aharonov-Casher effect,'’ measurement of 
the equality of proton and electron charges, 
and a precise measurement of the momentum 
of a photon. This latter measurement should 
produce a new high precision value for the 
fundamental constants Na^. 

• interferometers for atoms and molecules will 
offer more accurate ways to measure intrinsic 
properties of these particles, like their polariza¬ 
bility. They will also open up new areas of 
study, such as measurements of the 'index of 
refraction' of a gas for a particle beam which 
passes through it. 

Our interferometer consists of three 400 nm period 
transmission gratings, mounted 0.66 m apart on 
separate translation stages inside the vacuum 
envelope. During operation, the O'" and 1” order 
beams from the first grating strike the middle 
grating (which is 140 pm wide) where they are 
diffracted in the 1** and —1” orders so that they 
converge at the third grating. At the second 
(middle) grating the beams have widths of 30 pm 
(FWHM) and are separated by 27 pm. The first 
two gratings form an interference pattern in the 
plane of the third grating, which acts as a mask to 
sample this pattern. The detector, located 0.30 m 
beyond the third grating, records the flux trans¬ 
mitted by the third grating. Figure 7 shows the 
design of the interferometer. 

The data necessary to determine the interferometer 
phase contrast are acquired by modulating the 
position of one grating relative to the other two 
and simultaneously recording the signal from the 
atom counting electronics as well as the signal 
from an optical interferometer used to measure the 
relative position of the gratings. After removing 
data obscured by noise spikes from the hot wire, 
the atom count rate data are averaged into bins 
according to relative grating position, resulting in 
the fringe pattern shown in figure 8. 


* D.W. Keith, C.R. Ekstrom, Q.A. Turchette, and D.E. Pritchcard, Phys. Rev. Lett. 66: 2693 (1991). 

’ D.W. Keith and M.J. Rooks, J. Vac. Sc/. Zechnol. B 9: 2846 (1991). 

* F. Riehle, T. Kisters, A. Witte, J. Helmcke, and J. Borde, Phys. Rev. Lett. 67: 177 (1991); M. Kasevich and S. 
Chu, Phys. Rev. Lett. 67: 181 (1991). 

* F. Riehle, T. Kisters, A. Witte, J. Helmcke, and J. Borde, Phys. Rev. Lett. 67:177 (1991). 

’0 M. Kasevich, S. Chu, Phys. Rev. Lett. 67:181 (1991). 

” Y. Aharonov and A. Casher, Phys. Rev. Lett 53: 319 (1984). 
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Figure 7. Our current atom interferometer with laser interferometer stabilization system (Not to scale). 


The peak-to-peak amplitude of our interference 
signal is 70 Hz, which enables us to determine the 
interferometer phase to a precision of 0.1 rad in 1 
min. The excellent long-term stability of our posi¬ 
tion stabilization system provides measured atom- 
interferometer phase drift of less than 0.1 rad over 
10 minutes. 

The key component of our interferometer is the set 
of three matched transmission diffraction gratings 
which we constructed during two visits at the 
National Nanofabrication Facility (NNF) at Cornell 
University. The first trip resulted in the develop¬ 
ment of a new process for fabricating atom optics. 
The process allows fabrication of precisely posi¬ 
tioned openings in thin silicon nitride membranes 
mounted in silicon frames. The pattern created in 
the membrane is determined by an electron beam 
writer, making the process quite versatile. This 
process was used to create the diffraction gratings 
used in the interferometer. In addition, several 
zone plates (atom lenses) were also built and suc¬ 
cessfully demonstrated later. During our second 
trip, we devised ways to reduce the electron beam 
writing time. This decreased thermal drift during 
the writing period assured higher overall accuracy 
and also increased our overall productivity. We 
made a wide variety of diffraction gratings with 
various heights and periods between 100 and 300 
nm as well as an assortment of single and double 
slits. 

We have recently started construction of an inter¬ 
action region that will allow us to perform several 
experiments including measurement of the 
Aharonov-Casher phase shift, the electric polariza¬ 
bility of the ground state of sodium, and the index 


of refraction for sodium passing through noble gas 
atoms. 

The scientific future of atoms interferometers looks 
bright: (1) Atom beam sources are inexpensive 
and intense relative to other particle 
beams/sources (e.g., neutrons, electrons), (2) 
several techniques have now been demonstrated 
to make interferometers for them, and (3) the 
atoms which may be used in them come with a 
wide range of parameters such as polarizability, 
mass, and magnetic moment. This assures the 
applicability of these instruments to a wide range 
of measurements of both fundamental and prac¬ 
tical interest. 



PaitaiiaaB 


Figure 8. Interference signal from 400 seconds of 
data. Hot wire detector background of 40 Hz is sub¬ 
tracted. The solid line is a least squares fit to the data 
with a 400 nm period. The error bars are purely statis¬ 
tical and slightly underestimate the super-Poissonian 
character of the hot wire background. 
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2.4 Cooling and Trapping 
Neutral Atoms 

Sponsors 
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Our current objective is to develop an intense 
source of slow atoms which will be used in 
studies of cold collisions, atom optics and atom 
interferometry. Slow atoms sources typically have 
fluxes of a few times 10^ atoms/sec. Improve¬ 
ments by two to three orders of magnitude seem 
feasible and are crucial for many experiments. 

Experiments with dense samples of cold neutral 
atoms promise new exciting discoveries in basic 
and applied physics. Due to the considerably 
reduced thermal motion of atoms they are ideal for 
high resolution spectroscopy and for more accu¬ 
rate atomic frequency standards. 

• Collisions of ultra-cold atoms in such samples 
are characterized by a long deBroglie wave¬ 
length and dominated by weak long-range 
interactions. Since the collision duration for 
slow atoms greatly exceeds the radiative 
decay time, stimulated and spontaneous 
radiative transitions can take place during the 
collision. Slow collisions are therefore radi¬ 
cally different from fast collisions studied so 
far and will become an exciting new field of 
atomic physics.’* 

• High density samples of atoms open possibil¬ 
ities for observing quantum collective effects 
such as Bose-Einstein condensation and col¬ 
lectively enhanced or suppressed radiative 
decay. 

In 1991, we designed and built two sources of 
cold atoms based on two different principles; a 
“Diffuse Light Slower' and an 'Inverted Zeeman 
Slower.' We also completed the data analysis and 
modeling of the laser spectroscopy and laser 
cooling experiments performed in a magnetic 
trap.” In addition, one of us (K. Helmerson) com¬ 
pleted his doctoral dissertation.’* 

The following discusses in more detail the pro¬ 
gress made in the three areas mentioned above. 


” P.S. Julienne and J. Vigu, 'Cold Collisions of Ground- and Excited-state Alkali-metal Atoms,' Phys. Rev. A 44: 
4464 (1991). 

T3 K. Helmerson, A. Martin, and D.E. Pritchard, 'Laser and rf Spectroscopy of Magnetically Trapped Neutral Atoms,' 
submitted to J. Opt Soc. Am. B; K. Helmerson, A. Martin and D.E. Pritchard, 'Laser Cooling of Magnetically 
Trapped Neutral Atoms," submitted to J. Opt Soc. Am. B. 

K. Helmerson, Laser Cooling and Spectroscopy of Magnetically Trapped Atoms, Ph.D. diss.. Dept, of Physics, MIT, 
1991. 
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2.4.1 Slowing Atoms with Diffuse 
Light 

We have demonstrated a new technique for 
slowing and brightening atomic beams which uses 
isotropic monochromatic laser light detuned to the 
red of the atomic resonance. The atoms compen¬ 
sate for the changing Doppler shift as they decel¬ 
erate by absorbing photons at a variable angle in 
accordance with the Doppler resonance condition. 
The use of isotropic light and the automatic angle 
selection method of compensating for the 
changing Doppler shift distinguish diffuse light 
slowing from other light slowing schemes realized 
so far, and is technically simpler to implement. 

In our experiments isotropic light was generated 
by shining laser light into a tube of diffusely 
reflecting material around the atomic beam (figure 
9). Sodium atoms were slowed from approximately 
300 m/s to below 100 m/s resulting in a contin¬ 
uous slow beam with a flux of greater than 10* 
atoms/s.''* This performance is comparable with 
other slowing schemes and is limited by the light 
losses in the diffuse reflector. We are currently 
working on ways to reduce the loss in the diffuse 
reflector well below its current 1 per cent. 

Since isotropic light cools all components of the 
velocity, a diffuse light slower produces a slow 
beam with a smaller divergence. Other advantages 
of diffuse light slowing are the simple experimental 
realization, tolerance to laser jitter, and absence of 
an intense slowing laser beam collinear with the 
slowed atoms which can interfere with subsequent 
experiments. 


2.4.2 Inverted Zeeman Slower 

In a Zeeman slower the changing Doppler shift as 
an atom slows is compensated by the Zeeman shift 
in an inhomogeneous magnetic field. This method 
has the advantage of producing a continuous 
beam of slow atoms and featuring practically 
unlimited velocity capture range. The original 
implementation of Zeeman slowing, however, 
encountered difficulties in producing beams of 
atoms with velocities lower than 200 m/s. The 
major problems were off-resonant slowing of 
atoms after they left the slower and substantial 



Figure 9. Schematic diagram of the experimental 
setup. A thermal sodium beam passes through the bore 
of a tube of diffusely reflecting material which is filled 
with laser light via six equally spaced small aperatures. 
Fluorescence-excitation spectra produced by laser 
beams perpendicular and at 45° to the atomic beam are 
used to obsen/e the velocity distribution of the atoms as 
they emerge from the diffuse light cavity. 

spreading of the slow atomic beam due to trans¬ 
verse heating. 

To solve the first problem we chose the 
I F=2,m=2 > -♦ I F=3,m=3 > cycling transition to 
slow sodium atoms. The transition requires a mag¬ 
netic field increasing along the atoms trajectory 
(Inverted Zeeman slower). The magnetic field 
drops off rapidly after reaching its maximum thus 
quickly shifting atoms out of resonance and 
reducing off-resonant slowing. It seems therefore 
possible to perform subsequent atomic optics 
experiments without undue interference from the 
intense slowing laser beam. 

We have built and tested an inverse Zeeman 
slower which allows for optical access to the slow 
atoms inside. Transverse laser cooling will increase 
the brightness of the resulting slow atom beam. 
Experiments are still in progress but preliminary 
results indicate that the performance is better than 
any other slower demonstrated so far. We have 
slowed about 50 per cent of the thermal beam 
down to velocities between 40 m/s and 100 m/s. 
After optimization we intend to use the slower for 
efficient loading of an atom trap. 


W. Ketterle, A.G. Martin, M. Joffe and D.E. Pritchard, 'Slowing Atoms with Diffuse Laser Light,' postdeadline 
paper. Seventh Interdisciplinary Laser Conference (ILS), Monterey, California, September 22-26,1991. 
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2.4.3 Laser Cooling and 
Spectroscopy of Magnetically 
Trapped Neutral Atoms 

During 1991, we concluded the data analysis and 
modeling of the experimental results obtained pre¬ 
viously in a superconducting magnetic trap. This 
work can be divided into two parts: 

1. Laser and rf spectroscopy of the trapped 
atoms: 

Laser absorption, fluorescence and rf spectra 
of the trapped atoms were modeled and inter¬ 
preted to extract the energy distribution of the 
trapped atoms. The results of our model show 
that the atoms, when loaded into the trap, had 
an effective temperature of 50 mK, compa¬ 
rable to the actual trap depth. We have 
shown how the rf spectra can give the energy 
distribution of the trapped atoms without prior 
knowledge of the shape of the trapping field 
potential. For our experiments, the energy 
distribution of the trapped atoms was well 
described by a Boltzman distribution, trun¬ 
cated at the energy corresponding to the 
highest energy trapped atom.'* 

2. Laser Doppler cooling of the trapped atoms: 

Our experiments have shown also how simple 
one-dimensional optical molasses can reduce 
the temperature of the trapped atoms to ~ 1 
mK. The motion transverse to the cooling 
laser is cooled by the coupling between the 
different translational degrees of freedom pro¬ 
vided by the trap. The cooling process was 
modeled, and the theoretical results show 
good qualitative agreement with the exper¬ 
imental data. The measurements indicated, 
and the model confirms, that for cooling rates 
that are low compared to the trap coupling 
time, ultimate temperatures a few times the 
Doppler limit can be achieved.'^ 
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2.5 Precision Mass 
Spectroscopy of Ions 
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In 1991, we continued our program to substan¬ 
tially improve our precision mass measurement 
experiment. These improvements should allow us 
to reach a precision of about 10“” in our mass 
measurements of individual atomic and molecular 
ions, the next step toward our ultimate goal of a 
few parts in 10'^. This capability will allow us to 
do a variety of experiments which address issues 
in both fundamental and applied physics, 
including: 


'* K. Helmerson, A. Martin, and D.E. Pritchard, ‘Laser and rf Spectroscopy of Magnetically Trapped Neutral Atoms," 
submitted to J. Opt. Soc. Am. B. 

K. Helmerson, A. Martin and D.E. Pritchard, ‘Laser Cooling of Magnetically Trapped Neutral Atoms,' submitted to 
J. Opt Soc. Am. B. 


174 RLE Progress Report Number 134 




Chapter 2. Basic Atomic Physics 


• The - ^He'^ mass difference, important in 
ongoing experiments to determine the electron 
neutrino rest mass: 

• Determination of excitation and binding ener¬ 
gies of atomic and molecular ions by 
weighing the small decrease in energy. 
Am = Ebind/^^' 

• Determination of Avogadro's number Na, by 
weighing y-rays—its accurate determination 
would permit the replacement of the “artifact" 
mass standard by an atomic mass standard: 
and 

• Improvement of many traditional applications 
of mass spectroscopy by orders of magnitude 
improvement in both accuracy and sensitivity. 

Our experimental approach is to measure ion 
cyclotron resonance on a single molecular or 
atomic ion in a Penning trap, a highly uniform 
magnetic field with axial confinement provided by 
weaker electric fields. We monitor the ion's oscil¬ 
lation along the magnetic field lines by detecting 
the currents induced in the trap electrodes. 
Working with only a single ion is essential because 
space charge from other ions leads to undesired 
frequency shifts. This work in trapping and preci¬ 
sion resonance draws on techniques developed by 
Hans Dehmelt at the University of Washington and 
Norman Ramsey at Harvard, for which they shared 
the 1989 Nobel Prize. 

We have developed techniques for driving, 
cooling, and measuring the frequencies of all three 
normal modes of Penning trap motion. Thus we 
can manipulate the ion sition reproducibly to 
within 30 microns of the center of the trap, cor¬ 
recting for electrostatic shifts in the cyclotron fre¬ 
quency to great accuracy. We use a rr-pulse 
method to coherently swap the phase and action 
of the cyclotron with the axial modes.’® Therefore, 
although we detect only the axial motion directly. 


we can determine cyclotron frequency by meas¬ 
uring the phase accumulated in the cyclotron 
motion in a known time interval (figure 10). 

In the past two years, we have built an entirely 
new Penning trap and detector, including a 
higher-Q resonant circuit and quieter RF SQUID, 
these have improved our signal-to-noise ratio by a 
factor of two. We have also constructed a new 
highly stable DC electric field supply, and added a 
pressure regulator to the liquid helium bath of our 
superconducting magnet to help stabilize the field, 
in addition, we have begun construction of ion 
optics and an external ion source that will allow us 
to make the ions in a discharge at room temper¬ 
ature and then load them into the trap. This will 
eliminate the problem of residual neutral gas in the 
trap when using volatile species such as hydrogen 
and helium, and should permit rapid cycling 
between two different ions which will dramatically 
increase measurement precision. 

With all the improvements to the system, we 
foresee being able to resolve the existing 10 eV 
discrepancy between measurements of 
m(®H) - m(®He) in the next year. In addition, we 
plan to demonstrate a classical squeezing tech¬ 
nique which should reduce the thermal fluctu¬ 
ations of our measurement by a factor of three to 
five. After that, we plan to continue development 
of techniques to measure two ions of different 
mass simultaneously. The two-ion technique in 
combination with the various improvements made 
over the last two years should lead to precision in 
the range of 10"”.’* 

Publications 

Cornell, E.A., K.R. Boyce, D.L.K. Fygenson, and 

D.E. Pritchard. 'Two Ions in a Penning Trap: 

Implications for Precision Mass Spectroscopy." 

Phys. Rev. A. Submitted for publication. 


’® E.A. Cornell, R.M. Weisskoff, et al., Phys. Rev. A 41: 312 (1990). 

’9 E.A. Cornell, Mass Spectroscopy Using Single ton Cyclotron Resonance, Ph.D diss.. Dept, of Physics, MIT, 1990. 
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Time between drive pulse and n-pulse tsecl 


Figure 10. For each plotted point, we perform the following experiment: The initially cold ion is pulsed into a 
cyclotron orbit of known initial phase and then allowed to evolve *in the dark for an indicated amount of time, t. 
Then a pulse is applied which exchanges cyclotron and axial motions, bringing the ion's cyclotron action arid phase 
into the axial mode. As the ion's axial motion rings down, its phase is detected. The appropriate multiple of 360 is 
added, and a line is fitted to the points. The slope of the line is the frequency difference between the frequency 
generator and the trap cyclotron frequency. 
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Active MMW images ere also generated using 
models and software from related projects. Prelim¬ 
inary simulation of 35 GHz sensor scene reveals 
the lack of depth information and resolution, but 
shows good weather penetration. Future work is 
required in the area of postprocessing in order to 
achieve better images. 

1.3 Multilayer Media and 
Superconducting Electronics 

Sponsors 
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Lam, Michael C. Moldoveanu, Hongsing Lee, Ilya 
Lyubomirsky 

A method for the calculation of current distribu¬ 
tion, resistance, and inductance matrices for a 
system of coupled superconducting transmission 
lines having finite rectangular cross-section is pre¬ 
sented. These calculations allow accurate charac¬ 
terization of both high-Te and low-Te 
superconducting strip transmission lines. For a 
single stripline geometry with finite ground planes, 
the current distribution, resistance, inductance, and 
kinetic inductance are calculated as a function of 
the penetration depth for various film thicknesses. 
These calculations are then used to determine the 
penetration depth for Nb, NbN, and YBa 2 Cu 307 -« 
superconducting thin films from the measured 
temperature dependence of the resonant frequency 
of a stripline resonator. The calculations are also 
used to convert measured temperature dependence 
of the quality factor to the intrinsic surface resist¬ 
ance as a function of temperature for an Nb strip¬ 
line resonator. 

The frequency-dependent resistance and induc¬ 
tance of uniform transmission lines are calculated 
with a hybrid technique that combines a cross- 
section coupled circuit method with a surface inte¬ 


gral equation approach. The coupled circuit 
approach is most applicable for low-frequency cal¬ 
culations. while the integral equation approach is 
best for high frequencies. The low-frequency 
method consists in subdividing the cross section 
of each conductor into triangular filaments, each 
with an assumed uniform current distribution. The 
resistance and mutual inductance between the fila¬ 
ments are calculated, and a matrix is inverted to 
give the overall resistance and inductance of the 
conductors. The high-frequency method expresses 
the resistance and inductance of each conductor in 
terms of the current at the surface of that con¬ 
ductor and the derivative of that current normal to 
the surface. A coupled integral equation is then 
derived to relate these quantities through the dif¬ 
fusion equation inside the conductors and 
Laplace's equation outside. The method of 
moments with pulse basis functions is used to 
solve the integral equations. An interpolation 
between the results of these two methods gives 
very good results over the entire frequency range, 
even when few basis functions are used. Results 
for a variety of configurations are obtained and are 
compared with experimental data and other 
numerical techniques. 

A spectral domain dyadic Green's function formu¬ 
lation defining the fields inside a multilayer chiral 
medium due to arbitrary distribution of sources is 
developed. The constitutive parameters and the 
chirality of each layer are assumed to be different. 
The fields are obtained in terms of electric and 
magnetic type dyadic Green's functions. The sin¬ 
gular behavior of these dyadic Green's functions in 
the source region is taken into account by 
extracting the delta function singularities. The 
fields in any layer are obtained in terms of appro¬ 
priately defined global reflection and transmission 
matrices. 

In order to understand the physical meaning of 
rational reflection coefficients in one-dimensional 
inverse scattering theory for optical waveguide 
design, we have studied the relation between the 
poles of the transverse reflection coefficient and 
the modes in inhomogeneous dielectrics. By using 
a stratified medium model, it is shown that these 
poles of the reflection coefficient have a 
one-to-one correspondence to the discrete modes, 
which are the guided and leaky modes. The radi¬ 
ation modes have continuous real values of trans¬ 
verse wave numbers and are not represented by 
the poles of the reflection coefficient. Based on 
these results, applications of the Gel'fand-Levitan- 
Marchenko theory to optical waveguide synthesis 
with the rational function representation of the 
transverse reflection coefficient are discussed. 

We developed an inversion algorithm referred to as 
the Renormalized Source-Type Integral Equation 
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approach. The objective of this method is to over¬ 
come some of the limitations and difficulties of the 
iterative Born technique. It recasts the inversion, 
which is nonlinear in nature, in terms of the sol¬ 
ution of a set of linear equations; however, the 
final inversion equation is still nonlinear. The 
derived inversion equation is an exact equation 
which sums up the iterative Neuman (or Born) 
series in a closed form and; thus, is a valid repre¬ 
sentation even in the case when the Born series 
diverges; hence, the name Renormalized Source- 
Type Integral Equation (STIE) Approach. This 
renormalized STIE approach has been applied to 
underground soil moisture inversion and formation 
inversion in the borehole geometry. 

Several recent papers discuss techniques for 
reducing the parasitic circuit capacitance of 
millimeter-wave planar air bridge type GaAs 
Schottky barrier diodes by completely etching 
away the semi-insulating material under the active 
region and transferring the device to a second 
carrier that is better suited, both electrically and 
mechanically, to the intended high frequency 
application. This technique and a similar proce¬ 
dure are being explored as a means of producing a 
semi-integrated subharmonically pumped two- 
diode mixer for use in 640 GHz. For this applica¬ 
tion, a planar antiparallel diode pair is incorporated 
into the center of a pair of low pass microstrip 
filters formed on an insulating GaAs substrate. An 
appealing substitute for millimeter and 
submillimeter wavelength applications is fused 
quartz which has an order of magnitude lower loss 
tangent, one-third the dielectric constant, and is 
much more mechanically rigid than GaAs. Using 
known techniques, it is possible to thin an entire 
GaAs wafer to a thickness on the order of 1-3 
microns and then transfer the remaining active 
layer onto quartz. An interesting question that 
arises as a result of this procedure is the effect of 
the ultrathin GaAs layer on the filter characteristics. 
A three-dimensional finite difference time domain 
(FO-TD) method is used to calculate the response 
of the filter both with and without the GaAs layer.* 
It is shown that even a very thin layer of high 
dielectric constant material can have a noticeable 
effect on the filter response. The FO-TO analysis 


agreed fairly well with measured response and 
could be used for subsequent design. 

In the last few years there has been an increasing 
interest in the interaction of electromagnetic fields 
with chiral media. A chiral medium is a subset of 
the wider class referred to as bianisotropic. Such 
media are characterized by linear constitutive 
relations that couple the electric and magnetic field 
vectors by three independent scalars. A rigorous 
formulation deals with the problem of radiation of 
electric and magnetic sources embedded in a 
layered lossless chiral medium. The fields are 
obtained in terms of dyadic Green's functions of 
electric or magnetic type represented in two- 
dimensional spectral-domain form. First, the 
spectral domain dyadL Green’s functions of elec¬ 
tric and magnetic types for an unbounded chiral 
medium are derived. The singular behavior of the 
various dyadic Green's functions in the source 
region is investigated and taken into account by 
extracting the delta function singularities. Finally, 
introducing global upward and downward reflec¬ 
tion and transmission matrices, the dyadic Green's 
functions in any layer of the stratified chiral 
medium are derived. 

The growth of epitaxial Fe film on GaAs was first 
demonstrated using molecular beam epitaxy. 
Deposition of such films by ion-beam sputtering 
was reported and their application to microwave 
filters has previously been discussed. A theoretical 
analysis is performed to the analysis of a tunable 
band-stop filter using epitaxial Fe films on [011] 
GaAs, in which the Fe layer is part of a microstrip 
line that runs along either a [100] or a [110] 
direction, the easy and hard directions of 
magnetization of the Fe film, respectively.* The 
analysis indicates that peak attenuation should 
occur at the ferromagnetic resonance (FMR) fre¬ 
quency and be proportional to the length of the 
microstrip line, and inversely proportional to the 
substrate thickness. 

A rigorous approach to the problem of radiation of 
electric or magnetic sources in a stratified arbi¬ 
trarily magnetized linear plasma.''* The fields are 
obtained in terms of dyadic Green’s functions of 
electric or magnetic type represented in the spec- 


* P.H. Siegel, J.E. Oswald, R.G. Dengler, O.M. Sheen, and S.M. Ali, 'Measured and Computed Performance of a 
Microstrip Filter Composed of Semi-Insulating GaAs on a Fused Quartz Substrate,' i£££ Microwave Guided Wave 
te«. 1 (4): April (1991). 
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tral domain. First, the dyadic Green's functions for 
an unbounded arbitrarily magnetized plasma is 
derived. The formulation is considerably simplified 
by using the KOB system of coordinates in con¬ 
junction with the Fourier transform. This leads to 
compact and explicit expressions for the dyadic 
Green's functions. The distributional singular 
behavior of the various dyadic Green's functions in 
the source region is investigated and taken into 
account by extracting the delta function singular¬ 
ities. Finally, the dyadic Green's function in any 
arbitrary layer is obtained in terms of appropriately 
defined global upward and downward reflection 
and transmission matrices. The field expressions 
for any arbitrary distribution of sources or linear 
antennas can be obtained by performing a convo¬ 
lution integral over the volume of the antenna 
weighted by the current density on the antenna. 

The input impedance of a microstrip antenna con¬ 
sisting of two circular microstrip disks in a stacked 
configuration driven by a coaxial probe is investi¬ 
gated.'^ A rigorous analysis is performed using a 
dyadic Green's function formulation where the 
mixed boundary value problem is reduced to a set 
of coupled vector integral equations using the 
vector Hankel transform. Galerkin's method is 
employed in the spectral domain where two sets of 
disk current expansions are used. One set is based 
on the complete set of orthogonal modes of the 
magnetic cavity, and the other employs Chebyshev 
polynomials with the proper edge condition for 
disk currents. An additional term is added to the 
disk current expansion to properly model the 
current in the vicinity of the probe/disk junction. 
The input impedance of the stacked microstrip 
antenna including the probe self-impedance is cal¬ 
culated as a function of the layered parameters and 
the ratio of the two disk radii. Disk current dis¬ 
tribution and radiation patterns are also presented. 
The calculated results are compared with exper¬ 
imental data and shown to be in good agreement. 

A method for the analysis of complex frequency- 
dependent signal interconnections and planar cir¬ 
cuits terminated with nonlinear load is 
investigated. The frequency-dependent portion of 
the system is first analyzed to get the scattering 
parameters and hence the impulse responses. Two 
different methods, depending on the complexity of 


the system, are used to get the scattering para¬ 
meters: the three-dimensional finite-difference 
time-domain method and an analytical procedure. 
The nonlinear convolution equations governing the 
overall system are then derived and solved numer¬ 
ically. The transient responses of a pair of coupled 
dispersive microstrip lines, a corner discontinuity, 
and a microstrip switching circuit are presented. 

We analysed three-dimensional (3-D) multichip 
module (MCM) interconnects, in this technology, 
the vertical interconnects consist of small con¬ 
ductor plated via holes etched by a photo lithog¬ 
raphy system. The via dimensions are in the same 
order as the microstrip and stripline widths to 
reduce the transmission line discontinuities. Two 
3-D transitions are investigated: microstrip-via- 
stripline and microstrip-via-90 degree stripline. The 
finite-difference time-domain with nonuniform grid 
is applied to the analysis of these transitions. Elec¬ 
tric field distribution and pulse propagation under 
the strip and the strip line are presented. The scat¬ 
tering parameters for various cases are calculated 
and compared. Geometrical effects such as dif¬ 
ferent conductor extensions on top of the vias and 
different hole sizes in the reference plane are also 
investigated. 

A method for calculation of the current distribu¬ 
tion, resistance, inductance matrices for a system 
of coupled superconducting transmission lines 
having finite rectangular cross section is devel- 
oped.'^ These calculations allow accurate charac¬ 
terization of both high-Te and low-Te 
superconducting transmission lines. For a single 
stripline geometry with finite ground planes, the 
current distribution, resistance, and kinetic induc¬ 
tance are calculated as a function of the pene¬ 
tration depth for various film thicknesses. These 
calculations are then used to determine the pene¬ 
tration depth for Nb, NbN, and YBa 2 Cu 307 -» 
superconducting thin films from the measured 
temperature dependence of the resonant frequency 
of a stripline resonator. The calculations are also 
used to convert measured temperature dependence 
of the quality factor to the intrinsic surface resist¬ 
ance as a function of temperature for an Nb strip¬ 
line resonator. 

The calculations of resistance and inductance are 
used for the determination of surface impedance of 
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YBa 2 Cu 307 -x thin films using a stripline 
resonator.''^ The surface impedance as a function 
of frequency from 1.5 to 20 GHz, as function of 
temperature from 4K to the transition temperature 
(90 K), and as a function of the RF magnetic field 
from zero to 300 Oe, is obtained. At low temper¬ 
atures the surface resistance of the films shows a 
very week dependance on the magnetic field up to 
225 to 250 Oe. At 77 K, the surface resistance is 
proportional to the square of the field. The pene¬ 
tration depth shows a much weaker dependence 
on the field than does the surface resistance. 

The Ginzburg-Landau (GL) theory is used to 
predict the nonlinear behavior in a supercon¬ 
ducting stripline resonator as a function of input. 
A method for calculating the nonlinear inductance 
and the fractional change in the resonant fre¬ 
quency of a stripline resonator is presented. By 
solving the GL equations inside the supercon¬ 
ducting strip, the spatial variation of the number 
density of superconducting electrons, and hence 
the spatial variation of the magnetic penetration 
depth are obtained for different values of input 
current. First, an infinite parallel plate transmission 
line is considered where the one-dimensional GL 
equations are solved. The two-dimensional case 
of strjpline is then considered. Nonlinear induc¬ 
tances are calculated as functions of input current 
for different superconducting strip lines. Compar¬ 
isons of the calculated resonance frequency shift 
with measurements for YBajCujO?.* stripline reso¬ 
nators show excellent agreement. 

A full-wave numerical analysis is applied to accu¬ 
rately characterize superconducting transmission 
lines embedded in a layered media. An integral 
equation formulation is developed by using a 
spectral domain Green's function for stratified 
media. Galerkin's method with roof top basis func¬ 
tions for the electric field distribution inside the 
superconductor is then employed to solve for the 
complex propagation constant. The thickness of 
the superconducting film is arbitrary in this anal¬ 
ysis, and the formulation rigorously accounts for 
the anisotropy of the superconducting film. The 
propagation characteristics of a superconducting 
microstrip transmission line with a thin dielectric 
buffer layer are investigated. A superconducting 
stripline configuration with an air gap is also 
studied. 

A full-wave analysis is applied to accurately char¬ 
acterize normal and superconducting coplanar strip 
lines (CPS) and coplanar waveguides (CPW). A 


volume integral equation formulation which rigor¬ 
ously considers the finite thickness, conductivity, 
and anisotropy of the strips, is developed. The 
phase and attenuation constants for copper CPS 
and CPW at 77 K and room temperatures are cal¬ 
culated and compared to experimental results 
obtained by short-pulse propagation measure¬ 
ments. 
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A branching model is proposed for the remote 
sensing of vegetation. The frequency and angular 
response of a two-scale cylinder cluster are calcu¬ 
lated to demonstrate the significance of vegetation 
architecture. The results indicate that it is neces¬ 
sary for theoretical remote sensing models to take 
into account the architecture of vegetation which 
plays an important role in determining the 
observed coherent effects. A two-scale branching 
model is implemented for soybeans with its 
internal structure and the resulting clustering 
effects considered. Furthermore, at the scale of 
soybean fields, the relative locations of soybean 
plants is described by a pair distribution function. 
The polarimetric backscattering coefficients are 
obtained in terms of the scattering properties of 
soybean plants and their pair distribution. Theore¬ 
tical backscattering coefficients are evaluated 
using a hole-corrected pair distribution function. 
Backscattering coefficients calculated are in good 
agreement with extensive data collected from 
soybean fields. It is found that the hole-correction 
approximation, which prevents two soybean plants 
from overlaying each other, greatly improves the 
agreement between the model and these three 
data sets near normal incidence, by introducing 
destructive interference at small angles of inci¬ 
dence. 
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